Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC42.6% 118.3 / 19.0 68V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (37/2);  new base C (0/29);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 3.38e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378658‑4378952
                                                                                                                                                       |                                                                                                                                                       
cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgcc                                                                                                                                                   <  2:513656/150‑1 (MQ=255)
cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgc                                                                                                                                                    <  1:480145/149‑1 (MQ=255)
 gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccac                                                                                                                                                 >  1:664239/1‑151 (MQ=255)
 gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccac                                                                                                                                                 >  1:1248232/1‑151 (MQ=255)
   gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccg                                                                                                                                               >  1:470753/1‑151 (MQ=255)
    cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccgc                                                                                                                                              >  2:1421082/1‑151 (MQ=255)
        aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCAt                                                                                                                                           >  2:918562/1‑150 (MQ=255)
        aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCAt                                                                                                                                           >  1:493411/1‑150 (MQ=255)
                  ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCt                                                                                                                                >  2:1325856/1‑151 (MQ=255)
                  ccaACCCCGCCAGAAATTGTTTTGCTCCCCAATCCAAACTCAAACAAGTTTTGCTGTGTGAGGGGTTGTTTTTGGGCGATGCCCGCATACCCCAAGAAAAAGAAATATGTTTTCCCGTAATAAATTCCCGCCACCGCCATACCCGGGATc                                                                                                                                 <  1:305522/150‑1 (MQ=255)
                       cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACACGGGATCTTttgt                                                                                                                           >  2:694972/1‑151 (MQ=255)
                        aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          >  2:475731/1‑151 (MQ=255)
                         cGCCAGTAAATTTTTTGCTGACCGAACCTAACTAAAACAACGTTTGCTGTGTGAGGGGCTGTTTTTGGGCGATGCCCGCATACCCCCAGTAAAAGTAATAGGGTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                         <  1:704657/151‑1 (MQ=255)
                          gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        >  2:485619/1‑151 (MQ=255)
                           ccATTAAGTGTTTGGTGGCCCAAACCAAACCCAAAAATGGTTTGTTGTGGGAGGGTCTTTTTTTGGGCGATGCCCCCCCAGCCCAAGGAAAAATAAAAGGGTTTCCCAAGAAAATTCCCCGCCACCGCCATACCCGGGATCt                                                                                                                                <  2:268391/142‑1 (MQ=255)
                           caattAATTGTTTTGCTGCCCGACCCTACCTAAAACATCTTTTGCTGTGTGGGGGGTTTTTTTTGGGCGATGTCCGCATACCCCAGGGAAAAGAAATAGGGTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                       <  1:1242584/147‑1 (MQ=255)
                             aGTAATTTTTTGGCTGACCGACCCTACCTAAAACACCTTTTGCTGGGTGGCGGGCTGTTTTTTGGCGATGTCCGCATACCCCAGGTAAAAGTAAAAAGTTTTCCCCTGAAAATTCCCCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          <  1:481326/146‑1 (MQ=255)
                                aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:675946/1‑150 (MQ=255)
                                aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  1:848505/1‑150 (MQ=255)
                                 aaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  >  1:813280/1‑150 (MQ=255)
                                 aaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  >  2:645598/1‑150 (MQ=255)
                                     ttttGCGGCCCAAACCTACCACAAACAACGTTTGTTGTGTGCCGGGCTTTTTTTGGGCGATCCCCGAATAGCCCAAGTAAAAGTAATAAGTTTTCCCTTAAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                              <  1:1236449/150‑1 (MQ=255)
                                     tttgGCTGACCAAACCAACCCCAACAAACGTTTGTTGTGTGACGGGCTTTTTTTTGGCGATGCCCCCACACCCCCAGAAAAAGTAATAAGTTTTCCCTTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                               <  2:36369/149‑1 (MQ=255)
                                          cTGCCCGAACCTACCTAAAACAACGTTTGCTGTGTGACGGGTTTTTTTTGGGCGTTGCCCGCATACCCCCAGGAAAAGAAAAAGGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         <  1:461555/150‑1 (MQ=255)
                                           tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                       >  2:884076/1‑151 (MQ=255)
                                           tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                       >  1:987547/1‑151 (MQ=255)
                                                gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                  >  1:1056126/1‑151 (MQ=255)
                                                gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                  >  1:1401366/1‑151 (MQ=255)
                                                  aCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                                 >  2:1261243/1‑150 (MQ=255)
                                                     tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  1:218422/1‑151 (MQ=255)
                                                     tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  1:40936/1‑151 (MQ=255)
                                                      aaCTCAAACACCGTTTGCGGTGTGACGGGTTGTTTTTTGGCGATGCCCGCACACCCCCGGGAAAATTAATAGGTTTTCCCTAAAAAAATTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             <  2:35462/150‑1 (MQ=255)
                                                      aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  2:383415/1‑150 (MQ=255)
                                                       acaaaaccacggTTTGTTGTGTGAGGGGCTGTTTTTTGGCGATGCCCCCATCGCCCAGGGAAAAGAAAAAAGTTTTCCCCTAAAAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                            <  2:619440/139‑1 (MQ=255)
                                                        ctaaaaCAACGTTTGCTGTGTGGGGGGTTTTTTTTGGGCGATGCCCGCACACCCCCAGTAAAAAAAAAATGTTTTCCAGAAAAAATTTCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                            <  2:1199858/146‑1 (MQ=255)
                                                           aaaCAACGTTTGCTGTGGGCCGGGCTGTTTTTTGGCGATGTCCCCATACCCCCAGGAAAAGAAATAAGTTTCCCCTTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          <  1:331136/148‑1 (MQ=255)
                                                                  gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                 >  2:509188/1‑150 (MQ=255)
                                                                   tttctggggggaggggCTGTTTTTGGGGGATGCCCCAATACCCCAAGAAAAAGTAAAAAGTTTTCCCTTAAAAATTACCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                  <  1:347561/135‑1 (MQ=255)
                                                                   tttGTTGTGGGACGGGCTGTTTTTGGGCGTTGTCCCCATACCCCCAGAAAAATTAATAAGTTTTCCCTTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                 <  2:1429019/149‑1 (MQ=255)
                                                                    ttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                              >  2:1230179/1‑151 (MQ=255)
                                                                     tgttgtgtgAGGGGCTTTTTTTTGGCGTTGTCCGCATCCCCCAAGAAAAAGAAAAAAGTTTCCCCTTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                              <  1:1291855/147‑1 (MQ=255)
                                                                      tgttgggggaggggTGTTTTTTGGGGTTTGTCCGCAGAGCCCAGAGAAAAGTAATAGTTTTCCCCTTAAAAATTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                               <  2:1394155/139‑1 (MQ=255)
                                                                         gtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         >  2:641111/1‑151 (MQ=255)
                                                                           gtggggggttttttttgGGGGATGCCCACATCCCCCCAGAAAATGTATTTAGTTTTCCCGTAAAAATTTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             <  2:806637/133‑1 (MQ=255)
                                                                                 ggCTTCTTTTTGGCCGTTCTCCACATACCCCGAGAAAAAGTAAAGAGTTTTCCTCGAATAATTCACCCCTACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                               <  1:1448375/138‑1 (MQ=255)
                                                                                        ttttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                           >  2:468017/1‑150 (MQ=255)
                                                                                         tttGTCGTTGTCCCCACACCCCAGGGAAATGTAATAAGTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgt                                                          <  2:50895/150‑1 (MQ=255)
                                                                                         tttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         >  2:436311/1‑151 (MQ=255)
                                                                                         tttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         >  1:797500/1‑151 (MQ=255)
                                                                                            ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                        >  1:717054/1‑149 (MQ=255)
                                                                                            ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                        >  1:536197/1‑149 (MQ=255)
                                                                                                 gagccccGCATACCCCCAGGAAAAGTAATAGGTTTTACCTTAAAAAATTCCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                        <  1:959159/141‑1 (MQ=255)
                                                                                                 cgatccccccacaccccAGGGAAAAGAAAAAAGTTTTCCCCTGAAAATTCCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                     <  2:746843/137‑1 (MQ=255)
                                                                                                      ccGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                              >  1:891984/1‑151 (MQ=255)
                                                                                                          atccccccAGAAAAAAAAATAAGTTTTCCCCTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            <  1:186016/145‑1 (MQ=255)
                                                                                                             gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       >  1:406760/1‑151 (MQ=255)
                                                                                                             gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       >  1:1010776/1‑151 (MQ=255)
                                                                                                              ccccAGTAAAAGTAAAATGTTTTCCCTAAAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                          <  2:343184/147‑1 (MQ=255)
                                                                                                              ccccAGAAAAAGAAAAAAGTTTTCCCTTAATAATTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                          <  1:918585/147‑1 (MQ=255)
                                                                                                                ccgcgggaaaataaataaGTTTGCCCTTGATAAATTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       <  2:540640/137‑1 (MQ=255)
                                                                                                                  agttAAATTAATAGGTTTCCCCCTAATAAATCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                         <  1:672210/141‑1 (MQ=255)
                                                                                                                  aGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    >  1:1371817/1‑149 (MQ=255)
                                                                                                                         aaataaaaaGGTTTTCCCTTAAAAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                               <  1:593871/147‑1 (MQ=255)
                                                                                                                                 aGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                     >  2:1393329/1‑151 (MQ=255)
                                                                                                                                     ttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAACGTCGttt                 >  2:601984/1‑151 (MQ=255)
                                                                                                                                       ccccTAAAAAATCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                <  1:160370/149‑1 (MQ=255)
                                                                                                                                             gATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa           >  2:1040788/1‑150 (MQ=255)
                                                                                                                                              aTAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                 <  2:17827/143‑1 (MQ=255)
                                                                                                                                              aTAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa           <  1:444466/149‑1 (MQ=255)
                                                                                                                                                aaaTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa           <  1:645607/147‑1 (MQ=255)
                                                                                                                                                  aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGTAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt     >  1:1433421/1‑151 (MQ=255)
                                                                                                                                                  aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      >  1:776071/1‑150 (MQ=255)
                                                                                                                                                  aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      >  2:1339332/1‑150 (MQ=255)
                                                                                                                                                     tACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt   >  2:425859/1‑151 (MQ=255)
                                                                                                                                                       accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:216855/1‑151 (MQ=255)
                                                                                                                                                       accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGGTTTGAACATAGGGTCt   >  2:475895/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378658‑4378952

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