Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913119,3960CT24.1% 102.2 / 12.7 55S39F (TCT→TTT) ampEprotein AmpE
Reads supporting (aligned to +/- strand):  ref base C (10/31);  new base T (13/0);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.03e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGC  >  NC_000913/119149‑119634
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
aCTGATTGATTGCTATCCGGATATGGCGAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGAGCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGGTAATTGTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGACCATCGTCTGGAAGCATTCGTTCGTCGGGTGAAACATTTTtctct                                                                                                                                                                                                                                                 <  2:263431/251‑1 (MQ=255)
  tggttgattgCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTtctctcc                                                                                                                                                                                                                                               >  1:403161/4‑250 (MQ=255)
         ttgCTAACCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCCTCGTCTTGAAGCGTTCTTTCGGCGGGAGGAACAATTTTTTCTCGGtcgcgc                                                                                                                                                                                                                                        >  2:206569/1‑249 (MQ=255)
              aTCCGGATATCGCTAAAAACATGACGTGCCTTTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTTATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGGTAATTTTCGATCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTGCTTTCGGCTGGTGAAACATTTTTCTCTCGGGCGCACGTTa                                                                                                                                                                                                                                    <  1:541417/250‑1 (MQ=255)
                          cTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGGCCGTTAGGCCTGACCattt                                                                                                                                                                                                                      >  1:186745/1‑250 (MQ=255)
                            aaaaaCATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCatta                                                                                                                                                                                                                      <  1:442048/249‑1 (MQ=255)
                             aaaaCATTACGGGCCATTGTGATATTGCGCCGGATCGGAAAGCCGATCCCGGTCCGGCATTTTATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTGGAAGCGTTCTTTCGGCGGGTGAAACATGTTTCTCTCGGGCGCACGTCAGGCATGACCATTATTg                                                                                                                                                                                                                  <  2:37025/251‑1 (MQ=255)
                                 cATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCCTGACCATTATTGCGAt                                                                                                                                                                                                              >  1:123855/1‑251 (MQ=255)
                                     gcGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGc                                                                                                                                                                                                          <  2:55842/250‑1 (MQ=255)
                                       gggCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCCTTATTGCGAttggggg                                                                                                                                                                                                        >  1:445167/1‑245 (MQ=255)
                                               gtgATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTTGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACttttt                                                                                                                                                                                                <  1:87639/251‑1 (MQ=255)
                                                  atatTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATTACCATTATTGCGATGGGCGTGACTTTTTTAc                                                                                                                                                                                             <  1:298855/251‑1 (MQ=255)
                                                       gcgcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGtt                                                                                                                                                                                         <  1:162816/250‑1 (MQ=255)
                                                         gcCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACGTTTGGCATGCCCATTATTGCGATGGGCGTGGCTTTTTTACTGtttcg                                                                                                                                                                                      >  1:386715/1‑248 (MQ=255)
                                                            ggATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCCTTATTTCGTTGGGCGTGACTTTTTTACTTTTACccgc                                                                                                                                                                                  >  1:521127/1‑251 (MQ=255)
                                                               tCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTCTCGGGCGCACGCTATGCCTGACCCTTATTGCGATGGGCGTGGCTTTTTTACTGTTACGCGCAtt                                                                                                                                                                              >  1:14496/1‑251 (MQ=255)
                                                                cGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACACGCAtt                                                                                                                                                                              <  1:143520/250‑1 (MQ=255)
                                                                     aaCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGGTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAggg                                                                                                                                                                        <  1:237973/251‑1 (MQ=255)
                                                                      cccGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATTACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCTCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTTCAGGGa                                                                                                                                                                       <  2:317189/250‑1 (MQ=255)
                                                                      aCCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGa                                                                                                                                                                       <  2:140067/251‑1 (MQ=255)
                                                                                    tGCATTAGATTGGACACGGTTTCGTGTGCTGGTCAGCAAGGCGACAACATGACGCTATTTCCAACCTTACTGGGGTTAATTTTCGAGCGCCTGTTTAAGTTTGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGCCACCTTTTTTTCTCGGGGGCACGGTACGCATGACCCTTTTCTGGATGGGCGTGTCTTTTTTACTGTTACCCGCCTGTCAGGGAGTAtttgtcaccgt                                                                                                                                                         >  2:79645/1‑242 (MQ=255)
                                                                                        tttGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGAACACATGACGCTATTTACAACCTTACTGGCGTTAATTTTCGCGCGCCTGTTTACGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCAAGTTAGGCCTGACCATTATTGCGATGGGCGCGACTTTTTTAATGTTACGCGCATTGAAAGGAGGATTGTTCAACGTTccc                                                                                                                                                     >  1:32598/1‑251 (MQ=255)
                                                                                         ttGATTGGGCACTGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTccc                                                                                                                                                     <  2:345992/250‑1 (MQ=255)
                                                                                           gATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACg                                                                                                                                                  <  1:397321/251‑1 (MQ=255)
                                                                                                           gtgtGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGGTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGGAACATTTTTTTCTCCGGCGGCCGGTAGGCATGAACATTATTGCGACTGGCGTGACTTTTTTACTGATCAGCGCATTGCAGGGAGTATTGTTTAACGTTCCCACGCTCATGGTGTGgcagc                                                                                                                                  >  1:293634/1‑248 (MQ=255)
                                                                                                             gtgCTGGTCAGCAAGGAGACAAAATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGGTAGGCATGTCCATTTTTGCGGTGGGGGTGACTTTTTTTCTTTTACGCGGATTTCAGGGGGGATTGTTTTAAGTTTCCACGCAACTGGTGTGgccggtt                                                                                                                                >  1:147221/1‑246 (MQ=255)
                                                                                                              tgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATGACCcttatttc                                                                                                                                                                                                                 >  1:417333/1‑169 (MQ=255)
                                                                                                              tgCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGc                                                                                                                                                                                                                 <  2:417333/171‑1 (MQ=255)
                                                                                                                  ggTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTg                                                                                                                            <  1:15485/250‑1 (MQ=255)
                                                                                                                     cAGCAAGGAGCCAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTATTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCACGTTAGGCATTAACATTATTTCGATGGGGGTGACTTTTTTGCTGTTTCGCGCCTTGCAGGGAGTATTGTTTAACGTTCCCCCGCTACTGGTGTGGCTGCTTGTTGGttt                                                                                                                        >  2:66396/1‑251 (MQ=255)
                                                                                                                     cAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGttt                                                                                                                        <  1:390140/251‑1 (MQ=255)
                                                                                                                                caacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTg                                                                                                             <  2:46601/251‑1 (MQ=255)
                                                                                                                                 aacaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTTTCTCGGGCGCCCGTTAGGACTTAACCTTCTTGCGATTGGGGATACTTTTTTTTTTTTACGCGCATTTCAGGGAGGAATTTTTAACGTTCCCCCGCCAATGGGGTGGgtgctgtatgggttgtggttgtttgg                                                                                                            >  2:370767/1‑231 (MQ=255)
                                                                                                                                  acaTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgc                                                                                                           <  2:386715/251‑1 (MQ=255)
                                                                                                                                   caTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGgcg                                                                                                          <  2:403161/251‑1 (MQ=255)
                                                                                                                                   caTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGCGAAAAATTTTTTTCTTGGGCGCACGGTAAGCCTGACCCTTATTGCgc                                                                                                                                                                                                               >  1:425873/1‑151 (MQ=255)
                                                                                                                                   caTGACGCTATTTACAACCTGACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGa                                                                                                                                                                                                               <  2:425873/152‑1 (MQ=255)
                                                                                                                                    aTGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGGGGGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCCACTGGTGTGGCTGCTGGTTGGGTTGCTGTGTGTTGgcgc                                                                                                         >  1:357620/1‑251 (MQ=255)
                                                                                                                                      gACGCTATTTACAACCTTACTGGTGTTAAGTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCa                                                                                                        <  2:147221/250‑1 (MQ=255)
                                                                                                                                        cGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTTGTGTGGCGGCTGATTGGTTTGCTGTGTATTGGCGCAGGt                                                                                                     <  1:80332/251‑1 (MQ=255)
                                                                                                                                            aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTaa                                                                                                   <  1:476518/249‑1 (MQ=255)
                                                                                                                                            aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAg                                                                                                 <  1:482738/251‑1 (MQ=255)
                                                                                                                                            aTTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTaa                                                                                                   >  1:473680/1‑249 (MQ=255)
                                                                                                                                                  acccTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTGCTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTc                                                                                           <  1:380705/249‑1 (MQ=255)
                                                                                                                                                         cTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGGGGCACCTTTTTATCTCCGGGGCCCGGTAGGCCTGACCCTTTTTGCGCTGCGGGTGACTTTTTTACTGTTTACCGCCTTTCAGGGGGTATTGTTTACCGTTCCCCCTCTACTGGTGTGGCCGCTGATTGGTTTGTTGTGTATTTGCGCTGGTGAAAATCGTCTTccttt                                                                                    >  2:426595/1‑247 (MQ=255)
                                                                                                                                                             tgtTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTgg                                                                                                            <  1:386144/223‑1 (MQ=255)
                                                                                                                                                              gtTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATg                                                                               <  1:483639/251‑1 (MQ=255)
                                                                                                                                                                   ttttCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACCATTTTCTCTCTGGCGCACGATAGGCATGACCATTATTGCGATGGGCGTGGCTTTTTTACTGTTACGCGCCTTGCAGGGAGGATTGTTCAACTTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTTCTGTGTTTTGGCGCCGGTAAAGTTCGTCTTTATTTTCATGCTTAt                                                                          >  1:395185/1‑251 (MQ=255)
                                                                                                                                                                      tCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCt                                                                        <  2:158232/250‑1 (MQ=255)
                                                                                                                                                                       cGAGCGCCTGTTTAGTGTGGGCGAGCGCTGGCAGCTTGATCATCGTCGTGAAGCGTGCGTGCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGGTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGGTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTg                                                                       <  2:357620/250‑1 (MQ=255)
                                                                                                                                                                           cgcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACCTTTTTCTCTCGGGCGCACGGTTGGCATGACCATTATTGCGGTGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAAAGTTCCCACGCTACTGGTGTGGCTGCTGGTTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAgc                                                                  >  1:409197/1‑251 (MQ=255)
                                                                                                                                                                           cgcCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTTTTTCGGGCGCACCGTAGGCATGCCCATTATTGCGCTGGGCGTGACTTTTTTTCTGTTACCGCCATTTCAGGGGGGATTTTTCAACGTTTCCACGCTACTTGTGTGGCTGCTGATTGGTTTGCTGGGTTTTGGCGCAGGGAAAGGTCCTCTTCATTAttatttttatttggccgc                                                                  >  2:235833/1‑234 (MQ=255)
                                                                                                                                                                                gTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCtt                                                             <  2:123855/251‑1 (MQ=255)
                                                                                                                                                                                          gggCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTTTTGCAATGGGGGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGGTTGCTGGGGATTGGCGCAGGGAAAGTTCGTTTTCATTATTATGCTTATCTGACAGCTGCTTAACGTAAtt                                                    >  1:397813/1‑250 (MQ=255)
                                                                                                                                                                                                                    gTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTAATTGGt                                                                                                                          <  1:149463/155‑1 (MQ=255)
                                                                                                                                                                                                                        tGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGGGTATTTTTTAACGTTCCCACGCTACTGGGGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGGAAAGTTCGTCTTCATTATCCTGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCg                      >  1:404138/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTCTTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGc  >  1:525406/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
ACTGATTGATTGCTATCCGGATATCGCTAAAAACATGACGGGCCATTGTGATATTGCGCCGGATCGGAAAACCGATCCCGGTCCTGCATTTGATTGGGCACGGTTTCGTGTGCTGGTCAGCAAGGAGACAACATGACGCTATTTACAACCTTACTGGTGTTAATTTTCGAGCGCCTGTTTAAGTTGGGCGAGCACTGGCAGCTTGATCATCGTCTTGAAGCGTTCTTTCGGCGGGTGAAACATTTTTCTCTCGGGCGCACGTTAGGCATGACCATTATTGCGATGGGCGTGACTTTTTTACTGTTACGCGCATTGCAGGGAGTATTGTTCAACGTTCCCACGCTACTGGTGTGGCTGCTGATTGGTTTGCTGTGTATTGGCGCAGGTAAAGTTCGTCTTCATTATCATGCTTATCTGACAGCTGCTTCACGTAATGATAGCCATGCCCGTGCCACGATGGCTGGCGAACTCACCATGATTCACGGC  >  NC_000913/119149‑119634

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