Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,5140CT21.1% 127.2 / 11.3 57E196K (GAA→AAA) cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (16/29);  new base T (12/0);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.30e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGA  >  NC_000913/2889274‑2889764
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
aGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGc                                                                                                                                                                                                                                                  <  1:7793/251‑1 (MQ=255)
         cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTcc                                                                                                                                                                                                                                         <  2:152491/251‑1 (MQ=255)
         cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCGGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTcc                                                                                                                                                                                                                                         <  1:431753/251‑1 (MQ=255)
                   gCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGGAGGGCTGGGCGAGGATCGGTTCTTTATCAGTAGTGGGGGCTTTTTTCTGGTCGAGCCAGATCTCCGCATAcgcgc                                                                                                                                                                                                                               >  1:363634/1‑251 (MQ=255)
                     gggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTGATTTTCGCGGCAGGGAGGTCTGGCCGGGGGTCGGTTCTTCTTCAGTAGTGGCGACTTTTTCCTTGTCGAGGCAGATCTCCGccttagcccgg                                                                                                                                                                                                                             >  2:159418/1‑241 (MQ=255)
                                      ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTGTATTTTAGCGGCAGGGAGGGCTGGGCGCGGATTCGGTTTTTATTTGGTGTGGGGCGCTTTTTCTGGTCGCGCCAGATCTCCGCCTACCCGCGGGGACGGGGCAACAGa                                                                                                                                                                                                              >  2:321667/1‑249 (MQ=255)
                                             gtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGa                                                                                                                                                                                                      <  2:469836/250‑1 (MQ=255)
                                                    ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtcttc                                                                                                                                                                                               <  1:397397/250‑1 (MQ=255)
                                                      aaaGCCCCCCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCTTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCttcttc                                                                                                                                                                                            <  2:476801/251‑1 (MQ=255)
                                                                 cccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTCTCCTGGTCGATCCAGATCTCCGCATACGCGCGGGGTCGAGTCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTa                                                                                                                                                                                 <  2:461378/251‑1 (MQ=255)
                                                                          aGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATTCCTACCGCGGTTTTTAATTTACGCGGCAAGGAGGTCTGGGCGAGGGTCGGTTTTTCAACAGTAATGGCGGCTTTTTTCTGGGCGAGCCAAACCTCCGCATACGACCGGGTGCGCGGCAACAAATGCTCCGCAATCTTCTTCccccctccctacgcctccgc                                                                                                                                                                        >  1:69481/1‑235 (MQ=255)
                                                                            gTTAAAGCCCAACAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTTAATTTTCGCGGCAGGGAGGGCTGGCCGAGGGTCGGGTCTTTATCAGTAGTGGGGGCTTTTTTCTGGTCGGGCCCGATCTCCCCATACGCGCGGGGACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACCCGTACGCTTCcgccg                                                                                                                                                                      >  1:141078/1‑249 (MQ=255)
                                                                              tAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCGCTACCGCGGGTTTTAATTTACGCGGGAGGTAGGGCTGGCCGAGGAGCGGTTCTTCTTTCGGAGAGGCGCCTTTTTCCTGCGCTGGCCAGATTCCCCCATACGCCCGGGGACAAGGCAACAGATCCTCCGCAATCTTCTTTGCCCACCCGTACGCTTCCTCGTGc                                                                                                                                                                    >  1:315009/1‑251 (MQ=255)
                                                                               aaaGCCCCCCAGCTGGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTGCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCGTCGCCCACTCGTACGCTGCCGCGTGc                                                                                                                                                                    <  2:390678/250‑1 (MQ=255)
                                                                               aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGc                                                                                                                                                                    <  2:469291/250‑1 (MQ=255)
                                                                                    ccaccaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTCCCGGGGTTTTGAATTTTCGCGGCCGGGAGGTCTGGCCGAGGGTCGGGTCTTCTTCAGTAGTGGGGACTTTTTTCTGGGCGCGCCAGATATCCGCATAAGCGCGGGGACGAGGGCACAAATACTCAGAAATCTTTTTCTCCCACTCGTACCCTTCCGCGTGCAGCTGc                                                                                                                                                              >  1:352634/1‑251 (MQ=255)
                                                                                            ttGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGGAGGTCTGGCCGAGGATCGGTTCTTTATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGAGGCAACAAATGCTCAGAAATCTTCTTCGCCCCCTCGTACGCTTCCGCGGGCAGCTGCGACTCGt                                                                                                                                                       >  1:181528/1‑250 (MQ=255)
                                                                                             tGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGAGCCCAACCGGGGTTTTTAATTTTCGCGGCAGGGAGGGCTGGGCGGGGGTCGGGTTTTTCTCAGGAGTGGGGGCTTTTTTCTGGTGGAGCCCGATCTCCGCCTTCTCGCGGGGAGGAGGCAACACATGCTCAGAAATCTTTTTCTCCCCCCCGTACGCCTCCGCGCGCAGCGCCGACTCggtg                                                                                                                                                     >  2:59145/1‑248 (MQ=255)
                                                                                                 gTTTTCGGCGATCGCCCCGAAGTTCATGGCGTTGGCGTGCAGAGCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGt                                                                                                                                                 <  2:45398/251‑1 (MQ=255)
                                                                                                  ttgtCGGCGATCGCCACGAAGTGCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTTGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGTCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGt                                                                                                                                                 <  1:295196/247‑1 (MQ=255)
                                                                                                    ttCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATAGCGTTCTGGGGCGGGAGCACGACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGtta                                                                                                                                               <  2:134378/250‑2 (MQ=255)
                                                                                                      cGGCGATCGCCACGCAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGGAGGGCTGGCCGAGGATCGGGTCTTCATCAGTAGTGGCGGCTTTTTTCTGGGCGGGCCAGATCTCCCCATAAGCGCGGGGACGAGGCAACAGATGGTCAGAAATCTTCTTCGCCCACTCGTACGCCTCCCCGTGCAGCTGCGCCTCGGAAggggtgg                                                                                                                                              >  1:122337/1‑245 (MQ=255)
                                                                                                       ggCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCATTCTGTGGCGGGATCACTACCGCGGATTTTTATTTTCTCGGCAGGTAGGCTTGGCCGGGGGTCGGTTTTTCATCAGTAGTGGCGACTTTTTTCTGGGCGAGACAGATCTCCGCATACGCGCGGGTACCCGGGAACAAATTCTCATAAATCTTCCTTGCCCCCCCGTTCGCTTTCGCGGGGAGTTGCGACTCGTAAAGGTTCgcgg                                                                                                                                           >  1:473925/1‑248 (MQ=255)
                                                                                                        gCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGt                                                                                                                                          <  1:425548/251‑1 (MQ=255)
                                                                                                         cGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTg                                                                                                                                         <  1:159418/251‑1 (MQ=255)
                                                                                                          gATCGCCACGAAGTTCCTGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGGGGGTTTGACTTTTCGCGGCAGGGAGGTCTGGCCGAGGATCGGTTATTCATCAGGAGTGGCGACTTTTTTCTGGTCGCGCAAAATCTCCGCCTACGCGCGGGGACGAGGCAACAAAAGCTCAGAAAACTTCTTTGCCCACTCGTACGCTTCCCCGGGCAGCTGCGCCTCGTAAGGGGTCGGGGTgg                                                                                                                                        >  1:181939/1‑250 (MQ=255)
                                                                                                              gCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGCGGGTTTTAATTTACGCGGCAGGGAGGTCTGGCCGAGGATCGGTTTTTCATCAGTAGGGGCGGCTTTTTTCTGGGCGAGCCAGATCTCCGCATACGCGCGGGGAAGCGGCAACAGAAGCTCAGAAATCTTTTTCGCCCACTCCTTCGCTTCCGCGTGCCGCTGCGCCTCGTTAGGGGTCGAgggggaggg                                                                                                                                    >  1:378916/1‑244 (MQ=255)
                                                                                                                   gAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTGGCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTGCGCCCACTCGTACGCGTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGt                                                                                                                               <  2:187453/251‑1 (MQ=255)
                                                                                                                   gAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGTGGTTTTTGATTTACCCGGCAGGGAGGTCTGGCCGGGGATCGGGTCTTCCTCAGGAGTGGCGGCTTTTTTCCGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGCGGCAAAAAATGCTCAAAAATATTCTTCGCCCACTCGTACGCTTCCGCCGGCAGCTGCGGCCCGCAAGGGTTCGCGGGGCagcgttcgt                                                                                                                               >  1:410962/1‑244 (MQ=255)
                                                                                                                   gAAGTGCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGt                                                                                                                               <  2:472264/251‑1 (MQ=255)
                                                                                                                      attCATGTCGTTGGCCGGCCAGTCCGATATGTTGTGTGGCCGGATCCCTACCGTGGGTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTAc                                                                                                                            <  2:12908/250‑1 (MQ=255)
                                                                                                                                      gTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAATTTACGCGGGCGGTAGGTCTGGCCGAGGATCGGGTCTTCATCAGTAGGGGCGACTTTTTTCTGGGCGAGCCCGATCTCCGCATAAGCGCGGGGACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCCCTCGTACGCTTCCGCGTGCCGGTGCGGCTCGGAAGGGTTCGAGGGGCAGAGGAAGTTACGGTTCATGTCGtt                                                                                                               >  1:63396/1‑249 (MQ=255)
                                                                                                                                       ggCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGc                                                                                                            <  2:230004/250‑1 (MQ=255)
                                                                                                                                         cAGATCGATATCGTTCTGTGGCGGGATCCATACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGCGGGTCGGGTCTTCATCAGTAGTGGCGGCTTTTTTCTGGTCGAGGCCGATCTCCGCCTACGCGCGGGGACGCGGCAACCGATGCTCAGAAAACTTCTTCGCCCCCTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGGTCGGGGGGCAGAGGACGTTACGGGTCAttgcgtttgctg                                                                                                          >  1:387021/1‑240 (MQ=255)
                                                                                                                                           gATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTc                                                                                                        <  1:321667/251‑1 (MQ=255)
                                                                                                                                           gATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTc                                                                                                        <  1:420959/251‑1 (MQ=255)
                                                                                                                                            aTCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTc                                                                                                                                                                                                                                                      <  1:246460/108‑1 (MQ=255)
                                                                                                                                                atatCGTGCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTGAGCTGTCGCCAg                                                                                                   <  2:227605/251‑1 (MQ=255)
                                                                                                                                                                      cTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgc                                                                             >  1:312735/1‑251 (MQ=255)
                                                                                                                                                                        aCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGCCCGACCGCGTgc                                                                             >  1:18773/1‑249 (MQ=255)
                                                                                                                                                                               ttttGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTATTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTgg                                                                    >  1:211014/1‑251 (MQ=255)
                                                                                                                                                                                  tGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                                                                 >  1:183307/1‑251 (MQ=255)
                                                                                                                                                                                   gAATGTACGCGTCAGGTAGGTCGGTCCGCGGATCGGTTCTTCAGCAGTAGTGGCGACTTTCTCCTGTTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATTTCGGTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                                                                 <  2:68111/250‑1 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              >  2:435968/1‑251 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTAGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              >  1:457329/1‑251 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGGCTGGCCGAGGATCGGTTCTTCATCAGTAGGGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              <  2:298901/251‑1 (MQ=255)
                                                                                                                                                                                        tACGCGGCAGGTAGGCCTGGCCGAGGATCGGTTATTCTTCCGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGAACGTTACGGTTCATGTCGTTGGCTGTCGCCAGCGCATCAAGCCCGACCGGGTGCAGAACCTGGTGCCCCgtt                                                           >  2:485626/1‑249 (MQ=255)
                                                                                                                                                                                                aggtaggtCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                                                   >  1:125505/1‑251 (MQ=255)
                                                                                                                                                                                                aggtaggtCTGGCAGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                                                   >  2:228189/1‑251 (MQ=255)
                                                                                                                                                                                                  gtaggtCTGGCCGAGGATCGGTTCTTCAGCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGTGCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCt                                                 <  1:46158/251‑1 (MQ=255)
                                                                                                                                                                                                    aggtCTGGCCGAAGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                                                <  2:118994/250‑1 (MQ=255)
                                                                                                                                                                                                        cTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                                            >  1:114680/1‑250 (MQ=255)
                                                                                                                                                                                                                     cGGTTCTTCATCAGCAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg                                >  2:252590/1‑250 (MQ=255)
                                                                                                                                                                                                                             tcatcaGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                         <  1:460780/250‑1 (MQ=255)
                                                                                                                                                                                                                              catcaGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaaa                        <  2:363634/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aCTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGt         >  1:433425/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ttttCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTc        <  1:92666/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccTGGTCGAGAAAGATCTCCGAATACGCGCGGGTACGCGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGAGCAGCTGCGACTCGTAAGGGTTCGCGGTGCAGAGTACGTTACGGTTCATGTCGTTCGCTGCCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAAGTCTGGCGGTTGGTCAGGCgc  >  2:437450/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
AGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGA  >  NC_000913/2889274‑2889764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: