Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC43.7% 147.6 / 20.4 88V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (44/5);  new base C (0/38);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.45e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378653‑4378951
                                                                                                                                                           |                                                                                                                                                      
ccccTCGAGCAATAGCGTCGCCACCAGGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑a                                                                                                                                                        >  2:934501/1‑150 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑a                                                                                                                                                        <  2:910619/150‑1 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc‑g                                                                                                                                                    <  1:1180369/151‑1 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc‑g                                                                                                                                                    >  1:1293238/1‑151 (MQ=255)
      gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc‑gccac                                                                                                                                                >  2:117394/1‑151 (MQ=255)
         cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc‑gccaccgc                                                                                                                                             <  1:424296/151‑1 (MQ=255)
         cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc‑gccaccgc                                                                                                                                             >  1:1320931/1‑151 (MQ=255)
                  cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCgg                                                                                                                                    >  1:248766/1‑151 (MQ=255)
                   gCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCggg                                                                                                                                   >  2:1076031/1‑151 (MQ=255)
                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGAt                                                                                                                                 >  1:750403/1‑151 (MQ=255)
                           gCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTttg                                                                                                                           >  2:1260669/1‑151 (MQ=255)
                            cACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTttgt                                                                                                                          >  2:752758/1‑151 (MQ=255)
                              cGCCAGAAATTGTTTTGCTGACCGAACC‑TACCTCAAACACCTTTTGTTGTGGGACGGGCTGTTTTTGGGCGATGTCCGCATACCCCCAGGTAAAGTAATAAGT‑TTTACCC‑TG‑AAAAA‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTttg                                                                                                                           <  2:346453/148‑1 (MQ=255)
                                ccATTAAATTTTTTGCTGCCCGACCC‑TACCTAAAACAGCTTTTGCTGTGTGCGGGGCTGTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGTAA‑AAGG‑TTTTCCCTTA‑ATAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                      <  2:196966/151‑1 (MQ=255)
                                ccAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                      >  1:516425/1‑151 (MQ=255)
                                   gTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                    >  2:922752/1‑150 (MQ=255)
                                     aatttttttGCTGCCCAACCC‑TAACTAAAACACCTTTTGCTGTGTGGCGGGCTGTTTTTGGGCGTTGCCCGCAT‑GCCCCCGGGAAAAAAAAAAGG‑TTTTCCC‑CT‑AAAAA‑TTCCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                   <  2:1201896/144‑1 (MQ=255)
                                     aattGTTTTGCTCACCGAACC‑AAACACAAACAGCTTTTGCTGTGTGGCGGGTTGTTTTTTGGGGATCCCCCCATCCCCCCAGAAAAAATAATAAGT‑TTTCCCT‑TA‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                    <  1:376252/145‑1 (MQ=255)
                                     aaaTTTTTTGCTGCCCAAACT‑TACCACAAACATCTTTTGTTGTTGGGCGGTTTTTTTTTTGGGGATGCCCGCATCCCCCCAGGTAAAGTAATAAGT‑TTTCCCT‑TG‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                   <  1:1149160/149‑1 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                       >  1:1279921/1‑145 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                  >  2:321888/1‑150 (MQ=255)
                                       aTGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                 >  1:14636/1‑149 (MQ=255)
                                          ttttGCGGCCCAACCC‑AAACTAAAACAACTTTTGCTGTGTGAGGGGCTTTTTTTGGGCGATGTCCGCATACCCCCAGGAAAAGTAATAGGT‑TTTCCCC‑TG‑AAAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                              <  2:12043/149‑1 (MQ=255)
                                           tttGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                           >  2:72600/1‑151 (MQ=255)
                                                tGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                      >  2:1076358/1‑151 (MQ=255)
                                                tGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                      >  2:880238/1‑151 (MQ=255)
                                                cgACCCAACCGTACCTCAAAAACCTTTTGTTGGGTGGGGGGCTTTTTTTGGGCGATTCCCGCACCCCCCAGGGTAAAGTAA‑AATG‑TTTTCCC‑CT‑AAAAA‑TT‑CCCGGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        <  1:492706/149‑1 (MQ=255)
                                                   ccGACCC‑AAACTAAAACAACGTTTGCGGTGTGAGGGTCTGTTTTTGGGCGATGCCCTCATGCCCCAAGAAAAAGTAATAAGT‑TTTCCCC‑TA‑ATAAT‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                        <  1:970622/146‑1 (MQ=255)
                                                   ccGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                   >  1:1016186/1‑151 (MQ=255)
                                                    cGACCC‑AAACTCAAACAACTTTTGTTGTGGGCGGGGTTTTTTTTGGGCGATGCCCGCATACCCCCAGAAAAAGTAA‑AAGG‑TTTTCCCATG‑AAAAA‑TC‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                       <  2:412064/146‑1 (MQ=255)
                                                      accc‑TAACAAAAACAACGTTGTGTGTGTGACGGGTTGTTTTTGGGTTATGCCAGTATCCCCCATGGTAAAGAAAAAATG‑TTTCCCC‑AA‑ATAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                                <  1:694024/149‑1 (MQ=255)
                                                           tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                            >  1:1154054/1‑151 (MQ=255)
                                                            acccaaaaCACCTTTTGTTTTGGGGCGGGCTGTTTTTGGGTGATGTCCGCAT‑GCCCCAAGGAAAATAAAAAGG‑TTTTCCCTTA‑AAAAA‑TT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                            <  2:161490/145‑1 (MQ=255)
                                                                  accactgTTTGCTGTGGGCGGGGCTTTTTTTTGGCGATCCCCGAATACCCCCAGAAAAATAAATAATT‑TTTCCCT‑TA‑ATAAA‑TT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     <  1:622395/145‑1 (MQ=255)
                                                                  accaacGTTTGTTGTGTGAGGGTTTGTTTTTTGCCAATGCCCGAATCCCCCTAGGAAAAGTAATAGGT‑TTTCCCT‑TA‑AAAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     <  2:228860/149‑1 (MQ=255)
                                                                  aacatcTTTTGCTGGGTGGCGGGCTTTTTTGTGGCGATCCCCAAAAAGCCCCTGAAAAAGTAATAAGT‑TTCCCCC‑TG‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                     <  1:350072/151‑1 (MQ=255)
                                                                   acaagtttggCTGTGTGACGGCTTGTTTTTTGGGAATGCCCGAATCGCCCCAGGAAAATTAA‑AAGG‑TTTTCCCTTG‑AAAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                         <  1:586222/137‑1 (MQ=255)
                                                                    caacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                   >  1:899669/1‑151 (MQ=255)
                                                                     aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                  >  1:659307/1‑151 (MQ=255)
                                                                       cGTTTGCTGTGTGACGGGCTGCTTTTTGGGGATGTCCGCATAGACCCAGGCAAAGTACTAAGG‑CTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                 >  1:1156246/1‑150 (MQ=255)
                                                                        gTTTGGTGTGTGGCGGTCTGTTTTTTGGCTTTGCCCGCATGGCCCGGGAAAAAGTAATATGG‑TTTCCCC‑TA‑ATAAA‑TC‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                 <  2:631059/149‑1 (MQ=255)
                                                                          ttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                             >  2:640886/1‑151 (MQ=255)
                                                                          ttGCTGTGAGGCGGGCTTTTTTTTGGGAATGCCCGAATAGCCCCTGTTAAAGTAATAAGG‑TTTCCCG‑TG‑AAAAT‑TG‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                              <  2:1273314/150‑1 (MQ=255)
                                                                          tgggtggggggcGGGCTCTTTTTTGGCGATCCCCGCTGACCCCCAGAAAATGTAATGTGG‑TTTCCCT‑TG‑AAAAA‑TC‑GCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                <  1:496982/137‑1 (MQ=255)
                                                                           tggtgtgtgAGGGGCTTCTTTTTGGCGATGCCCGAATACCCCAGGGTAAATAAATAGGG‑TTTCCCC‑TA‑ATAAT‑TT‑ACC‑GCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                              <  2:190504/146‑1 (MQ=255)
                                                                           tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                            >  2:328761/1‑151 (MQ=255)
                                                                              tgtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                         >  2:1053273/1‑151 (MQ=255)
                                                                               gtgtgACGGGTTTTTTTTGGGCGATGTCCGAATAGCCCCAGAAAAAGTAAAAGGT‑TTTCCCC‑TA‑AAAAT‑TT‑CCC‑GCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                              <  1:1245198/145‑1 (MQ=255)
                                                                                   ggcggttttttttttGGCGATGCCCGCATACCCCCAGGAAAAGAAATAGTG‑TTTCCCC‑TA‑AAAAT‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                    <  1:617205/141‑1 (MQ=255)
                                                                                   gcgggttttttttttGGTTATCTCCCCATACCCCCTGAAAAGATAATATTT‑TTTCCCT‑TG‑ATAAA‑TT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                      <  1:634045/139‑1 (MQ=255)
                                                                                   gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                    >  2:947691/1‑151 (MQ=255)
                                                                                    aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                   >  2:750523/1‑151 (MQ=255)
                                                                                     cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  1:1006924/1‑151 (MQ=255)
                                                                                      ggtcttttttttGGCGATGTCCGCAGACCCCCAGATAAAGTAATAAGT‑TTTACCT‑TA‑AAAAT‑TT‑CCC‑GCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                    <  1:293111/141‑1 (MQ=255)
                                                                                      ggtcttttttttGGCGATGCCCGCATACCCCCAGTAAAAGAAATAAGG‑TTTCCCC‑TG‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                           <  1:695254/134‑1 (MQ=255)
                                                                                      gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  2:1278652/1‑150 (MQ=255)
                                                                                      gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  2:1130516/1‑150 (MQ=255)
                                                                                           gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                             >  2:1252255/1‑150 (MQ=255)
                                                                                             ttttgGGCAAGGTCCTCAGAGCCCCAGTAAAACTAAAATGT‑TTTACCT‑TA‑AAAAT‑TT‑CCC‑GCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                           <  1:757512/150‑1 (MQ=255)
                                                                                                 ggggaatgccccaaaaGCCCCAGG‑AAAAAAAAAAAG‑TTTTCCC‑TT‑AAAAA‑TTCCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       <  1:66665/135‑1 (MQ=255)
                                                                                                  ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTCCCC‑TG‑ATAAA‑TT‑ACC‑GCCACCCCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                      <  2:930810/150‑1 (MQ=255)
                                                                                                   gCGATGTCCGCATACCCCAAGGTAAA‑TAAAAAGGTTTTCCCC‑TG‑AAAAA‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                    <  2:1190235/151‑1 (MQ=255)
                                                                                                   gCGATGCCCGCATCCCCCCGGGAAAAGTAATAA‑G‑TTTTCCC‑TGAATAAA‑TC‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                    <  2:1331774/151‑1 (MQ=255)
                                                                                                        gTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTAAGGGGCAGCAAATGTGGAGCAag                                                 >  1:288610/1‑149 (MQ=255)
                                                                                                          ccgaatgccccgAGGAAAATTAAAAAGT‑TTCCCCC‑TA‑AAAAA‑TT‑CCC‑GCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                             <  2:470773/143‑1 (MQ=255)
                                                                                                          ccGAATACCCCCAGTTAAAATAAT‑AGG‑TTTACCC‑CT‑ATAAATTT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                <  2:957610/148‑1 (MQ=255)
                                                                                                           cGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                            >  1:1173808/1‑151 (MQ=255)
                                                                                                              aTAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          >  2:321983/1‑150 (MQ=255)
                                                                                                              aTAGCCCCAGGAAAAGTAATAGGT‑TTTCCCC‑TA‑ATAAA‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          <  2:1090649/150‑1 (MQ=255)
                                                                                                                      aGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                  >  2:219392/1‑150 (MQ=255)
                                                                                                                       ggaaaaGAAATAGGG‑TTTCCCC‑GA‑AAAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                <  1:52462/148‑1 (MQ=255)
                                                                                                                       ggTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                  >  1:166236/1‑149 (MQ=255)
                                                                                                                          aaaGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                              >  1:816196/1‑150 (MQ=255)
                                                                                                                           aaattaaaaGG‑TTTTCCCTTG‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          <  1:916834/131‑1 (MQ=255)
                                                                                                                               aatagGT‑TTTCCCC‑TA‑AAAAA‑TC‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                        <  1:1168906/151‑1 (MQ=255)
                                                                                                                               aatagGT‑TTTACCC‑TA‑AAAAA‑TC‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                        <  1:433620/151‑1 (MQ=255)
                                                                                                                               aataaGT‑TTTCCCC‑TA‑ATAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                           <  1:1330040/148‑1 (MQ=255)
                                                                                                                               aataaGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                          >  1:1141658/1‑149 (MQ=255)
                                                                                                                                    gg‑TTTCCCT‑TA‑AAAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                    <  2:777868/150‑1 (MQ=255)
                                                                                                                                    gg‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAgcggcg                   >  1:1096069/1‑151 (MQ=255)
                                                                                                                                       tttACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                 >  2:523438/1‑151 (MQ=255)
                                                                                                                                        ttACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGctt                >  1:660230/1‑151 (MQ=255)
                                                                                                                                         tttcccttaaaaaa‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                    <  2:221886/137‑1 (MQ=255)
                                                                                                                                           ccc‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg              >  2:855045/1‑150 (MQ=255)
                                                                                                                                           ccc‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGCGGTAATTGTGCGATGCACAATATCGTTGCTTTGTTGAGGGGCAGCAAATGTGGCGCAAGAGGCGGTAATTAATAATGCGCAAAGGGCCGtt                 >  1:1061390/1‑147 (MQ=255)
                                                                                                                                           ccc‑TA‑ATAAA‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg              <  2:807464/150‑1 (MQ=255)
                                                                                                                                            cc‑TT‑ATAAT‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg              <  2:209154/149‑1 (MQ=255)
                                                                                                                                                  aaaaa‑TT‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa            <  2:841327/145‑1 (MQ=255)
                                                                                                                                                  aaaaa‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa          <  2:248800/147‑1 (MQ=255)
                                                                                                                                                  aTAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt         >  2:938964/1‑150 (MQ=255)
                                                                                                                                                  aTAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt         >  2:240812/1‑150 (MQ=255)
                                                                                                                                                   tAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt         >  2:1091845/1‑149 (MQ=255)
                                                                                                                                                    aaa‑TT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAg       <  2:353426/150‑1 (MQ=255)
                                                                                                                                                    aaa‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg      >  2:237660/1‑151 (MQ=255)
                                                                                                                                                         t‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt  >  2:1216137/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378653‑4378951

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