Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,6260GA23.6% 114.7 / 12.6 55F180F (TTC→TTTycfTinner membrane protein YcfT
Reads supporting (aligned to +/- strand):  ref base G (35/7);  new base A (0/13);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 5.34e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCC  >  NC_000913/1174376‑1174876
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           
tACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                                                                                                                                             >  2:15755/1‑251 (MQ=255)
  ccGAATCGTTGCTCGTATTTTTAAAAGAGCTTTATGAGCACCACAACCGATACCAGCGACAACAACAGCGAGATATTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTGAACACAGGTCATTATGGTTGCGCCATACCAGGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCGCGGCGTGGTAACGaa                                                                                                                                                                                                                                                           <  2:424898/251‑1 (MQ=255)
    gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTNTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                                                                                                                                                         >  2:208425/1‑251 (MQ=255)
    gAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAt                                                                                                                                                                                                                                                         >  1:150765/1‑251 (MQ=255)
         gTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAAcc                                                                                                                                                                                                                                                    >  1:162993/1‑251 (MQ=255)
           tGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCTTGGGAACGAAATTAACCGc                                                                                                                                                                                                                                                  >  1:22017/1‑251 (MQ=255)
                  aTTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTc                                                                                                                                                                                                                                           >  2:54659/1‑251 (MQ=255)
                               cTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaac                                                                                                                                                                                                                               >  1:114051/1‑250 (MQ=255)
                                   aTGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGc                                                                                                                                                                                                                          >  1:51395/1‑251 (MQ=255)
                                     gATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCTAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                                                                                                                                                         >  1:155190/1‑250 (MQ=255)
                                              aaTCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGc                                                                                                                                                                                                                >  1:431092/1‑250 (MQ=255)
                                                          cgacaccggaggcGAGATAGAAAACAACAAGGCACCGACCGCAAGAACGGAAAGCAAAGAAGCCATCAGCAGAGGGCGGCGCAACGGCCCCTCTTTAACACAGGTCATTATTGTTGCGCAAACCCAGGCGCCAAGGCTGTAATAAGGCAAATGGCGGATCACACTGTTCATCCCCCCCCCCGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                                                                                                                                                   <  1:29440/240‑1 (MQ=255)
                                                          cGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACNAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                                                                                                                                                   >  2:212780/1‑251 (MQ=255)
                                                           gACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGCCGGCt                                                                                                                                                                                                  >  2:361185/1‑251 (MQ=255)
                                                           cacaacaGCAGCGAGATAGAAAACACCCAGGCACCGCCCGCAAAAACGGTAAGCAAAGAAGCCACCAGCAGATGGCGGCGCACCGCCACCTTTTAACCACAGGTCTTTATTGTTCCGCCAAACCATCCGCCAAGGCTGTAAAAAGGAAAATTGCGAATCACCCTGTTCATCCCCCCCCCGGGCGTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                                                                                                                                                  <  1:146640/250‑1 (MQ=255)
                                                               acaGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGACATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                                                                                                                                               >  2:305701/1‑250 (MQ=255)
                                                                  gcagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                                                                                                                                           >  1:332898/1‑251 (MQ=255)
                                                                   cggggagaaaaaaaacaccaagcccccgaccgCCAGAAGGGAAAGCAAAAAAGCCAACAGAAGAGGGCGGCACACCGCCACCTTTATACCACAGTTCATTTTTGTTGCGCAAAACCAGGCGCAAAGGCTGTAATAAGGCAAAGTGCGGACCACACTCTTCATCCCCCACCAGGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                                                                                                                                          <  2:244140/226‑1 (MQ=255)
                                                                   cagcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                                                                                                                                          >  2:378774/1‑251 (MQ=255)
                                                                    agcGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGCACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                                                                                                                                                         >  1:76283/1‑251 (MQ=255)
                                                                        agaTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACAt                                                                                                                                                                                     >  1:161293/1‑251 (MQ=255)
                                                                          aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGTAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAACTTTTACATAc                                                                                                                                                                                   >  1:465800/1‑251 (MQ=255)
                                                                          aTAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGAAACGGCACCTCTTTAACACCGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTAAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAc                                                                                                                                                                                   >  1:465820/1‑251 (MQ=255)
                                                                            agaaaaCACCCAGGCCCCGACCGCCAAAACGGTAAGCAAAGAAGCCACCAGCAGATGGCGGCGCAACGGCACCTCTTAACCACAGGTCTTTTTTGTTGCGCCAACCCATGCGCCAAGGCTGAAAAAAGGCAAATGGCGGACCCCACTTTTCTTCCCCCCCCCCGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                                                                                                                                                 <  1:378774/248‑1 (MQ=255)
                                                                                 acaacCAGCCCCCGCCCGCCAGAACGTTAAGCAAAAAAGCCATCAGCAGATGGCGGCGCACCGGCCCCTCTTTAACACAGGTTATTATTGTTGCGCAAAACCATGCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCAGGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaaa                                                                                                                                                                            <  1:208425/251‑1 (MQ=255)
                                                                                   aacCAGGCACCGACCGCCAGAACGGTAAGCAAAGCAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACgcaaa                                                                                                                                                                           >  2:360963/1‑246 (MQ=255)
                                                                                                   caaaaagggaaaacAAAGAACCCACCGGCAGGTGGCGGCGCACCGGCCCCTCTTAAACAAAGGTATTATTTGTTGCGCAAACCCATGCGCCAGGGCTGTAATAAGGCAAATGGCGGACCACCTTGTTCTTCCCCCACCACGGCGGGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                                                                          <  2:401636/238‑1 (MQ=255)
                                                                                                   caaGAACGGAAAGCAAAGAACCCATCAGCAGAGGGCGGCGCAACGGCCCCTCTTAAAAACAGGTCATTATTGTGGCGCCAAACCATGCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACACTGTTCACCCCCCCCCACGGGGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                                                                                                                                          <  2:161293/249‑1 (MQ=255)
                                                                                                               gCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCAACTCTTTAACACAGGAATTTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTAATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCCGACGGCTAAAAATTTTACATACCACGAAAAAGACAATTAACGCCTACATATACCACaac                                                                                                                                              >  2:464879/1‑250 (MQ=255)
                                                                                                                aaaagacgccaccagtagtgggcggcGAAACGGCACCTTTTAACAACAGGTTTTTTGTGTGGCGCAAAACCATGCCCCAGGGCTGTAAAAGGTCAATTTCCGGATCACTCTTTTCTTCCCCCCCCACGGCGGGGAAACAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGc                                                                                                                                             <  1:373255/231‑1 (MQ=255)
                                                                                                                     aagtcatcacatgtgggctgcccAACGCCACCTCTTAAAAACAGGTTATTATTGTTGCGCCAATCCCTGCGCCAAGTTTGAAAAAAGCAAATTTGCGATCCCCCCTGTCCATCCCCCCCCACGGCGTGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTATCGCATACAGATACCACAAGCTGGTg                                                                                                                                        <  2:325332/229‑1 (MQ=255)
                                                                                                                         cATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg                                                                                                                                    >  1:129647/1‑251 (MQ=255)
                                                                                                                             agcggATGGCGGCCCAACGCCCCCTTTTTAACCCAGGTAATTTTTGTGGCCCAAACCCATCCCCCAGGGCTGTAATAAGGCAATTGGCGGATCCCCCTGTTCATTCCCCCCCTCGGCGTGGGAAGAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt                                                                                                                                <  1:56608/247‑1 (MQ=255)
                                                                                                                                agagggcgcgacagcgtcccctttttAACACNGGACATTATGGTTGCGCCAAACCACGCGCAAAGGCTAAAAGAAACAAAGTGGCGGATCCCACTTTTCCTCCCCCACCACGGGGGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat                                                                                                                             <  2:204858/228‑1 (MQ=255)
                                                                                                                                    ggcggcGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAAACAGATACCACAAGCGGGTGCTGGCGGCGatcatc                                                                                                                         >  1:207870/1‑251 (MQ=255)
                                                                                                                                       ggcgCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCCCAAGCTGGTGCTGGCGGTGATCAAcccc                                                                                                                      >  1:212281/1‑250 (MQ=255)
                                                                                                                                               ggCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCTGCAGGa                                                                                                              >  1:361057/1‑250 (MQ=255)
                                                                                                                                                 cACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAAcc                                                                                                           >  1:136902/1‑250 (MQ=255)
                                                                                                                                                  ccctttttAACACAGTTCATTATTGTTGCGCCACACCAGCCCCCAGGGCTTTAATAAGGCAAATTGCGAACCCAACTGTTCATTCCCCCCCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTc                                                                                                          <  1:54659/246‑1 (MQ=255)
                                                                                                                                                          taaacaagggTCATTATGGTTGCCCCAAACAAGGCCCCAAGGGTTAAAAAAGGAAAATTGAGGACCCCACTGTTCCTCCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCAcc                                                                                                       <  2:332898/240‑1 (MQ=255)
                                                                                                                                                                    gtatttaggttggcgcCAAACCAGGCCCCAGGGCTGAAAAAAGGAAAAGTGGGCATCCCCCGGTTCATCCCCCCCCCGGCCGTGGAAACGAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGc                                                                                           <  2:134994/239‑1 (MQ=255)
                                                                                                                                                                     attttttgttgCGCTAACCCACGCGCAAAGGCTGTAATAAGGCAAATGGCGGACCCCCCTGTTCATCCCCCACCGCGGGGGGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCa                                                                                          <  2:155190/245‑1 (MQ=255)
                                                                                                                                                                      attattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTTTTCTTCCCCCCCCGCGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCAt                                                                                         <  1:361185/250‑1 (MQ=255)
                                                                                                                                                                       ttattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCCTCCCGTGCAGGAACTCACCGGGGGAATCGGCATa                                                                                        >  1:282191/1‑250 (MQ=255)
                                                                                                                                                                        tattGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATTCCGGGCAAGGCACCACCGGCGGGATCGGCATaa                                                                                       >  2:25822/1‑250 (MQ=255)
                                                                                                                                                                                  gcCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGGGGGATCGGCATAAGCGGCATTgg                                                                             >  1:278031/1‑251 (MQ=255)
                                                                                                                                                                                   caaaaCCAGGCGACAAGGCTATAAAAAGGAAAATTGCGTATCCCAATTTTCAGTCCCCCCCGTGGCGGGTAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGa                                                                            <  1:307805/249‑1 (MQ=255)
                                                                                                                                                                                    cAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGAAAATCCCGTGCACGAAATCACCGGGGGAAACGGGATaaaggacatttggg                                                                           >  2:341871/1‑239 (MQ=255)
                                                                                                                                                                                                ccAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAAATCCCCCGGGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTa                                                                 >  1:397692/1‑250 (MQ=255)
                                                                                                                                                                                                      cTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCCTGCAGGAACTCAACGGTGGAATCGGCATAAGCGGCATTGGAGGCATTTCTTAAATcgcg                                                          >  1:409543/1‑251 (MQ=255)
                                                                                                                                                                                                        ccgatcgggCAAATGGCGGACCACACTGTTCATCCCCCCCCCCGGCGGGGGAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCt                                                        <  2:173059/244‑1 (MQ=255)
                                                                                                                                                                                                               ggCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGGGCAGGAAATCACCGGTGGGATCGGCAAAAGCGGCATTGGAGGCATTTCTTAAATCGCGCt                                                        >  1:173059/1‑244 (MQ=255)
                                                                                                                                                                                                                   aaaTTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGGAGGAACTCACCGGGGGAATCCGCATAAGCAGCATTAGAGGCATTGCTTAAATCGCGCTCAGGTgccc                                               >  1:58139/1‑249 (MQ=255)
                                                                                                                                                                                                                                              ccccACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCCCCGGTGGAAACGGCATAAGAGGCATTGGAGGGATTGCTTAAATCGCGCTCACGTGCCCGCCACTGGTTCAgcgaccataacgc                      >  1:264228/1‑241 (MQ=255)
                                                                                                                                                                                                                                                    cccggcagTGGGAACGAAATAACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCATGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCGGGTGCCAGCCACTGGTTCAGCGCCCTGAACgccagcc                 <  1:134774/244‑1 (MQ=255)
                                                                                                                                                                                                                                                             gggAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATAATCCCGTGCAGGAACTCACCGGGGGAATCGGCCGAAACGGGATTGGGAGCATTGCTTAAAACGCGCTCAGGGGCCAGCCAATGGTTCAGCGCACCTAACGCCAGACACttgaccc        >  1:338445/1‑245 (MQ=255)
                                                                                                                                                                                                                                                                   gAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGGGATCATCCCGGGCAGGAACACACAGGGGGGAACGGCATAAGCGGCATTGGGGGCATTGCTTAAATCGCGCTCAGGTGGCAGCCACAGGTTCAGCGAACTTAAAgccccgatcggaagagcggcgg  >  2:159396/1‑232 (MQ=255)
                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                           
TACCGAATCGTTGCTCGTATTGATAAAACAGCTTCATGATCACCACAATCGATACCAGCGACAACAGCAGCGAGATAGTAAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATCATCCCGTGCAGGAACTCACCGGTGGAATCGGCATAAGCGGCATTGGAGGCATTGCTTAAATCGCGCTCAGGTGCCAGCCACTGGTTCAGCGCACTTAACGCCAGCCACTGCACCACGCCCC  >  NC_000913/1174376‑1174876

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: