Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC40.7% 172.3 / 18.7 91V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (46/8);  new base C (0/37);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 5.07e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952
                                                                                                                                                     |                                                                                                                                                       
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTac                                                                                                                                                        >  2:728897/1‑151 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTa                                                                                                                                                         <  2:447609/150‑1 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGCTTACCCTGATAAATTac                                                                                                                                                        >  2:737338/1‑151 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccg                                                                                                                                                     >  1:147654/1‑151 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgcc                                                                                                                                                   <  1:1426271/150‑1 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgcc                                                                                                                                                   <  2:947458/150‑1 (MQ=255)
           aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccgcc                                                                                                                                             >  2:897820/1‑150 (MQ=255)
           aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccgcc                                                                                                                                             >  1:421678/1‑150 (MQ=255)
             aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTCCCGCCACCGCCAt                                                                                                                                           <  2:103764/150‑1 (MQ=255)
              gCGTCGCCCCCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATAc                                                                                                                                         <  1:447353/151‑1 (MQ=255)
                    ncaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGNNCGNNTAGCCCCAGGTAAAGTAATAAGGTTNNCCCTGATAAATTACCGCCACCGCCATACCCggg                                                                                                                                    >  1:105850/2‑150 (MQ=255)
                    ccaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCggg                                                                                                                                    >  2:961101/1‑150 (MQ=255)
                      cccaacccccccAATAATTTTTTTGCTGACCGAACCTAACACAAACACCGTTTGCTGTGGGACGGGCTGCTTTTTGGCGATGTCCGCATACCCCCAGAAAAAGAAAAAAGGTTTCCCTTAATAATTTCCCGCCACCGCCATACCCGGGATc                                                                                                                                 <  1:1192169/141‑1 (MQ=255)
                      cccaAGCCCGCCAGTAAATGTTTTGCTGACCGACCCAAACTAAAACACCTTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCCCATACCCCCAGGTAAAGAAATAAGGTTTCCCTTAATAAATTCCCGCCACCGCCATACCCGGGATc                                                                                                                                 <  1:298538/150‑1 (MQ=255)
                       ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCCTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCt                                                                                                                                >  2:236826/1‑151 (MQ=255)
                           gCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTttg                                                                                                                            >  1:430470/1‑151 (MQ=255)
                             ccGCCAGTAAATTTTTTGCTGCCCGAACCTAACTCAAACACCGTTTGCTGTGTGACGGGTTGTTTTTGGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGGTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          <  2:1349667/150‑1 (MQ=255)
                               gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        >  2:1028686/1‑151 (MQ=255)
                               gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        >  1:342510/1‑151 (MQ=255)
                                  aGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                     >  1:913905/1‑151 (MQ=255)
                                    tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:594763/1‑151 (MQ=255)
                                     aatttttttGCGGACCGACCCAAACTCAAACACCTTTTGTTGTGTGAGGGGTTGTTTTTGGGCGATGCCCGCATAGCCCAAGATAAAGTAATAAGTTTTCCCTTGAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   <  1:595198/145‑1 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:603824/1‑150 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:1007000/1‑150 (MQ=255)
                                     aaaTGTTTTGCTGACCAAACCTAACTAAAACAACGTTTGCTGTGTGGCGGGCTTTTTTTGGGCGATCCCCCCATAGCCCCAGGTAAAGTAATAAGTTTTCCCCTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  <  1:331413/151‑1 (MQ=255)
                                     aaaTGTTTGGCTGACCAAACCAAACTAAAAAAGGTTTTGCTGTGGGGCGGGTTGTTTTTGGGGGATGTCCGAATACCCCGGGAAAAAGAAAAAGGGTTTCCCTTGATAAATGACCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  <  2:939298/151‑1 (MQ=255)
                                         tttttGCTGACCGACCCTAACTAAAACAACTTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCAGGAAAAGTAATAAGTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                               <  2:1140360/149‑1 (MQ=255)
                                              gCTGCCCGACCCTACCCCAAACAGCTTTTGCTGTGTGAGGGGTTGTTTTTGGGCAATGCCCGCATAGCCCCAGTAAAATTAATAATTTTTCCCGTAAAAATTACCGGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         <  1:1416459/151‑1 (MQ=255)
                                              gCTGACCGACCCTAACTCAAACAACTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGTTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                          <  2:587227/150‑1 (MQ=255)
                                                 gACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                      >  2:868091/1‑151 (MQ=255)
                                                     gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                  >  1:1319340/1‑151 (MQ=255)
                                                       aCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                                 >  2:658873/1‑150 (MQ=255)
                                                           aacaaaaaCACCTTTTGCTGGGTGACGGGCTGTTTTTTGGCGATGTCCGCATACCCCCAGTAAAATTAAAAAGTTTTCCCTTAATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                            <  2:490911/146‑1 (MQ=255)
                                                           aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  1:200723/1‑150 (MQ=255)
                                                            acaaaaaCAACTTTTGCTGTGTGAGGGGCTTTTTTTTGGCGATGTCCGCATACCCCAAGGTAAAGAAAAAGGGTTTCCCCTGAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           <  2:41971/147‑1 (MQ=255)
                                                            aCTAAAACAACTTTTGTTGTGTACCGGGTTTTTTTTTGGCGTTGCCCCCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTGAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           <  1:1353355/151‑1 (MQ=255)
                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          >  1:700516/1‑150 (MQ=255)
                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          >  1:882515/1‑150 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:843279/1‑151 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:237938/1‑151 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:331346/1‑151 (MQ=255)
                                                               aaaaCACCGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGGTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        <  2:1167699/150‑1 (MQ=255)
                                                                aaccaccttttGCTGTGTGAGGGGCTTTTTTTTGGCTATGTCCCCATCCCCCAAGGAAAATAAATAAGGTTTCCCCTGAAAAATCCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                       <  2:943690/143‑1 (MQ=255)
                                                                 aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                      >  1:20649/1‑151 (MQ=255)
                                                                  acagcGTTTGCTGGGTGACGGGCTGCTTTGGGGCGATGGCGGAATCGCCCCGGGAAAAATAATAAGTTTTCCCCTAATAAATGCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     <  2:1138151/151‑1 (MQ=255)
                                                                       gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                 >  1:727929/1‑150 (MQ=255)
                                                                          tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             >  2:492126/1‑151 (MQ=255)
                                                                          tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             >  2:181206/1‑151 (MQ=255)
                                                                           gCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                            >  2:963712/1‑151 (MQ=255)
                                                                                gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         >  1:351145/1‑149 (MQ=255)
                                                                                    ggggtttttttttGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   <  2:1129024/143‑1 (MQ=255)
                                                                                    ggggttttttttgGGCGATGTCCGCATAGCCCAAGGAAAAGTAAAAGGTTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   <  2:1006024/143‑1 (MQ=255)
                                                                                     gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   >  2:1006913/1‑150 (MQ=255)
                                                                                       gcttttttttGGCGTGGTCCGAATAGCCCCGGGAAAAGAAATAAGGTTTCCCCTAATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                                <  2:1400060/146‑1 (MQ=255)
                                                                                         ttttttttGGCGATGCCCGCATACCCCCAGAAAAATTAATAAGTTTTCCCCTGATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               <  1:773313/147‑1 (MQ=255)
                                                                                         tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               >  2:1345765/1‑150 (MQ=255)
                                                                                          gttttttGGCGATGCCCGCATACCCCGAGGTAAAGTAATATGTTTTCCCCTAATAATTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               <  1:222620/147‑1 (MQ=255)
                                                                                             tttgGGCGTTGCCCGCATAGCCCCGGGAAAAGTAATAAGTTTTCCCCTAATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                           <  2:52838/150‑1 (MQ=255)
                                                                                               ttGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGGTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         <  2:1376375/150‑1 (MQ=255)
                                                                                               tgggCGATGTCCGCATAGCCCCAGGTAAAGTAAAAAGGTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         <  2:811784/148‑1 (MQ=255)
                                                                                                 ggCGATGTCCGCATAGCCCAAGGTAAAGTAATAAGGTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAg                                                      <  2:493826/151‑1 (MQ=255)
                                                                                                  gCGATGCCCGCATAGCCCCAGGTAAAGTAATAAGGTTTCCCTTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                     <  2:546994/151‑1 (MQ=255)
                                                                                                  gCGATGCCCGCATACCCCCAGGTAAAGAAATAAGTTTTCCCCTAATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                     <  2:1215714/151‑1 (MQ=255)
                                                                                                    gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                    >  2:41310/1‑150 (MQ=255)
                                                                                                    aaggTCCGCATACCCCGGGAAAAAATAATAGGTTTTCCCATTAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                   <  2:788448/148‑1 (MQ=255)
                                                                                                     aTGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                    >  1:344417/1‑149 (MQ=255)
                                                                                                        tCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCACCAAATGTGGAGCAAGAgg                                               >  1:911255/1‑151 (MQ=255)
                                                                                                        cccacATACCCCCAGGAAAAGTAATAGGTTTTCCCTTAATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                               <  1:1055220/147‑1 (MQ=255)
                                                                                                         ccGCATAGCCCCAGGAAAAGTAATAGGTTTTCCCCTGATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                              <  2:534510/151‑1 (MQ=255)
                                                                                                           gCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            >  2:1451575/1‑151 (MQ=255)
                                                                                                             aTAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                           >  1:1397297/1‑150 (MQ=255)
                                                                                                               aGCCCGAGGAAAATAAAAAAGTTTTCCCCTAAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                        <  1:401913/151‑1 (MQ=255)
                                                                                                               aGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                        >  2:400238/1‑151 (MQ=255)
                                                                                                                gCCCCAGGAAAATAAATAGGTTTTCCCCTGATAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       <  2:1309857/151‑1 (MQ=255)
                                                                                                                 ccccAGGAAAAGAAAAAAGTTTTCCCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                      <  2:1039390/151‑1 (MQ=255)
                                                                                                                  cccAGGAAAAGAAATAGGGTTTCCCTTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                     <  2:175217/151‑1 (MQ=255)
                                                                                                                   caggttAAAGTAATAGGTTTTCCCTTGAAAATTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    <  2:407546/146‑1 (MQ=255)
                                                                                                                   caagatAAAGAAATAAGTTTTCCCCTAATAAATAACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    <  1:742548/146‑1 (MQ=255)
                                                                                                                      ggaaaaaTAATAAGTTTTACCCTGATAAATTACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                 <  1:803426/148‑1 (MQ=255)
                                                                                                                      ggTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                   >  1:1319197/1‑149 (MQ=255)
                                                                                                                      gaaaaaGTAAAAAGGTTTCCCTGGATATTTTCCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                 <  2:319886/148‑1 (MQ=255)
                                                                                                                       gTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:553353/1‑151 (MQ=255)
                                                                                                                       gTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:321778/1‑151 (MQ=255)
                                                                                                                            gTAATAGGTTTTCCCCTAATAAATTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                           <  1:475903/151‑1 (MQ=255)
                                                                                                                            gTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:361272/1‑149 (MQ=255)
                                                                                                                            gTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:190828/1‑149 (MQ=255)
                                                                                                                             taataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:1289951/1‑148 (MQ=255)
                                                                                                                              aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                           >  2:353206/1‑149 (MQ=255)
                                                                                                                               ataaGGTTTCCCCTGATAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                        <  2:438556/151‑1 (MQ=255)
                                                                                                                                  aGGTTTANCTTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGNNANNCACAATATCGTTGATNNNNNNNGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                     <  2:105846/151‑1 (MQ=255)
                                                                                                                                       tccccTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                 <  1:1418795/148‑1 (MQ=255)
                                                                                                                                        aCCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg               >  1:780918/1‑151 (MQ=255)
                                                                                                                                             gATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  2:119403/1‑151 (MQ=255)
                                                                                                                                             gATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  2:334236/1‑151 (MQ=255)
                                                                                                                                              aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGTACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATCAGGCGCAGAGCGTCGTTTTGAACAt          >  1:876075/1‑150 (MQ=255)
                                                                                                                                                aaaTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       <  1:737345/151‑1 (MQ=255)
                                                                                                                                                aaaTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       >  2:339245/1‑151 (MQ=255)
                                                                                                                                                 aaTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      <  2:353707/151‑1 (MQ=255)
                                                                                                                                                  aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt     >  1:70872/1‑151 (MQ=255)
                                                                                                                                                     accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:172281/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952

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