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.

AGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGA  >  NC_000913/2889274‑2889764
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
aGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGc                                                                                                                                                                                                                                                  <  1:7793/251‑1 (MQ=255)
         cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTcc                                                                                                                                                                                                                                         <  2:152491/251‑1 (MQ=255)
         cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCGGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTcc                                                                                                                                                                                                                                         <  1:431753/251‑1 (MQ=255)
                   gCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGGAGGGCTGGGCGAGGAT‑CGG‑TTCTTTATCAGTAGTGGGG‑GCTTTTTTCTGGTCGAGCCAGATCTCCGCATAcgcgc                                                                                                                                                                                                                               >  1:363634/1‑251 (MQ=255)
                     gggCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTGATTTTCGCGGCAGGGAGGTCTGGCCGGGGGT‑CGG‑TTCTTCTTCAGTAGTGGCG‑ACTTTTTCCTTGTCGAGGCAGATCTCCGccttagcccgg                                                                                                                                                                                                                             >  2:159418/1‑241 (MQ=255)
                                      ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTGTATTTTAGCGGCAGGGAGGGCTGGGCGCGGATTCGG‑TT‑TTTATTTGGTGTGGGGCGCTTTTTC‑TGGTCGCGCCAGATCTCCGCCTACCCGCGGGGACGGGGCAACAGa                                                                                                                                                                                                              >  2:321667/1‑249 (MQ=255)
                                             gtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGa                                                                                                                                                                                                      <  2:469836/250‑1 (MQ=255)
                                                    ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtcttc                                                                                                                                                                                               <  1:397397/250‑1 (MQ=255)
                                                      aaaGCCCCCCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCTTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCttcttc                                                                                                                                                                                            <  2:476801/251‑1 (MQ=255)
                                                                 cccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTCTCCTGGTCGATCCAGATCTCCGCATACGCGCGGGGTCGAGTCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTa                                                                                                                                                                                 <  2:461378/251‑1 (MQ=255)
                                                                          aGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATTCCTACCGCGGTTTTTAATTTACGCGGCAAGGAGGTCTGGGCGAGGGT‑CGG‑TTTTTCAACAGTAATGGCG‑GCTTTTTTCTGGGCGAGCCAAACCTCCGCATACGACCGGGTGCGCGGCAACAAATGCTCCGCAATCTTCTTCccccctccctacgcctccgc                                                                                                                                                                        >  1:69481/1‑235 (MQ=255)
                                                                            gTTAAAGCCCAACAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTTAATTTTCGCGGCAGGGAGGGCTGGCCGAGGGT‑CGG‑GTCTTTATCAGTAGTGGGG‑GCTTTTTTCTGGTCGGGCCCGATCTCCCCATACGCGCGGGGACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACCCGTACGCTTCcgccg                                                                                                                                                                      >  1:141078/1‑249 (MQ=255)
                                                                              tAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCGCTACCGCGGGTTTTAATTTACGCGGGAGGTAGGGCTGGCCGAGGAG‑CGG‑TTCTTCTTTCGGAGAGGCG‑CCTTTTTCCTGCGCTGGCCAGATTCCCCCATACGCCCGGGGACAAGGCAACAGATCCTCCGCAATCTTCTTTGCCCACCCGTACGCTTCCTCGTGc                                                                                                                                                                    >  1:315009/1‑251 (MQ=255)
                                                                               aaaGCCCCCCAGCTGGCCGTTTGCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTGCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCGTCGCCCACTCGTACGCTGCCGCGTGc                                                                                                                                                                    <  2:390678/250‑1 (MQ=255)
                                                                               aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGc                                                                                                                                                                    <  2:469291/250‑1 (MQ=255)
                                                                                    ccaccaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTCCCGGGGTTTTGAATTTTCGCGGCCGGGAGGTCTGGCCGAGGGT‑CGG‑GTCTTCTTCAGTAGTGGGG‑ACTTTTTTCTGGGCGCGCCAGATATCCGCATAAGCGCGGGGACGAGGGCACAAATACTCAGAAATCTTTTTCTCCCACTCGTACCCTTCCGCGTGCAGCTGc                                                                                                                                                              >  1:352634/1‑251 (MQ=255)
                                                                                            ttGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGGAGGTCTGGCCGAGGAT‑CGG‑TTCTTTATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGAGGCAACAAATGCTCAGAAATCTTCTTCGCCCCCTCGTACGCTTCCGCGGGCAGCTGCGACTCGt                                                                                                                                                       >  1:181528/1‑250 (MQ=255)
                                                                                             tGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGAGCCCAACCGGGGTTTTTAATTTTCGCGGCAGGGAGGGCTGGGCGGGGGT‑CGGGTTTTTC‑TCAGGAGTGGGG‑GCTTTTTTCTGGTGGAGCCCGATCTCCGCCTTCTCGCGGGGAGGAGGCAACACATGCTCAGAAATCTTTTTCTCCCCCCCGTACGCCTCCGCGCGCAGCGCCGACTCggtg                                                                                                                                                     >  2:59145/1‑248 (MQ=255)
                                                                                                 gTTTTCGGCGATCGCCCCGAAGTTCATGGCGTTGGC‑GTGCAGAGCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGt                                                                                                                                                 <  2:45398/251‑1 (MQ=255)
                                                                                                  ttgtCGGCGATCGCCACGAAGTGCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTTGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGTCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGt                                                                                                                                                 <  1:295196/247‑1 (MQ=255)
                                                                                                    ttCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATAGCGTTCTGGGGCGGGAGCACGACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGtta                                                                                                                                               <  2:134378/250‑2 (MQ=255)
                                                                                                      cGGCGATCGCCACGCAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGGAGGGCTGGCCGAGGAT‑CGG‑GTCTTCATCAGTAGTGGCG‑GCTTTTTTCTGGGCGGGCCAGATCTCCCCATAAGCGCGGGGACGAGGCAACAGATGGTCAGAAATCTTCTTCGCCCACTCGTACGCCTCCCCGTGCAGCTGCGCCTCGGAAggggtgg                                                                                                                                              >  1:122337/1‑245 (MQ=255)
                                                                                                       ggCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCATTCTGTGGCGGGATCACTACCGCGGATTTTTATTTTCTCGGCAGGTAGGCTTGGCCGGGGGT‑CGG‑TTTTTCATCAGTAGTGGCG‑ACTTTTTTCTGGGCGAGACAGATCTCCGCATACGCGCGGGTACCCGGGAACAAATTCTCATAAATCTTCCTTGCCCCCCCGTTCGCTTTCGCGGGGAGTTGCGACTCGTAAAGGTTCgcgg                                                                                                                                           >  1:473925/1‑248 (MQ=255)
                                                                                                        gCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGt                                                                                                                                          <  1:425548/251‑1 (MQ=255)
                                                                                                         cGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTg                                                                                                                                         <  1:159418/251‑1 (MQ=255)
                                                                                                          gATCGCCACGAAGTTCCTGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGGGGGTTTGACTTTTCGCGGCAGGGAGGTCTGGCCGAGGAT‑CGG‑TTATTCATCAGGAGTGGCG‑ACTTTTTTCTGGTCGCGCAAAATCTCCGCCTACGCGCGGGGACGAGGCAACAAAAGCTCAGAAAACTTCTTTGCCCACTCGTACGCTTCCCCGGGCAGCTGCGCCTCGTAAGGGGTCGGGGTgg                                                                                                                                        >  1:181939/1‑250 (MQ=255)
                                                                                                              gCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGCGGGTTTTAATTTACGCGGCAGGGAGGTCTGGCCGAGGAT‑CGG‑TTTTTCATCAGTAGGGGCG‑GCTTTTTTCTGGGCGAGCCAGATCTCCGCATACGCGCGGGGAAGCGGCAACAGAAGCTCAGAAATCTTTTTCGCCCACTCCTTCGCTTCCGCGTGCCGCTGCGCCTCGTTAGGGGTCGAgggggaggg                                                                                                                                    >  1:378916/1‑244 (MQ=255)
                                                                                                                   gAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTGGCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTGCGCCCACTCGTACGCGTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGt                                                                                                                               <  2:187453/251‑1 (MQ=255)
                                                                                                                   gAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGTGGTTTTTGATTTACCCGGCAGGGAGGTCTGGCCGGGGAT‑CGG‑GTCTTCCTCAGGAGTGGCG‑GCTTTTTTCCGGTCGAGCCAGATCTCCGCATACGCGCGGGGACGCGGCAAAAAATGCTCAAAAATATTCTTCGCCCACTCGTACGCTTCCGCCGGCAGCTGCGGCCCGCAAGGGTTCGCGGGGCagcgttcgt                                                                                                                               >  1:410962/1‑244 (MQ=255)
                                                                                                                   gAAGTGCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGt                                                                                                                               <  2:472264/251‑1 (MQ=255)
                                                                                                                      attCATGTCGTTGGCCGGCCAGTCCGATAT‑GTTGTGTGGCCGGATCCCTACCGTGGGTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTAc                                                                                                                            <  2:12908/250‑1 (MQ=255)
                                                                                                                                      gTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTGAATTTACGCGGGCGGTAGGTCTGGCCGAGGAT‑CGG‑GTCTTCATCAGTAGGGGCG‑ACTTTTTTCTGGGCGAGCCCGATCTCCGCATAAGCGCGGGGACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCCCTCGTACGCTTCCGCGTGCCGGTGCGGCTCGGAAGGGTTCGAGGGGCAGAGGAAGTTACGGTTCATGTCGtt                                                                                                               >  1:63396/1‑249 (MQ=255)
                                                                                                                                       ggCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGc                                                                                                            <  2:230004/250‑1 (MQ=255)
                                                                                                                                         cAGATCGATATCGTTCTGTGGCGGGATCCATACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGCGGGT‑CGG‑GTCTTCATCAGTAGTGGCG‑GCTTTTTTCTGGTCGAGGCCGATCTCCGCCTACGCGCGGGGACGCGGCAACCGATGCTCAGAAAACTTCTTCGCCCCCTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGGTCGGGGGGCAGAGGACGTTACGGGTCAttgcgtttgctg                                                                                                          >  1:387021/1‑240 (MQ=255)
                                                                                                                                           gATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTc                                                                                                        <  1:321667/251‑1 (MQ=255)
                                                                                                                                           gATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTc                                                                                                        <  1:420959/251‑1 (MQ=255)
                                                                                                                                            aTCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTc                                                                                                                                                                                                                                                      <  1:246460/108‑1 (MQ=255)
                                                                                                                                                atatCGTGCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTGAGCTGTCGCCAg                                                                                                   <  2:227605/251‑1 (MQ=255)
                                                                                                                                                                      cTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgc                                                                             >  1:312735/1‑251 (MQ=255)
                                                                                                                                                                        aCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGCCCGACCGCGTgc                                                                             >  1:18773/1‑249 (MQ=255)
                                                                                                                                                                               ttttGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTATTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTgg                                                                    >  1:211014/1‑251 (MQ=255)
                                                                                                                                                                                  tGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                                                                 >  1:183307/1‑251 (MQ=255)
                                                                                                                                                                                   gAATGTACGCGTCAGGTAGGTCGGTCCGCGGAT‑CGG‑TTCTTCAGCAGTAGTGGCG‑ACTTTCTCCTGTTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATTTCGGTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                                                                 <  2:68111/250‑1 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              >  2:435968/1‑251 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTAGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              >  1:457329/1‑251 (MQ=255)
                                                                                                                                                                                     aTTTACGCGGCAGGTAGGGCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGGGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                              <  2:298901/251‑1 (MQ=255)
                                                                                                                                                                                        tACGCGGCAGGTAGGCCTGGCCGAGGAT‑CGG‑TTATTCTTCCGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTCCGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGAACGTTACGGTTCATGTCGTTGGCTGTCGCCAGCGCATCAAGCCCGACCGGGTGCAGAACCTGGTGCCCCgtt                                                           >  2:485626/1‑249 (MQ=255)
                                                                                                                                                                                                aggtaggtCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                                                   >  1:125505/1‑251 (MQ=255)
                                                                                                                                                                                                aggtaggtCTGGCAGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                                                   >  2:228189/1‑251 (MQ=255)
                                                                                                                                                                                                  gtaggtCTGGCCGAGGAT‑CGG‑TTCTTCAGCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGTGCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCt                                                 <  1:46158/251‑1 (MQ=255)
                                                                                                                                                                                                    aggtCTGGCCGAAGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                                                <  2:118994/250‑1 (MQ=255)
                                                                                                                                                                                                        cTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                                            >  1:114680/1‑250 (MQ=255)
                                                                                                                                                                                                                     cGG‑TTCTTCATCAGCAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg                                >  2:252590/1‑250 (MQ=255)
                                                                                                                                                                                                                             tcatcaGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                         <  1:460780/250‑1 (MQ=255)
                                                                                                                                                                                                                              catcaGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaaa                        <  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)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          
AGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGC‑GTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAT‑CGG‑TTCTTCATCAGTAGTGGCG‑ACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGA  >  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: ‑