Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6050TC35.2% 101.5 / 14.9 55G187G (GGA→GGGycfTinner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base T (32/3);  new base C (0/19);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 8.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTC  >  NC_000913/1174356‑1174842
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
cagcagGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATACCACTGTTCCTTc                                                                                                                                                                                                                                                   >  2:8994/1‑251 (MQ=255)
 agcagGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAACTTGCGGATCACACTGTTCATTcc                                                                                                                                                                                                                                                  >  2:713025/1‑251 (MQ=255)
   cagGCGGGTGGAGCGCATACCGAATCGTTGCTCGTATGGATAAAACAGCTTCAGGATCACCCCAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGAGCACACGGTTCATTcccc                                                                                                                                                                                                                                                <  2:122959/251‑1 (MQ=255)
      gCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCATATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATGCCccacc                                                                                                                                                                                                                                             <  2:498929/251‑1 (MQ=255)
      gCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAAGCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                                                                                                                             <  2:41527/251‑1 (MQ=255)
                                      aTTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                                                                                                                          >  2:256579/1‑251 (MQ=255)
                                                  aagccgtgtgctctcccaatctccccAGCGACAACAGCAGCGAGATAGAAAACAACCGGCCCCCGCCCGCCAGAAGGGTAAGCAAAAAACCCATCAGCAGAGGGCGGCGCAACGGCCCCTCTTTAACACGGGTCATTATTGTTGCGCAAAACCATGCGCCAGGGCTTAAAAAAGGCAAATTGCGGATCACACTGTTATTCCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                              <  1:69769/227‑1 (MQ=255)
                                                    ttCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACCACCAGGCACCGACCGCCAGAACGGTAAGCAAAGACCCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                                                                                                             >  1:82291/1‑250 (MQ=255)
                                                            caccacAATCGCTACCAGCGACAACAGCAGCGCGATAGTACACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGACGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCCAGGCTGTAATAAGGCACATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                                                                                                                     >  2:455047/1‑250 (MQ=255)
                                                             accacAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                                                                                                                     >  1:340916/1‑249 (MQ=255)
                                                               cacaATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTg                                                                                                                                                                                 >  1:30076/1‑251 (MQ=255)
                                                                acaATCGATACCAGCGACAACCGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCCGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAgggg                                                                                                                                                                                >  2:366938/1‑248 (MQ=255)
                                                                 caATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                                                                                               >  1:16123/1‑251 (MQ=255)
                                                                  aaTCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                                                                                               >  2:550735/1‑250 (MQ=255)
                                                                     cGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGTGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcg                                                                                                                                                                           >  1:378003/1‑251 (MQ=255)
                                                                        tACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                              >  2:69769/1‑229 (MQ=255)
                                                                           cAGCAACAACAGCAGGGAAAAAGAAAACACCCAGGCCCCGCCCCCCAGAACGGAAAACAAAGAAGCCAACAGCAGAGGGGGGCCCAACGGCCCCTCTTAAACAAGGGTCATTATTGTGGCGCAAACCCACGCGCAAAGGTTGAAATAAGGAAAATGGCGGACCCCACTTTTCTTCCCCCCCCACGGGGGGGAAAAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                                                                                                                     <  2:560973/251‑1 (MQ=255)
                                                                              cGACAACAGCAGCGAGATCGTAACCAACCAGGCACCGACCGCCAGAACGGCAAGCAAAGAAGCCATCAGCAGCTGGCGGCGCAACGGCACCTCTTTAACACAGGTCCTTATTGTTGCGCCACACCATGCGCCAAGGCTGTCATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACACAGCAGTACAAACAAGGCGAATAGGGGCAGCGCCAGACGGc                                                                                                                                                                  >  2:127723/1‑251 (MQ=255)
                                                                                  aacaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                                                                                                               >  1:80879/1‑250 (MQ=255)
                                                                                         gcGAGATAGTAAACAACCAGGCACCGACCGCCAGACCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCCACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCACGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAACCAAGGCGAATAGTGGCAGGGCCCGACGGCTAAAAAtttt                                                                                                                                                        >  2:304030/1‑250 (MQ=255)
                                                                                           gagaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                                                                                                                                      >  1:548001/1‑250 (MQ=255)
                                                                                               aaGTAAACAACCAGGCACCCCCCGCCAGAACGGTAAGAAAAGAGCCCATCAGCAGATGGCGGCGCACCGCCACCTCTTAAACACAGGTCATTATTGTTGCGCCAAACCAGGCGCCAGGGCTGAAAAAAGGCAAATGGCGGACCACACTGTTCATCCCCCCCCGCGGCGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                                                                                                                 <  2:80879/250‑1 (MQ=255)
                                                                                                    aacaacCAGGCACCGCCCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCaaga                                                                                                                                             >  2:112670/1‑247 (MQ=255)
                                                                                                                     cgccgaaaCGGTAAGCAAAAAGCCCAGCAGCAGAGGGCGGCCCAACGGCTCCTCTTTAACACTGGTTATTTTTGTTGCGCCAAACCAGGCGCCAAGGCTTAAATAAGGAAAATTGCAGAACACCCTGTTCCTCCCCCGCCAGGCGGTGAGAAGAAAATCAACCGCCACTCTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCATACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                                                                                                           <  1:69703/245‑1 (MQ=255)
                                                                                                                       ccAGAAGGGAAAGAAAAGAAACCAACAGCAGGTGGGGGCGCAACGGCACCTCTTAAACACAGGTAATTTTTGTGGCCCAAAACCATGCGCCAGGGCTGTAAAAAGAAAAATTGCGGATCACCCTGTTCATCCCCCCCCACGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                                         <  1:8994/251‑1 (MQ=255)
                                                                                                                          aaaagggaaaaaaaaaaagcaagcagcagATGGGGGCCAAACGCCCCCTTTTAAACACGGTTCATTATTGTGGCCCAAAACCAGGCGCCAAGGTTGTAAAAAGGAAAATGCGGGATCACACTGTTCTGCCCCCCCCACGGCTGGGAAACGAATTAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                                                                                                                       <  1:612799/228‑1 (MQ=255)
                                                                                                                                aaaGCAAAAAAGCCACCGGCAGGTGGCGGCGCACCGGCACCTCTTTAAAACAGGTCATTATTGTGGCCCCAAACCATGCGCCAGGGCGTTAAAAAGGCAAATGGCGAATCACTCTGTTAATCCCCCCCCCCGGCGTGGAAACGAAATTTACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTGACATACCACGAAATAGACAATTAACGCATACAGATACcac                                                                                                                <  1:455047/250‑1 (MQ=255)
                                                                                                                                      aagaagCCAGAAGAAGAGGGGGGCGCAACGGCACCTCTATAACACAGATAATTTTTGTTGCCCCAACCCGTGCGCAAAGGCTGAAAAAAGACAAATGGCGGATCACCTTTTTCATCCCCCCCCAGGGCGGGGGAACAAAATTATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg                                                                                                          <  2:251056/251‑1 (MQ=255)
                                                                                                                                                  agcagATGGGGGAGGAGCGGCACCTCTTAAACACAGGTAATTTTGGTTCCGCAAAACCAGGCGCCAAGGCTGTAAAAAGGCAAATGGCGGATCACACTGTTCATCCCCCACCACGGGGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                                                                                               <  2:548001/251‑1 (MQ=255)
                                                                                                                                                  acgagATGGGGGCGCAACGGCCCCTTTTAAACACAGTTATTTATTGTTGCGCAAAACCATGCGCCAAGGCTGAAAAAAGGAAAATTGCGGATCCCACTTTTCATCCCCCCCCACGGCGTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACGTACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                                                                                               <  1:698648/248‑1 (MQ=255)
                                                                                                                                                      gagggcggcGCACCGGCCCCTCTTAAACAAAGGTCATTATTGTTGCGCAAACCCATGCCCAAAGGCTGTAAAAAGGAAAATGGCGGATCCCCCTTTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc                                                                                          <  1:256579/248‑1 (MQ=255)
                                                                                                                                                         ggcggcGCACCGCCCCCTCTTTAACACGGTTCATTATTGTGGCGCAAACCCATGCGCCAAGGCTGTAAAAGGGCAAATGGCGGATCACACTGTTCATCCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatcatc                                                                                       <  1:184559/251‑1 (MQ=255)
                                                                                                                                                                 aaCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGc                                                                                 >  1:484256/1‑249 (MQ=255)
                                                                                                                                                                              tAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTACTAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAACATCCCGTGCCGGAACTCACCgggg                                                                 >  1:272980/1‑249 (MQ=255)
                                                                                                                                                                              tAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGATGTCATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCCGCGCCAGACGGCTAAAAATTTTACATACCACGACATAGCCAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGAGATCATCCCGTGCAGGAACTCACCGGTg                                                                 >  1:491638/1‑251 (MQ=255)
                                                                                                                                                                                acacAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGACCGAACTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTggg                                                               >  1:339201/1‑250 (MQ=255)
                                                                                                                                                                                acacAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATcccctgacggaacccacccggggg                                                                >  2:93041/1‑230 (MQ=255)
                                                                                                                                                                                 cacaGTTCATTATTGTTGCGCAAACCCAGGCCCCAAGGCTGAAAAAGGGCAAATTGCGGATCACACTTTTCTTCCCCCCCCGCGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTc                                                                        <  1:147707/241‑1 (MQ=255)
                                                                                                                                                                                 cacaGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGGTCATTCCGTGGAGGAACTCaaccgtggg                                                               >  1:403069/1‑242 (MQ=255)
                                                                                                                                                                                 cacaGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGATCATCCCGTGCAAGAACTc                                                                        >  2:147707/1‑241 (MQ=255)
                                                                                                                                                                                   caGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCCGGAACTCACCggggggatt                                                            >  2:629679/1‑244 (MQ=255)
                                                                                                                                                                                           tattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGACTAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACAACAAGCAGGCGCTGGCGGAGAACATCCCGGGCAGGAACCTCCCGGGGGAACCGGCATa                                                      >  1:517866/1‑249 (MQ=255)
                                                                                                                                                                                           tattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGGGCTGGCGGGGCATCACCCGGGCAGGAAcacccccggggaaactgcataa                                                      >  2:732270/1‑229 (MQ=255)
                                                                                                                                                                                               gttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCTGCAGGa                                                                             >  1:475332/1‑223 (MQ=255)
                                                                                                                                                                                                ttgCGCCAAACCATGCGCCAAGGGTGTAAAAAGGCAAATTCCGGACCACACTGTTCATCCCCCACCACGGCGGGGAAACAAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGc                                                <  1:304030/250‑1 (MQ=255)
                                                                                                                                                                                                      aaaaCAATGCCCCAAGGCTGAAAAAAGGAAAATTGCGGAGCCCTTTGTCCTTCCCCCCCCCCGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAATCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAg                                         <  1:484581/250‑1 (MQ=255)
                                                                                                                                                                                                        aacctgacccaaggggtgaaaaaaGGCAAATTGCGGACCCCTCTGTTCTTCCCCCCCCAGGGGGTGAGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAGTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAGCATCCCGTGCAGGAACTCGCCGGTGGAATCGGCATTAGCGGCATTGGAGGc                                       <  2:82291/229‑1 (MQ=255)
                                                                                                                                                                                                           cATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCAATCCGGGCAGGAAACCACCGGGGGAATCGGCATAAGCGGCATTTGAGGCAtt                                    >  1:697572/1‑251 (MQ=255)
                                                                                                                                                                                                                         gaaaaaaGGAAAATTGCGGATCACCCTTTTCATCCCCCCCCACGGCGGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATcgcg                        <  1:629679/244‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaaaGGAAAATTGGGGATCCCCCTGTTCTTCCCCCCCCGCGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAAGCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTc                     <  1:233150/247‑1 (MQ=255)
                                                                                                                                                                                                                              aaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCATTGGAGGCCTTGCTTAAATCGCGCTCAggg                 >  1:100121/1‑250 (MQ=255)
                                                                                                                                                                                                                                 gaaaaTTGCGGATCCCACTGTTCATCCCCCCCCACGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgcc              <  1:550735/249‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTTGCGGATCACCCTGTACATCCCCCCCCCCGGGGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAGCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTgccag            <  1:112670/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGGAACGGCATAAGCGGCATTGGAGGCATTGCTTAAAACGAGCCCAGGCgccaccccccggttt  >  1:55752/1‑240 (MQ=255)
                                                                                                                                                                                                                                                acacTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTgggtc  >  1:237142/1‑246 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
CAGCAGGCTGGTGGAGCGCATACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTC  >  NC_000913/1174356‑1174842

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