Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 231,287 A→G M781V (ATG→GTG)  SAUSA300_RS01035 → type‑1 restriction enzyme R protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793231,2870AG100.0% 376.6 / NA 129M781V (ATG→GTG) SAUSA300_RS01035type‑1 restriction enzyme R protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (122/7);  total (122/7)

Only 100 of 132 total aligned reads displayed.

TAGCTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGATGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACAAA  >  NC_007793/231201‑231383
                                                                                      |                                                                                                
tAGCTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATAc                                                                                     >  1:3191409/1‑100 (MQ=255)
tAGCTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATAc                                                                                     >  1:4434211/1‑100 (MQ=255)
agccTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATAc                                                                                     >  1:5123818/4‑100 (MQ=255)
  gCTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGa                                                                                   >  1:1712648/1‑100 (MQ=255)
  ccTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGa                                                                                   >  1:3629442/2‑100 (MQ=255)
    tgtgTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGAtt                                                                                 >  1:3579134/1‑100 (MQ=255)
     gtgtACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTa                                                                                >  1:5631761/1‑100 (MQ=255)
     gtgtACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTa                                                                                >  1:1811965/1‑100 (MQ=255)
       gtACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAAt                                                                              >  1:3385126/1‑100 (MQ=255)
       gtACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAAt                                                                              >  1:2538632/1‑100 (MQ=255)
            ttCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGaa                                                                         >  1:1105855/2‑100 (MQ=255)
            aTCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGag                                                                         <  1:560865/100‑2 (MQ=255)
            aTCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGag                                                                         <  1:1822395/100‑2 (MQ=255)
             acAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATAAATGTGAat                                                                        >  1:1164921/2‑100 (MQ=255)
              cAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAatt                                                                       >  1:755922/1‑100 (MQ=255)
              cAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAatt                                                                       >  1:502478/1‑100 (MQ=255)
              cAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAatt                                                                       >  1:4148482/1‑100 (MQ=255)
                ctgtgctagagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattat                                                                     >  1:805825/8‑100 (MQ=255)
                aGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattat                                                                     >  1:391883/1‑100 (MQ=255)
                aGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattat                                                                     >  1:4444770/1‑100 (MQ=255)
                  aaaaaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatat                                                                   >  1:3820126/2‑100 (MQ=255)
                   aaaaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatt                                                                  >  1:925181/1‑100 (MQ=255)
                   aaaaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatt                                                                  >  1:3777688/1‑100 (MQ=255)
                    ttaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:4981181/3‑100 (MQ=255)
                    tattgagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:549065/5‑100 (MQ=255)
                    tatgcagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:1513398/6‑100 (MQ=255)
                    tagctagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:2135193/6‑100 (MQ=255)
                    tacagagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:5644225/4‑100 (MQ=255)
                    taaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:4688108/2‑100 (MQ=255)
                    ctgagagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:4419666/4‑100 (MQ=255)
                    ctcagagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:5075075/4‑100 (MQ=255)
                    ccaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:1019051/3‑100 (MQ=255)
                    aacgcagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:5768385/6‑100 (MQ=255)
                    aaaaGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatatta                                                                 >  1:3144913/1‑100 (MQ=255)
                     tattagCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAattatattat                                                                >  1:1917155/5‑100 (MQ=255)
                       ctggCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGa                                                              >  1:1847765/4‑100 (MQ=255)
                        ctgCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGaa                                                             >  1:3517694/3‑100 (MQ=255)
                        aggCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGaa                                                             >  1:1030808/3‑100 (MQ=255)
                         gtggACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATATGATTAATGTGAATTATATTATGAat                                                            >  1:5326449/4‑100 (MQ=255)
                         ggCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAat                                                            >  1:1661117/2‑100 (MQ=255)
                         agCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAat                                                            >  1:923317/1‑100 (MQ=255)
                         agCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAat                                                            >  1:2114739/1‑100 (MQ=255)
                         agCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAat                                                            >  1:3692896/1‑100 (MQ=255)
                          gCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAata                                                           >  1:135134/1‑100 (MQ=255)
                          gCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAata                                                           >  1:2038176/1‑100 (MQ=255)
                          ccGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAata                                                           >  1:110100/2‑100 (MQ=255)
                          cagACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAata                                                           >  1:1722299/3‑100 (MQ=255)
                           tgccGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATgcttat                                                          >  1:3593887/4‑95 (MQ=255)
                           gtgtggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:5080385/5‑100 (MQ=255)
                           ggACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:320896/2‑100 (MQ=255)
                           ctaCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:3226123/3‑100 (MQ=255)
                           caaCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:1995889/3‑100 (MQ=255)
                           cGACGGCTGGGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:702504/1‑100 (MQ=255)
                           cGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:2073093/1‑100 (MQ=255)
                           cGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:1517207/1‑100 (MQ=255)
                           acaCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatat                                                          >  1:2261038/3‑100 (MQ=255)
                            taCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:3041276/2‑100 (MQ=255)
                            taCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGCGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:2826989/2‑100 (MQ=255)
                            gccGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGGAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:242971/3‑100 (MQ=255)
                            gACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:2561586/1‑100 (MQ=255)
                            gACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:3966830/1‑100 (MQ=255)
                            gACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:2927464/1‑100 (MQ=255)
                            caCGGCTGAGAAAAATAAAGTATCCATTTTAAATGGTATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:2101336/2‑100 (MQ=255)
                            caCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAGATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:3359145/2‑100 (MQ=255)
                            caCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:3331781/2‑100 (MQ=255)
                            caCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAACAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:1808092/2‑100 (MQ=255)
                            aaCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatata                                                         >  1:504469/2‑100 (MQ=255)
                             gtggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:1171396/3‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGTTACGATTAATGTGAATTATATTATGAatatat                                                        >  1:3888461/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:5276064/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:5195920/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:3161782/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:4326480/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:5497125/1‑100 (MQ=255)
                             aCGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAatatat                                                        >  1:833060/1‑100 (MQ=255)
                              cGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATAtt                                                       >  1:718660/1‑100 (MQ=255)
                              cGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATAtt                                                       >  1:2800052/1‑100 (MQ=255)
                              cGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATAtt                                                       >  1:3916101/1‑100 (MQ=255)
                               ggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTg                                                      >  1:2443347/1‑100 (MQ=255)
                               ggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTg                                                      <  1:5026842/100‑1 (MQ=255)
                               ggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTg                                                      <  1:4237505/100‑1 (MQ=255)
                               ggCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTg                                                      <  1:784075/100‑1 (MQ=255)
                                   gagaAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGAc                                                  >  1:1838375/1‑100 (MQ=255)
                                                      ttctAAATGATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACaa                               <  1:3248561/97‑1 (MQ=255)
                                                              gcccttGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGaac                       >  1:4769120/5‑100 (MQ=255)
                                                              gATATTGATTTCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGaac                       >  1:1491654/1‑100 (MQ=255)
                                                                    actgtCGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAAcgtc                 >  1:152203/5‑100 (MQ=255)
                                                                         cGAAATAGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAAc            >  1:2715342/1‑100 (MQ=255)
                                                                               aGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGaa      >  1:586107/1‑100 (MQ=255)
                                                                               aGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGaa      >  1:1691964/1‑100 (MQ=255)
                                                                               aGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGaa      >  1:1261372/1‑100 (MQ=255)
                                                                               aGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGaa      >  1:5591045/1‑100 (MQ=255)
                                                                               aGAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGaa      >  1:4671021/1‑100 (MQ=255)
                                                                                gAAATGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAAc     >  1:2390747/1‑100 (MQ=255)
                                                                                aaaaTGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAAc     >  1:5638121/2‑100 (MQ=255)
                                                                                 aaaTGGTGCGTAGTGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACa    >  1:2924948/1‑100 (MQ=255)
                                                                                  aaTGGTGCGTAATGATACGATTAATGTGAATTATATTATGACTATATTGAGACAAATTGATCTTGAAGACAAAGCGCAACAACGTCGTAACCAAGAACaa   >  1:5380072/1‑100 (MQ=255)
                                                                                  aaTGGTGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACaa   >  1:4110094/1‑100 (MQ=255)
                                                                                   atggtgCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACaaa  >  1:4052935/1‑100 (MQ=255)
                                                                                   atggtgCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACaaa  >  1:25509/1‑100 (MQ=255)
                                                                                      |                                                                                                
TAGCTGTGTACGATCAAGTAAAAAGAGCGACGGCTGAGAAAAATAAAGTATCCATTTTAAATGATATTGATTTCGAAATAGAAATGATGCGTAATGATACGATTAATGTGAATTATATTATGAATATATTGAGACAAATTGATCTTGAAGACAAAGCGGAACAACGTCGTAACCAAGAACAAA  >  NC_007793/231201‑231383

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑