Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7270GT25.4% 149.1 / 10.1 67E257* (GAA→TAA) dmsAdimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base G (6/44);  new base T (17/0);  total (23/44)
Fisher's exact test for biased strand distribution p-value = 2.99e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTT  >  NC_000913/941486‑941976
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            
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)
                                                     gcggcgGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGGGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                                                                                                                                               <  2:190207/251‑1 (MQ=255)
                                                        gcggcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGGCAGAAATCCAATACCCGCATGATAATCATAGATCCGCGCCATACCGACCCCggg                                                                                                                                                                                            >  1:73074/1‑250 (MQ=255)
                                                           gcTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                                                                                                                                                         <  1:450845/251‑1 (MQ=255)
                                                            gtATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATTCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                                                                                                                                        <  2:424054/250‑1 (MQ=255)
                                                                                actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                                                                                                                                                    <  1:239323/251‑1 (MQ=255)
                                                                                   actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                                                                                                                                                 <  1:138450/251‑1 (MQ=255)
                                                                                      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAGCCCTGGAGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCAGCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTa                                                                                                                                                                <  2:214902/249‑1 (MQ=255)
                                                                                      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTTAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCCTCCTCGATCCGCGCTATACCGACACCCGGGCCGGGGGCGAAGATGAGGGGGTCCCTAtt                                                                                                                                                              >  1:25955/1‑251 (MQ=255)
                                                                                         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                                                                                                                                             <  2:127682/249‑1 (MQ=255)
                                                                                         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                                                                                                                                           <  1:29725/251‑1 (MQ=255)
                                                                                         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCGACACCGGGGCCGGGCGCGAAGAAGAGGGGAGcccctttcgt                                                                                                                                                           >  1:189210/1‑244 (MQ=255)
                                                                                         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATTTTTAACACGCCCGCCCGAAATCTAACGCCCGCATTATCCTCCTCGCTCCCCGCTATACCGACACCCGTGCCGGGGGCGGAGATgtgggggtcccttatccc                                                                                                                                                           >  1:53985/1‑234 (MQ=255)
                                                                                          cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTAAAAAGAACGCCAGAAATCTAAAGCCCGCATGGACATTATCGATCCCGGCTATACCGAAACCGGTGCCGGGCGCGGAGAGGAGGGGGTCCCTAtttcgc                                                                                                                                                          >  1:92066/1‑247 (MQ=255)
                                                                                                            tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCCGAAAACTAATGCCCCCATGATCATCCTCTGTCCCCCCCATACCGCCACCGGGGCCGGGGGCGAAAATGAGTGGGGCCCTATTCGTCCGGGGACAGATGCCGcgc                                                                                                                                       >  1:363897/1‑249 (MQ=255)
                                                                                                                    ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc                                                                                                                              <  1:259912/251‑1 (MQ=255)
                                                                                                                      aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGACAGGAAACTAATGGCCGCCTGGTCATCATCGAACCGCGCTATAACGACACCGGGGCCGGGGGCGGAGATGAGGGGATCCCCATTCCTCCGGGGAAAGATGCCGCACTGGTTaaagg                                                                                                                            >  1:178536/1‑248 (MQ=255)
                                                                                                                       ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTTGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc                                                                                                                              <  2:320060/248‑1 (MQ=255)
                                                                                                                       ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTTGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATTTTTAACAAGCACGCCCGAAATCTAATTCCCGCATGATCATCATCGATCCCCGCCATACCCACACCCGGGCCGGGGGCGGAGATGAGTGGGTCCCTATTTGTCCCGGGACAGATGCCGCCCTGGGTaaaaa                                                                                                                            >  1:320060/1‑247 (MQ=255)
                                                                                                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                                                                                                                                                 <  1:26693/210‑1 (MQ=255)
                                                                                                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc                                                                                                                      <  1:375828/251‑1 (MQ=255)
                                                                                                                                      cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATg                                                                                                            <  1:279417/251‑1 (MQ=255)
                                                                                                                                      cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTTTCTTGAACAGGCCCGCCCGAAAACTTATACCCCCCAGGTCCTCATCCAGCCGCGCTCTACCCACCCCCGTGCCGGGCGCGGAGATGAGTTGATCCCTCTTTGTCCGGGAACCGATGCCGCCCCGGGTAAACGTCTTGCGTACGTCATg                                                                                                            >  2:162206/1‑251 (MQ=255)
                                                                                                                                           ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGCGACCTACTATTTTTAACAGGCACGGCCGAAATCGAATACCCGCAAGATCATCATCTATTCGCGCCTTATCGACCACCGCGCCCGGCGCCAAGATGAGGGGATCCCCATTCGTACCGGGAAAGCTACCGCTCTGGGTAACGGTCTGGGGTAAGTCATGagta                                                                                                        >  2:91018/1‑247 (MQ=255)
                                                                                                                                            gCAACAGCCCGGCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACt                                                                                                      <  2:186916/251‑1 (MQ=255)
                                                                                                                                             cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTTGTCCGGGGAAAGATGCCGCCCTGGTTAAAGGTCTGGCGGACGTCATGAACACTg                                                                                                     >  1:157534/1‑251 (MQ=255)
                                                                                                                                                cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATATTGTACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATACATTCGCGCCATACCGAACCCGGCGCCGGGCGCGGAGATGAGGGGATCCCCCTTCGTCCGGGGACAACTGCCGCACTGGTTAACGGTCTGGCGGACGCCATGTTCAcctaa                                                                                                   >  2:137950/1‑246 (MQ=255)
                                                                                                                                                    ccGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCTATTCGTCCGGGGACAGATGCCGCCCTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAAccc                                                                                              >  1:82016/1‑250 (MQ=255)
                                                                                                                                                        cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAGCCtggt                                                                                           <  1:45388/250‑1 (MQ=255)
                                                                                                                                                         ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCCCGCCAGAAATCTAATGACCCCATGATCATCATCGATCCGCGCCATACCGACAACGGTGCCGGGGGCGAAGATGAGTGGATCCCTATTTTTCCGGGGACAGATGCCGCCCTGGTTAACGGGCTGGCGGACGTCCTTATTACTGAAAAACtggtgg                                                                                         >  1:242044/1‑251 (MQ=255)
                                                                                                                                                               cGAAAACAGTAAGCTGGTAGTGCTGTTTGGCAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCATCCTTTCCTTAACAGGGACGGCCGAAATCTTATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGGCGGGCGGGAAGATGGGTGGGTCCCTATTTGTCCGGGGAAAGATGCCGCCCTGGTTAACCGGCTGGCGTACGTCATGATTACTGAAAACCTGGTGGATCAgg                                                                                   >  1:368785/1‑251 (MQ=255)
                                                                                                                                                                           gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATa                                                                       <  1:251290/251‑1 (MQ=255)
                                                                                                                                                                                tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAACAAGCACGCCAGAAATATAATGCCCGCATGATCATCCTCGCTCCGCGCTATACCGACACCCGTGCCGGGCGCGGGGATGTGGGGGGCCCTATTAGTCCGGCAACAGAAGCCCCACTGGGTTAAGGGCTGGCGTACGGCATGATCACTGAAAAACTTGTGGAGCAGGGATTCCTCTATAAata                                                                   >  2:391370/1‑250 (MQ=255)
                                                                                                                                                                                 aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGaaaa                                                                                                 <  1:5070/219‑1 (MQ=255)
                                                                                                                                                                                    gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCCCATGATAAACATCGGTCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCTATTCGTTCGGGGACAGATGCCGCCCTGGTTTACGGTCTGGCGTACGGCCTGATCACTGAAAAACTGGGGGGTCAGGCATTTCTCGATAAATAtttcg                                                              >  2:446417/1‑248 (MQ=255)
                                                                                                                                                                                    gCTGTTTGGTAATAACCCTGGCAAACCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATGGATCCGCGCTATCCCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGCCCGGGAACAGATGCCGCAGTGGTTAACGGTCTGGCGTACGTCATTATCGCTAAAAACGGGGTGGACCAGGCATTCCTCGTTAAAGAttccg                                                              <  2:197233/251‑4 (MQ=255)
                                                                                                                                                                                       gTTTGGTAATATCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTg                                                           <  1:294340/251‑1 (MQ=255)
                                                                                                                                                                                          tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGGCTTTATATCTTTAAAAGGGACCCCAGAAATCTAATACCCGCCTGGTCATTATTGATACTCGCCATAACGACACCGGGGGCGGGGGCGGAGATGAGTGGATCCCTCTTTGTCCGGGGAAAGAATCCGCCCTGGTTAAAGGTCTGGCGGACGTCATGATCACTGGAAACCAGGTGGATCAGGGATTCCTCGATAAATATTGCGTTGGCt                                                        >  1:386141/1‑251 (MQ=255)
                                                                                                                                                                                           ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCt                                                        <  1:193411/250‑1 (MQ=255)
                                                                                                                                                                                           ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTCGTACTATCTTTTAAATGCACGCCAGAAAACTAATGCCCCCATGAGAATCATTGAAACGCGCTATATCGACACCCGGGTCGGGCGCGAAGGTGTGTGGGTCCCTATTCTTCCCGGAACAGAGGCCCCACCGGTTTACGGGCTGGCGGACGTCATGATAACTGAAAACCCGGTGGGTCAGGCATTCCCCGCTAAATAtttaggtggctc                                                       >  2:449329/1‑241 (MQ=255)
                                                                                                                                                                                             taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAAAAGGCACGCCCCAAAACTAAAACCCGCATGATTATCATCGATCCGCGCTATACCCACACCCGGGGCGGGCGCGAAGATGAGTGGATACCCATTCGTCCCGGGACAGAAGCCGCACTGGGTAAAGGGCTGGCGGACGTCCAGAGAACTGAAAAACTGGTGGATCAGGGATTTCTTGATAAAATTTTCGTTGGCTACg                                                     >  2:452691/1‑251 (MQ=255)
                                                                                                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACg                                                     <  1:227142/250‑1 (MQ=255)
                                                                                                                                                                                                    ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGaa                                              <  1:162206/251‑1 (MQ=255)
                                                                                                                                                                                                     cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTCTCTTTCACAAGCACGGCAGAAAACCAATGCCCGCATGATCAATATCGGTACCAGGTATACCGAAACCGGGGCCGGGGGCGAAGAAGAGGGGATCCCCATTCTGCCGGGAACAGACGCCGCCCTGGTTAACGGTCTGGCGTACGGCACGATCACTGTAAACCTAGGGGAGCAGGCATTCCTTCAGAAATATTGGGTTGTCTACGATTAGaa                                              >  1:260828/1‑250 (MQ=255)
                                                                                                                                                                                                        gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAGGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAAcc                                          <  2:445901/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAACAGGCACGCCCAAAATCTAATGACCGCATTATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGCCGCACTGGTTAAAGGTCTGGCGTACGTCCTGATCACTGAAAACCTGGGGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTg                                       >  1:328984/1‑251 (MQ=255)
                                                                                                                                                                                                            aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAGCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgc                                      <  2:242044/251‑1 (MQ=255)
                                                                                                                                                                                                                     gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTg                             <  2:25955/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATg                 <  1:446417/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATgg                <  1:137950/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTGGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCa             <  2:284455/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATaa       <  2:189210/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCGCTAGAAAGCt   <  1:296124/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAg     <  1:47605/248‑1 (MQ=255)
                                                                                                                                                                                                                                                tGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCtt  >  2:177295/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCt   <  2:344452/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            
GACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGTACGTCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAATGGCCACTATAAAGCTT  >  NC_000913/941486‑941976

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