Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG34.9% 121.4 / 20.2 64V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (1/40);  new base G (22/0);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 4.36e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGG  >  NC_000913/941464‑941962
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
ccGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                                                                                                                             <  2:12449/251‑1 (MQ=255)
    taagCTGGGCGGCACCATGACCCGCTCCTGGACGCCGGGAAATACCCTGGACGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACACCTCCGGGCAAATTGCGGGAGGTTTGAACTATACCTCCCGCGGCCGGGCACATGGCAACAGACCGTCGGGAGTCGATAACAGTAACCTGGGAGCGCTGTTTGGTTATCACCCTGGTGCAAAGCGAATGAGTGGTGGGGGGGTgggcg                                                                                                                                                                                                                                                         >  2:193411/4‑247 (MQ=255)
         tGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACtttatat                                                                                                                                                                                                                                                    >  1:284455/1‑246 (MQ=255)
                      gACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                                                                                                                        <  1:102748/250‑1 (MQ=255)
                       aCCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                                                                                                                                      <  1:437409/251‑1 (MQ=255)
                         ccGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                                                                                                                                    <  2:193650/251‑1 (MQ=255)
                                 ggCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTTCGGCGACTACTCCTCCGCGCAAATTGCGTAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGTCAACAGCCCTTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCTAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCt                                                                                                                                                                                                                            <  1:153822/251‑1 (MQ=255)
                                    cgccgGGAAATACCATGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACGCAGCCGCGCAAATGGCGGAAGGTTTGAAATATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGAGAGCGAAAAAAGTAAGCTGGTAGTGCTGGTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaag                                                                                                                                                                                                                         <  2:48393/251‑2 (MQ=255)
                                                    ggTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                                                                                                                         <  2:3009/251‑1 (MQ=255)
                                                      gcgcgcgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                                                                                                                       <  2:7224/250‑1 (MQ=255)
                                                        gcgcgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                                                                                                                                     <  2:385344/251‑1 (MQ=255)
                                                           cgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                                                                                                                                  <  1:353943/251‑1 (MQ=255)
                                                                gatgaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcgcg                                                                                                                                                                                               <  2:205191/249‑1 (MQ=255)
                                                                        gctgcGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                                                                                                                                       <  1:144641/249‑1 (MQ=255)
                                                                        gctgcGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                                                                                                                                       <  2:26378/249‑1 (MQ=255)
                                                                           gcggcgGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGGGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑Acc                                                                                                                                                                                 <  2:190207/251‑1 (MQ=255)
                                                                              gcggcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGGCAGAAATCCAATACCCGCATGATAATCATAGATCCGCGCCATACCGAC‑CCCggg                                                                                                                                                                              >  1:73074/1‑250 (MQ=255)
                                                                                 gcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGcc                                                                                                                                                                           <  1:450845/251‑1 (MQ=255)
                                                                                  gtATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATTCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCg                                                                                                                                                                          <  2:424054/250‑1 (MQ=255)
                                                                                                      actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                                                                                                                      <  1:239323/251‑1 (MQ=255)
                                                                                                         actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                                                                                                                                   <  1:138450/251‑1 (MQ=255)
                                                                                                            cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAGCCCTGGAGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCAGCGATCCGCGCTATACCGAC‑ACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTa                                                                                                                                                  <  2:214902/249‑1 (MQ=255)
                                                                                                            cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTTAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCCTCCTCGATCCGCGCTATACCGAC‑ACCCGGGCCGGGGGCGAAGATGAGGGGGTCCCTAtt                                                                                                                                                >  1:25955/1‑251 (MQ=255)
                                                                                                               ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                                                                                                                               <  2:127682/249‑1 (MQ=255)
                                                                                                               ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                                                                                                                             <  1:29725/251‑1 (MQ=255)
                                                                                                               ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCGAC‑ACCGGGGCCGGGCGCGAAGAAGAGGGGAGcccctttcgt                                                                                                                                             >  1:189210/1‑244 (MQ=255)
                                                                                                               ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATTTTTAACACGCCCGCCCGAAATCTAACGCCCGCATTATCCTCCTCGCTCCCCGCTATACCGAC‑ACCCGTGCCGGGGGCGGAGATgtgggggtcccttatccc                                                                                                                                             >  1:53985/1‑234 (MQ=255)
                                                                                                                cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTAAAAAGAACGCCAGAAATCTAAAGCCCGCATGGACATTATCGATCCCGGCTATACCGAA‑ACCGGTGCCGGGCGCGGAGAGGAGGGGGTCCCTAtttcgc                                                                                                                                            >  1:92066/1‑247 (MQ=255)
                                                                                                                                  tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCCGAAAACTAATGCCCCCATGATCATCCTCTGTCCCCCCCATACCGCC‑ACCGGGGCCGGGGGCGAAAATGAGTGGGGCCCTATTCGT‑CCGGGGACAGATGCCGcgc                                                                                                                         >  1:363897/1‑249 (MQ=255)
                                                                                                                                          ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAAc                                                                                                                <  1:259912/251‑1 (MQ=255)
                                                                                                                                            aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGACAGGAAACTAATGGCCGCCTGGTCATCATCGAACCGCGCTATAACGAC‑ACCGGGGCCGGGGGCGGAGATGAGGGGATCCCCATTCCT‑CCGGGGAAAGATGCCGCACTGG‑TTaaagg                                                                                                              >  1:178536/1‑248 (MQ=255)
                                                                                                                                             ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTTGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAAc                                                                                                                <  2:320060/248‑1 (MQ=255)
                                                                                                                                             ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTTGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATTTTTAACAAGCACGCCCGAAATCTAATTCCCGCATGATCATCATCGATCCCCGCCATACCCAC‑ACCCGGGCCGGGGGCGGAGATGAGTGGGTCCCTATTTGT‑CCCGGGACAGATGCCGCCCTGG‑GTaaaaa                                                                                                              >  1:320060/1‑247 (MQ=255)
                                                                                                                                                  ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                                                                                                                                   <  1:26693/210‑1 (MQ=255)
                                                                                                                                                  ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGc                                                                                                        <  1:375828/251‑1 (MQ=255)
                                                                                                                                                            cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATg                                                                                              <  1:279417/251‑1 (MQ=255)
                                                                                                                                                            cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTTTCTTGAACAGGCCCGCCCGAAAACTTATACCCCCCAGGTCCTCATCCAGCCGCGCTCTACCCAC‑CCCCGTGCCGGGCGCGGAGATGAGTTGATCCCTCTTTGT‑CCGGGAACCGATGCCGCCCCGG‑GTAAACGTCTTGCGTACGTCATg                                                                                              >  2:162206/1‑251 (MQ=255)
                                                                                                                                                                 ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGCGACCTACTATTTTTAACAGGCACGGCCGAAATCGAATACCCGCAAGATCATCATCTATTCGCGCCTTATCGACCACCG‑CGCCCGGCGCCAAGATGAGGGGATCCCCATTCGTACCGGGAA‑AGCTACCGCTCTGG‑GTAACGGTCTGGGGTAAGTCATGagta                                                                                          >  2:91018/1‑247 (MQ=255)
                                                                                                                                                                  gCAACAGCCCGGCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACt                                                                                        <  2:186916/251‑1 (MQ=255)
                                                                                                                                                                   cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCGAC‑ACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTTGT‑CCGGGGAAAGATGCCGCCCTGG‑TTAAAGGTCTGGCGGACGTCATGAACACTg                                                                                       >  1:157534/1‑251 (MQ=255)
                                                                                                                                                                      cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATATTGTACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATACATTCGCGCCATACCGAA‑CCCGGCGCCGGGCGCGGAGATGAGGGGATCCCCCTTCGT‑CCGGGGACAACTGCCGCACTGG‑TTAACGGTCTGGCGGACGCCATGTTCAcctaa                                                                                     >  2:137950/1‑246 (MQ=255)
                                                                                                                                                                          ccGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCAC‑ACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCTATTCGT‑CCGGGGACAGATGCCGCCCTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAAccc                                                                                >  1:82016/1‑250 (MQ=255)
                                                                                                                                                                              cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAGCCtggt                                                                             <  1:45388/250‑1 (MQ=255)
                                                                                                                                                                               ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCCCGCCAGAAATCTAATGACCCCATGATCATCATCGATCCGCGCCATACCGAC‑AACGGTGCCGGGGGCGAAGATGAGTGGATCCCTATTTTT‑CCGGGGACAGATGCCGCCCTGG‑TTAACGGGCTGGCGGACGTCCTTATTACTGAAAAACtggtgg                                                                           >  1:242044/1‑251 (MQ=255)
                                                                                                                                                                                     cGAAAACAGTAAGCTGGTAGTGCTGTTTGGCAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCATCCTTTCCTTAACAGGGACGGCCGAAATCTTATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGGCGGGCGGGAAGATGGGTGGGTCCCTATTTGT‑CCGGGGAAAGATGCCGCCCTGG‑TTAACCGGCTGGCGTACGTCATGATTACTGAAAACCTGGTGGATCAgg                                                                     >  1:368785/1‑251 (MQ=255)
                                                                                                                                                                                                 gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATa                                                         <  1:251290/251‑1 (MQ=255)
                                                                                                                                                                                                      tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAACAAGCACGCCAGAAATATAATGCCCGCATGATCATCCTCGCTCCGCGCTATACCGAC‑ACCCGTGCCGGGCGCGGGGATGTGGGGGGCCCTATTAGT‑CCGGCAACAGAAGCCCCACTGGGTTAA‑GGGCTGGCGTACGGCATGATCACTGAAAAACTTGTGGAGCAGGGATTCCTCTATAAata                                                     >  2:391370/1‑250 (MQ=255)
                                                                                                                                                                                                       aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGaaaa                                                                                   <  1:5070/219‑1 (MQ=255)
                                                                                                                                                                                                          gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCCCATGATAAACATCGGTCCGCGCTATACCGAC‑ACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCTATTCGT‑TCGGGGACAGATGCCGCCCTGG‑TTTACGGTCTGGCGTACGGCCTGATCACTGAAAAACTGGGGGGTCAGGCATTTCTCGATAAATAtttcg                                                >  2:446417/1‑248 (MQ=255)
                                                                                                                                                                                                          gCTGTTTGGTAATAACCCTGGCAAACCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATGGATCCGCGCTATCCCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGC‑CCGGGAACAGATGCCGCAGTGG‑TTAACGGTCTGGCGTACGTCATTATCGCTAAAAACGGGGTGGACCAGGCATTCCTCGTTAAAGAttccg                                                <  2:197233/251‑4 (MQ=255)
                                                                                                                                                                                                             gTTTGGTAATATCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTg                                             <  1:294340/251‑1 (MQ=255)
                                                                                                                                                                                                                tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTTATATCTTTAAAAGGGACCCCAGAAATCTAATACCCGCCTGGTCATTATTGATACTCGCCATAACGAC‑ACCGGGGGCGGGGGCGGAGATGAGTGGATCCCTCTTTGT‑CCGGGGAAAGAATCCGCCCTGG‑TTAAAGGTCTGGCGGACGTCATGATCACTGGAAACCAGGTGGATCAGGGATTCCTCGATAAATATTGCGTTGGCt                                          >  1:386141/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCt                                          <  1:193411/250‑1 (MQ=255)
                                                                                                                                                                                                                 ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCGTACTATCTTTTAAATGCACGCCAGAAAACTAATGCCCCCATGAGAATCATTGAAACGCGCTATATCGAC‑ACCCGGGTCGGGCGCGAAGGTGTGTGGGTCCCTATTCTT‑CCCGGAACAGAGGCCCCACCGG‑TTTACGGGCTGGCGGACGTCATGATAACTGAAAACCCGGTGGGTCAGGCATTCCCCGCTAAATAtttaggtggctc                                         >  2:449329/1‑241 (MQ=255)
                                                                                                                                                                                                                   taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAGGCACGCCCCAAAACTAAAACCCGCATGATTATCATCGATCCGCGCTATACCCAC‑ACCCGGGGCGGGCGCGAAGATGAGTGGATACCCATTCGT‑CCCGGGACAGAAGCCGCACTGG‑GTAAAGGGCTGGCGGACGTCCAGAGAACTGAAAAACTGGTGGATCAGGGATTTCTTGATAAAATTTTCGTTGGCTACg                                       >  2:452691/1‑251 (MQ=255)
                                                                                                                                                                                                                    aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACg                                       <  1:227142/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGaa                                <  1:162206/251‑1 (MQ=255)
                                                                                                                                                                                                                           cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTCTCTTTCACAAGCACGGCAGAAAACCAATGCCCGCATGATCAATATCGGTACCAGGTATACCGAA‑ACCGGGGCCGGGGGCGAAGAAGAGGGGATCCCCATTCTG‑CCGGGAACAGACGCCGCCCTGG‑TTAACGGTCTGGCGTACGGCACGATCACTGTAAACCTAGGGGAGCAGGCATTCCTTCAGAAATATTGGGTTGTCTACGATTAGaa                                >  1:260828/1‑250 (MQ=255)
                                                                                                                                                                                                                              gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAGGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAAcc                            <  2:445901/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAACAGGCACGCCCAAAATCTAATGACCGCATTATCATCATCGATCCGCGCTATACCGAC‑ACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGGACAGATGCCGCACTGG‑TTAAAGGTCTGGCGTACGTCCTGATCACTGAAAACCTGGGGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTg                         >  1:328984/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAGCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgc                        <  2:242044/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTg               <  2:25955/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATg   <  1:446417/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATgg  <  1:137950/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
CCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGAC‑ACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGT‑CCGGGAACAGATGCCGCACTGG‑TTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGG  >  NC_000913/941464‑941962

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