Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 6.0% intergenic (+28/‑488) sseA → / → ryfA 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3670GA 6.0% 300.6 / 5.1 99intergenic (+28/‑488)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (46/47);  new base A (4/2);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 101 total aligned reads displayed.

AACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTT  >  NC_000913/2653219‑2653512
                                                                                                                                                    |                                                                                                                                                 
aaCGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTCACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGa                                                                                                                                                  >  2:205609/1‑150 (MQ=255)
 aCGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGGATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGaa                                                                                                                                                 >  1:1142571/1‑150 (MQ=255)
      aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                                                                            <  2:368423/150‑1 (MQ=255)
       gCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGACAAACTGTACGACGGCGCATGGAGTGAATCGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAATTAAGAATTTTAAACTCAATAAAAAACCGCCGAATTTGGc                                                                                                                                           >  1:1669367/1‑150 (MQ=255)
        ccGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTAGGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACCGGCAATAAAAAACCGCCGAATTAGGCg                                                                                                                                          >  1:1222593/1‑150 (MQ=255)
        ccGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCg                                                                                                                                          >  2:780429/1‑150 (MQ=255)
         cGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                                                                         <  2:1619085/150‑1 (MQ=255)
          gTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                                                                        >  2:716423/1‑150 (MQ=255)
              ttttGTTAGCACTCGCGTCGCTGGTTTTGCCAATCTTGTATCTGTTCGACGGCGCATGTAGTGAATGGGGCTCGCGGGCTGTTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCTATAAAAAACCGGCGAATTTGGCGGttttt                                                                                                                                    <  2:475791/150‑1 (MQ=255)
               tttGTTCGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGTGCGCGCGGGCAGATTTACCGGTTGATCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttttt                                                                                                                                   <  1:827679/150‑1 (MQ=255)
               tttGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                                                                   <  1:370399/150‑1 (MQ=255)
                 tGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAt                                                                                                                                 <  2:515184/150‑1 (MQ=255)
                   ttAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATTGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                                                               <  2:1300402/150‑1 (MQ=255)
                   taaGCACTCGCGACGCTGGATGTGCCAATCGTGAAACTGTTCGTCTGAGCTTGGTGTGAATGGGGCGCGCGGGCTGATTTACCGGTTGTGCCAGTGTAATAAGTATTTTACAGGCAATAATAATCCGCCGAATTTGGCGGTTTTTTATTg                                                                                                                               <  2:1669367/148‑1 (MQ=255)
                     aGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                                                                             <  2:672336/150‑1 (MQ=255)
                     aGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGcc                                                                                                                             <  1:780429/150‑2 (MQ=255)
                     aGCACTCGCGACGCTGGATGCGCCAAACGTGACACGGTACGACGGCGCATGGAGTGAATGGGCCGCGCGGGCAGATGTACCGGTGGAGCCAGTGAAATAAGTATTTTACAGGCAATAAACAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                                                                             <  1:1514967/150‑1 (MQ=255)
                         cTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTc                                                                                                                         >  1:23992/1‑150 (MQ=255)
                           cgcgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTg                                                                                                                       <  1:1685339/150‑1 (MQ=255)
                             cgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGt                                                                                                                     >  2:1433401/1‑150 (MQ=255)
                                    ggATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGc                                                                                                              >  2:36563/1‑150 (MQ=255)
                                      ttgttCCAAACGGGAAACTGTACGACGGCGCCTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTTTTTTACAGGCAATCAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                                                                            <  2:121757/149‑1 (MQ=255)
                                           ccAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                                                                       >  2:838087/1‑150 (MQ=255)
                                            cAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                                                                      >  2:720294/1‑150 (MQ=255)
                                            cAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                                                                      >  1:192860/1‑150 (MQ=255)
                                               aCGTGACACTGTACGACGGCGCATGGAGTTCATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAAGTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagcag                                                                                                   <  2:1341044/150‑1 (MQ=255)
                                                 gTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                                                                 <  1:594709/150‑1 (MQ=255)
                                                 gTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                                                                 <  1:595277/150‑1 (MQ=255)
                                                   gAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTg                                                                                               <  2:1494243/150‑1 (MQ=255)
                                                    aaaCTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTATTTTATTGCTAGTCTGGTTCGCGGCCTATCCAGCAGGttca                                                                                              >  2:1453444/1‑148 (MQ=255)
                                                       cTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACtt                                                                                           >  2:1169889/1‑150 (MQ=255)
                                                         gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                         <  2:1622909/150‑1 (MQ=255)
                                                         gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                         >  1:1136497/1‑150 (MQ=255)
                                                         gTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                         >  2:1541922/1‑150 (MQ=255)
                                                         gTACGACGGCGCATGGAGTGAATGGGGCGAGCGGGCAGATTTACCGGTTGAGACAGTGAAATAAGTATTTTACAGGCAATAAAAAAACGCCGAATTTGGCGGTTTTTTATTTCTAGTCTTGTTCGAGGCCTTTCCAGCAGGTTGACTtgt                                                                                         >  2:825184/1‑150 (MQ=255)
                                                           tagacgGCGGATGGAGTGTATGGGTCTCGCGTTAAGATTTAACGGTTGAGCAATTGAAATAAGTATTTTACAGTCAATAAAAAACCGCCGATTTTGGGGGTTTTTTATTGATAGTTTGGTTGGCGGCCTTTCCAGCAGGTTGACTTgtgt                                                                                       <  2:144813/148‑1 (MQ=255)
                                                             gacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTa                                                                                     <  1:1541834/150‑1 (MQ=255)
                                                               cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                                                                   <  2:763579/150‑1 (MQ=255)
                                                                ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                                                                  <  1:534967/150‑1 (MQ=255)
                                                                  cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGa                                                                                <  1:962668/150‑1 (MQ=255)
                                                                    cATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                                                              >  2:71836/1‑150 (MQ=255)
                                                                        gAGTGAATGGGGCGCGCGGGCAGATTTACCGATTCAGCCTGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTCTTGCTAGTCTGGTTCGCGGCCTGTCCAGCAGGAGGACTTGTGTAACATGAGCAACg                                                                          >  1:260220/1‑150 (MQ=255)
                                                                         aGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                                                                         <  1:1272649/150‑1 (MQ=255)
                                                                            gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                                                                      >  2:160595/1‑150 (MQ=255)
                                                                            gAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                                                                      >  1:1065365/1‑150 (MQ=255)
                                                                             aaTGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGt                                                                     <  1:644462/150‑1 (MQ=255)
                                                                                ggggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                                                                  <  2:776912/150‑1 (MQ=255)
                                                                                ggggAGCGAGTGCAGATTTACCGGTTTAGCCAGTTAAATAAGTTTTTTACAGGCAATAAATAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                                                                  <  2:843335/150‑1 (MQ=38)
                                                                                 gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                                                                 <  1:674117/150‑1 (MQ=255)
                                                                                   gcgcgcGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                                                               <  2:599003/150‑1 (MQ=255)
                                                                                   gcgcgcGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTca                                                               <  2:1419294/150‑1 (MQ=255)
                                                                                    cgcgcgGGCAGATTTACCGGGTGAGCCAGTGAAATAAGTAATTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcac                                                              <  1:1459238/150‑1 (MQ=255)
                                                                                     gcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                                                             <  1:838087/150‑1 (MQ=255)
                                                                                        cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa                                                          >  1:606081/1‑150 (MQ=255)
                                                                                         gggCAGATTTACCTGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa                                                         <  1:952285/150‑1 (MQ=255)
                                                                                         gcgCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCATTCCACCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaa                                                         <  1:254272/148‑1 (MQ=255)
                                                                                            cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc                                                      <  2:1403772/150‑1 (MQ=255)
                                                                                                tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa                                                  >  1:1404361/1‑150 (MQ=255)
                                                                                                  tACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACt                                                >  2:1681190/1‑150 (MQ=255)
                                                                                                  tACCGGTTGAGACAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTACTCTGGTTCGCGGCCTTTCCAGCAGGTTGACATGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACt                                                >  2:273223/1‑150 (MQ=255)
                                                                                                      ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                            <  2:1219589/150‑1 (MQ=255)
                                                                                                      ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                            >  2:514232/1‑150 (MQ=255)
                                                                                                      ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                            <  2:192860/150‑1 (MQ=255)
                                                                                                      ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTCTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAcc                                            >  2:161697/1‑150 (MQ=255)
                                                                                                          gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCTAGGCCTTTCCAGCAGGTTGACTTGTGATACAAGAGCAACGCAGGAGCTTAACAGCAAAACAATACTCACAAGTa                                        >  1:1104351/1‑150 (MQ=255)
                                                                                                           aGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTaa                                       >  2:442091/1‑150 (MQ=255)
                                                                                                            gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAAc                                      >  1:1047207/1‑150 (MQ=255)
                                                                                                             ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAAct                                     >  2:470395/1‑150 (MQ=255)
                                                                                                             ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAAct                                     >  2:1275257/1‑150 (MQ=255)
                                                                                                               aGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAActca                                   >  2:1443968/1‑149 (MQ=255)
                                                                                                                  gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtt                                >  1:228626/1‑150 (MQ=255)
                                                                                                                    aaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt                              <  2:1583984/150‑1 (MQ=255)
                                                                                                                      tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGCTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCCTTTTg                            >  1:767801/1‑150 (MQ=255)
                                                                                                                        attATTTTACAGGCAATAAAAAACCGCCGAATTTTGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTc                          <  2:68646/148‑1 (MQ=255)
                                                                                                                        aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTc                          <  2:23992/150‑1 (MQ=255)
                                                                                                                        aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTc                          >  1:991646/1‑150 (MQ=255)
                                                                                                                         gTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCa                         <  1:720294/150‑1 (MQ=255)
                                                                                                                          tATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCATTCAAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCaa                        >  1:476489/1‑150 (MQ=255)
                                                                                                                           aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGTAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAg                       <  1:487733/150‑1 (MQ=255)
                                                                                                                           aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAg                       <  1:20954/150‑1 (MQ=255)
                                                                                                                            ttttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGc                      <  1:1541922/150‑1 (MQ=255)
                                                                                                                               ttcAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaa                   <  1:1054317/148‑1 (MQ=255)
                                                                                                                               tACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTGATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGCCTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaa                   >  1:1264758/1‑150 (MQ=255)
                                                                                                                                 cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACTTGAGCAAGGCAGGTGCGTCGCAGCAAAAAAATTCTAAGCAGTATCTATCTTTTTGGCAAGAAAAag                 >  1:172629/1‑150 (MQ=37)
                                                                                                                                 cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAag                 >  2:5826/1‑150 (MQ=255)
                                                                                                                                 cAGGCAATAAAAAACCGCCGAAGTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaatg                 >  1:71361/1‑148 (MQ=255)
                                                                                                                                   ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag               >  2:361415/1‑150 (MQ=255)
                                                                                                                                   ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag               >  2:257901/1‑150 (MQ=255)
                                                                                                                                     cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagagag             <  2:1648757/150‑1 (MQ=255)
                                                                                                                                     cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagagag             <  2:29383/150‑1 (MQ=255)
                                                                                                                                     cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagagag             >  1:1372832/1‑150 (MQ=255)
                                                                                                                                       aTAAAAAACCGCCAAATTCGGCGGTTGTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGtt           >  1:275705/1‑149 (MQ=255)
                                                                                                                                         aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAat         >  2:712809/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAat         >  1:1576595/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCGCAAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCATTCCAGCAGGTTGAATTGTGTTACATGAGCAAAGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTATTTTTGTAAATCAAAATAGAGTAat         >  2:285285/1‑150 (MQ=255)
                                                                                                                                            aaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAattat      >  1:1594245/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTg    <  2:407541/150‑1 (MQ=255)
                                                                                                                                              aCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGATGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTg    >  2:113995/1‑150 (MQ=255)
                                                                                                                                                cgccgAATTTGGCTGTTTTTTATTGCGAGTCTGGTTCGCGGCCTTTCCAGCAGGTTTAATTGTGTTACATTAGCAACTCAGGTGCTTCACAGCAAACCAATACTCACCAGTAAGACTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtt  <  2:838951/150‑1 (MQ=37)
                                                                                                                                                cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtt  <  1:596033/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTT  >  NC_000913/2653219‑2653512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑