Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC46.9% 158.0 / 41.7 113V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (55/5);  new base C (0/53);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 9.73e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGG  >  NC_000913/4378653‑4378947
                                                                                                                                                         |                                                                                                                                                   
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa                                                                                                                                                     <  1:1250421/150‑1 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccg                                                                                                                                                 <  2:808455/151‑1 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccg                                                                                                                                                 <  2:1474092/151‑1 (MQ=255)
  aatCGAGCAATAGAGTAGCCACCAAGCACGCCAGTAAATGTTTTGATGACAGAAAC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGACCAAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ctg                                                                                                                                                 >  2:667150/3‑149 (MQ=255)
   cTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccgc                                                                                                                                                >  1:1362210/1‑151 (MQ=255)
      gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccgccac                                                                                                                                             >  2:1221803/1‑151 (MQ=255)
        gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccgccacc                                                                                                                                            <  1:1218071/150‑1 (MQ=255)
         cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTa‑ccgccaccgc                                                                                                                                          <  2:1191163/151‑1 (MQ=255)
               cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATAcc                                                                                                                                    >  2:1180957/1‑151 (MQ=255)
                   gCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCggg                                                                                                                                >  2:392140/1‑151 (MQ=255)
                       ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCt                                                                                                                            >  2:269923/1‑151 (MQ=255)
                         aaGCCCGCAAGTAATTGTTTTGC‑GACCCAACCGAAACTAAAAAAATTTTTGCTGGGTGAGGGGTTTTTTTTTGGCAATGTCCGCATAGCCCCAGAA‑AAATTAAAAAGTTTTCCCT‑TG‑AAAAATCC‑CCCCCACCGCCATACCCGGGATCttt                                                                                                                          <  2:2089956/151‑1 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  2:449121/1‑151 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  2:1391735/1‑151 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  2:443089/1‑151 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  1:901206/1‑151 (MQ=255)
                           gCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttg                                                                                                                        >  1:978046/1‑151 (MQ=255)
                             aCGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  1:88481/1‑151 (MQ=255)
                             aCGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  1:593502/1‑151 (MQ=255)
                             aCGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  2:560904/1‑151 (MQ=255)
                              cGCCAGTAATTTTTGTGCGGACCAACCC‑TAACAAAAACATCTTTTGTTGTGTGAGGGGTTGTTTTTTGGCGATGCCCACATCCCCCCGGAA‑AAAAAAATAAGTTTTCCCA‑TA‑AAAAATCC‑CCCCCCCCGCCATACCCGGGATCtttt                                                                                                                         <  2:1564958/147‑1 (MQ=255)
                              cGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                     >  2:1549301/1‑151 (MQ=255)
                              cGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                     >  2:1754140/1‑151 (MQ=255)
                                ccAGTAAGTGTTTGGCTGACCGCCCC‑CACCAAAAACAGGTTTGGCTGGGGTCCGGGCTGTTTTTTGGCGAGGTCC‑TAACGGCCCAGGTAAAAATAATAGTTTTTCCCG‑AG‑ATAATTCC‑CCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      <  1:913940/148‑1 (MQ=255)
                                ccAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                   >  2:1106358/1‑151 (MQ=255)
                                 cAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  >  2:1227864/1‑151 (MQ=255)
                                    tAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               >  2:795958/1‑151 (MQ=255)
                                     atttGTTTGGCTCACCAACCC‑TACCCAAAACACCTTTTGCTGTGGGACGGGTTGTTTTTTGGCGATGCCCGAATACCCCCAGAA‑AAAGTAATAAGTTTTCCCT‑TA‑ATAATTTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               <  1:770893/147‑1 (MQ=255)
                                       aTGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                              >  2:516199/1‑149 (MQ=255)
                                          ttttGCTGCCCGAACC‑AAACTCAAACACCGTTTGCTGTGGGGGGGGCTGTTTTTTGGGGATTCCCGAATACCCCCGGAT‑AAAAAAAAAAGTTTTCCCT‑TA‑ATAATTAC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                           <  1:198604/149‑1 (MQ=255)
                                          ttttGCTGCCCAACCC‑TACCAAAAACACCTTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCATACCCCCAGGA‑AAAGAAATAAGTTTTCCCC‑TA‑AAAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                          <  1:1466637/150‑1 (MQ=255)
                                          ttttGCTGACCGAACC‑TAACTCAAACACGTTTTGTTGTGTGAGGGGTTGTTTTTTGGCGATGTCCGCATAGCCCAAGGA‑AAAGTAATAAGTTTTCCCT‑TA‑ATAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                          <  1:82277/150‑1 (MQ=255)
                                            ttGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                       >  1:1300533/1‑151 (MQ=255)
                                              gCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                      >  2:181337/1‑150 (MQ=255)
                                                tcaCCGACCT‑TAACAAAAACAATGTTTGCTGTGTGCCGGGCTTTTTTTTGGCCATCCCCCCAGACCCCAAGAT‑AAAAAAAAAAGGTTTCCCC‑GA‑AAAAATCC‑CCCCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                       <  1:1640219/145‑1 (MQ=255)
                                                tGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                   >  1:676335/1‑151 (MQ=255)
                                                 gACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                  >  2:2160021/1‑151 (MQ=255)
                                                  acataaCC‑TAACTCAAACAACTTTTGTTGTGAGACGGTCTTTTTTTTGGCTTTCCCCACAGACCCCGATGT‑AACATAAAAATGTTTACCT‑TG‑ATAATTCC‑ACGCAGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAttt                                                                                                 <  2:1546174/147‑2 (MQ=255)
                                                  acagAACC‑TAACTAAAACAACTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATGCCCGCACAGCCCCAGAA‑AAAGTAATATGTTTTCCCT‑TA‑AAAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                   <  2:122623/146‑1 (MQ=255)
                                                   ccaaccc‑TAACTAAAACACCGTTTGTTGTGTGAGGGTCTTTTTTTTGGCGTTTCCCGCATACCCCCAGGA‑AAAGTAATAAGTTTTCCCT‑TA‑ATAATTTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                    <  1:726020/142‑1 (MQ=255)
                                                      aaCC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  2:1808385/1‑151 (MQ=255)
                                                      aaCC‑TAACAAAAACAACTTTGGGTGTGGGAGGGGTTTTTTTTGGGCGATGCCCACACACCCCCGGAA‑AAAGTAAT‑AGGTTTTCCC‑TT‑AATAATTACCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                              <  2:662893/150‑1 (MQ=255)
                                                        cc‑TAACCAAAACAACGTTTGTTGTGTGGCGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGTA‑AAAGTAATAAGTTTTCCCT‑TG‑ATAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                               <  1:1362206/147‑1 (MQ=255)
                                                         c‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                           >  1:1895561/1‑150 (MQ=255)
                                                            accaaaaaCAACTTTTGCTTTGTGAGGGGTTTTTTTTTGGCTATGCCCGCATACCCCAAGAT‑AAAGTAATAGGTTTTCCCC‑AG‑ATAAATTC‑CCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                           <  1:1028686/143‑1 (MQ=255)
                                                            aacaaaaaCAACTTTTGCTGGGTGACGGGCTTTTTTTGGGCGATGCCCGCATACCCCCAGGA‑AAAGAAATAAGTTTTCCCT‑TA‑ATAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                        <  1:357121/146‑1 (MQ=255)
                                                             aCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       >  1:6557/1‑151 (MQ=255)
                                                              cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       >  1:2152846/1‑150 (MQ=255)
                                                                 aaaCACCTTTTGCTGTGTGGGGGGCTGTTTTTTGGGGTTGCCCGAATACCCCCAGAT‑AAAG‑AAAAAGGTTTTCCCTTA‑ATAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       <  1:1101247/147‑1 (MQ=255)
                                                                 aaaCACCTTTTGCTGTGGGACGGGTTTTTTTTTGGCGATGTCCCCATCCCCCCAGAA‑AAAGTAATAAGTTTTCCCT‑TA‑ATAATTCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                          <  2:593426/144‑1 (MQ=255)
                                                                 aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    >  2:2035928/1‑150 (MQ=255)
                                                                         tttGTTGTGTGGGGGGTTTTTTTTGGGGGATGTCCCCATCCCCCCAGAA‑AAAAAAAAAGGGTTTCCTT‑TA‑ATAAATGC‑CCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGa                                                                            <  2:1599690/150‑1 (MQ=255)
                                                                         tttGCTGTGTGACGGTCTGTTTTTTGGCGATGCCCGCATAGCCCGAGAT‑AAAGTAATAAGGTTTCCCC‑TG‑ATAATTTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             <  1:1665606/149‑1 (MQ=255)
                                                                          ttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                          >  1:2083983/1‑151 (MQ=255)
                                                                           tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                         >  2:1521349/1‑151 (MQ=255)
                                                                            gttgggtgAGGGGTTGTTTTTGGGCGTTGTCCGCATACCCCCGGGA‑AAAGAAAAAGGTTTTACCG‑TA‑ATAATTTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                         <  2:595122/145‑1 (MQ=255)
                                                                              tgtgtgGGGGGCTGTTTTTTGGGGATGTCCGCATACCCCCAGGA‑AAAGTAAAAAGTTTTCCCT‑TG‑AAAATTTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                          <  1:1011886/147‑1 (MQ=255)
                                                                                tgtgAGGGGTTTTTTTTTGGCGATCCCCGCATCCCCCCGGTA‑AAAGTAAAAAGTTTTCCCC‑TA‑AAAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  2:1035771/147‑1 (MQ=255)
                                                                                 gggcggggTTGTTTTTTGGCGATGTCCGCAACCCCCCAGAT‑AAAGTAATATGTTTTCCCT‑TA‑AAAATTTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                   <  2:972333/146‑1 (MQ=255)
                                                                                  ggACGGGTTTTTTTTTGGAAATGCCCCCACCCCCCAAGGA‑AAAGTAATGTTTTTCACCA‑TA‑TTAAATTA‑CCTCTCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     <  2:305401/147‑1 (MQ=255)
                                                                                    aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                >  2:732779/1‑151 (MQ=255)
                                                                                      gggCTTTTTTTTGGCGATGTCCGCATCGCCCCAGGT‑AAAGTAATAAGTTTCCCCC‑TA‑ATAATTTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                <  1:211237/149‑1 (MQ=255)
                                                                                      gggCTGTTTTTTGGCGATGTCCGCATACCCCCAGTT‑AAAAAAATAAGTTTTCCCC‑TG‑AAAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                   <  1:629962/146‑1 (MQ=255)
                                                                                      gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                               >  2:895373/1‑150 (MQ=255)
                                                                                        gctgctTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                            >  1:1145809/1‑151 (MQ=255)
                                                                                          tgtttttgGGCGATGCCCGCATAGCCCGAGGA‑AAAATAATAGGTTTTCCCC‑TA‑ATAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                             <  1:1838059/145‑1 (MQ=255)
                                                                                               tttGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                     >  2:2127755/1‑151 (MQ=255)
                                                                                                ttGGCGATCTCCGCATACCCCCAGAA‑AAAGTAATAAGTTTTCCCC‑TA‑ATAAATTC‑CCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                     <  1:2013085/150‑1 (MQ=255)
                                                                                                  ggCGATGTCCGCATACCCCCAGGA‑AAAGTAAAAAGTTTTCCCC‑TA‑ATAAATTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                     <  2:1479639/148‑1 (MQ=255)
                                                                                                   gCGATGCCCGCATACCCCAAGTA‑AAAGAAATAAGTTTTCCCC‑TA‑AAAATTCC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                 <  1:1411503/151‑1 (MQ=255)
                                                                                                        tcccGCACACCCCAAGGA‑AAAGAAATAGTTTTTCCCC‑AA‑ATAATTAC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  2:291417/147‑1 (MQ=255)
                                                                                                        gTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              >  2:164801/1‑149 (MQ=255)
                                                                                                          ccGCATACCCCCAGGA‑AAAGAAAAAAGGTTTCCCC‑TA‑AAAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  1:1955810/147‑1 (MQ=255)
                                                                                                              aTAGCACCAGTC‑AAAGTAATAAGGTTTACGC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGCCCGATGCACAATATCGTTGATTTGTTAAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGt                                      >  2:579581/1‑150 (MQ=255)
                                                                                                               tAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCgg                                       >  2:1024726/1‑148 (MQ=255)
                                                                                                                cgccccaagaaaaaGTAAAAATTTTTCCCG‑TG‑ATAATTTC‑CCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGta                                     <  2:1486965/139‑1 (MQ=255)
                                                                                                                  ccccAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGtaatt                                  >  2:1512026/1‑150 (MQ=255)
                                                                                                                  ccccAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGtaatt                                  >  1:983547/1‑150 (MQ=255)
                                                                                                                     cAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGtaattaat                               >  2:960755/1‑150 (MQ=255)
                                                                                                                        ggaaaataaaaaaGTTTTCCCT‑TA‑ATAAATTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATtaata                              <  1:1758202/138‑1 (MQ=255)
                                                                                                                         t‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAgg                           >  2:1955584/1‑150 (MQ=255)
                                                                                                                          aaaattatgatgttttCCCT‑TG‑ATAACTGC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGtaatt                                  <  2:2126318/130‑1 (MQ=255)
                                                                                                                           aaaaaaagaggttttCCCC‑TA‑AAAAATCC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGtaattaat                               <  2:2171086/133‑1 (MQ=255)
                                                                                                                           aaaaTAATAGGTTTTCCCC‑TG‑ATAAATTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAgg                           <  1:351791/149‑1 (MQ=255)
                                                                                                                           aaaaTAATAAGTTTTACCA‑TG‑AAAATTCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAg                            <  2:1978123/148‑1 (MQ=255)
                                                                                                                            aaGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGgcgc                        >  1:1522289/1‑151 (MQ=255)
                                                                                                                            aaGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGgcg                         >  2:885331/1‑150 (MQ=255)
                                                                                                                             aGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGgcgc                        >  1:332797/1‑150 (MQ=255)
                                                                                                                              gTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGgcgc                        >  2:445342/1‑149 (MQ=255)
                                                                                                                                aataaGTTTTCCCT‑TA‑ATAATTAC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCag                      <  2:1947033/149‑1 (MQ=255)
                                                                                                                                 ataaGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGc                   >  1:1140955/1‑151 (MQ=255)
                                                                                                                                  taaGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGcg                  >  1:1191262/1‑151 (MQ=255)
                                                                                                                                  ataagTTTTCCC‑TGAATAAATAC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGc                   <  2:195012/147‑1 (MQ=255)
                                                                                                                                     gttttCCCC‑TA‑AAAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGcgtc                <  1:139815/148‑1 (MQ=255)
                                                                                                                                      gTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGcgtcg               >  1:1257948/1‑150 (MQ=255)
                                                                                                                                        ttcccc‑TA‑AAAAATTC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGtt             <  1:1549875/147‑1 (MQ=255)
                                                                                                                                         ttcccttG‑ATAAATCC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGtt             <  1:2136315/144‑1 (MQ=255)
                                                                                                                                         tACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAGGCGGGGGAATTGTGCGATGCACAATATCGTTGATTTTTTGAGGGGCAGCAAAAGTGGGGCAAGAGGGGGGGAATTAA‑AAGGGGCAGAGCGTCGgttt           >  1:523587/1‑151 (MQ=255)
                                                                                                                                                    aaaTTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg  >  2:194552/1‑151 (MQ=255)
                                                                                                                                                      aTTC‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg  <  2:568135/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACC‑TAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑TG‑ATAAATTA‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGA‑GGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGG  >  NC_000913/4378653‑4378947

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: ‑