Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC49.0% 72.8 / 14.5 51V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (25/1);  new base C (0/25);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.09e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378657‑4378948
                                                                                                                                                      |                                                                                                                                                  
tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgcc                                                                                                                                               >  2:602256/1‑151 (MQ=255)
tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGATTACCCTGATAAATTaccgcc                                                                                                                                               >  2:759644/1‑151 (MQ=255)
 cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgcc                                                                                                                                               >  1:499622/1‑150 (MQ=255)
    gCAATAGCGTCGCCACCAACCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccacc                                                                                                                                            <  2:972362/150‑1 (MQ=255)
     cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccgc                                                                                                                                          >  2:683396/1‑151 (MQ=255)
            gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATAccc                                                                                                                                   >  1:814009/1‑151 (MQ=255)
            gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATAccc                                                                                                                                   >  2:734220/1‑151 (MQ=255)
            gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATAccc                                                                                                                                   >  1:391647/1‑151 (MQ=255)
             tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCg                                                                                                                                  >  2:141746/1‑151 (MQ=255)
                   ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCt                                                                                                                            >  2:204777/1‑151 (MQ=255)
                          cgcaataaattGTTTTGCTGCCCGACCCAAACAAAAACAATTTTTGTTGTGTGAGGGTTTGTTTTTGGGGGTTGCCCGCATACCCCCGGTAAAAGAAAAAAGTTTCCCCTTAAAAAATCCCCCCCCCCGCCATACCCGGGATCTTttgtt                                                                                                                      <  1:783327/140‑1 (MQ=255)
                          ccccATTAATTTTTTCGCTGCCCAACCCAAACAAAAAAAGCGTTGGTTGTGTGCGGGGGTTTTTTTTGGCGATCCCCGCACGCCCCCAGAAAAAGAAATAGGTTTTCCCGTAATAAATACCCGCCCCCGCCATACCCGGGATCTTttgttg                                                                                                                     <  2:692850/149‑1 (MQ=255)
                          cGCCAATAATTTTTTTCCTGCCCGAACATACCTCAAAAAACTTTTGCTGTGTGACGGGTTTTTTTTGGGCGATCTCCGCATCGCCCAAGATAAAGTAAAAGGTTTTCCCTTAATAAATTCCCGCCACCGCCATACCCGGGATCtttt                                                                                                                         <  1:387525/147‑1 (MQ=255)
                               gTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 >  1:707514/1‑150 (MQ=255)
                                    ttttttGCGGACCAAACCAACCAAAAACAACGTTTGCTGTGGGGGGGGTTGTTTTTTGGCGAGGCCCCCCCCCCCCCAGAAAAAGAAAAATTTTTTCCCATAATAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                             <  2:32283/147‑1 (MQ=255)
                                    tGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                           >  1:810593/1‑151 (MQ=255)
                                              accaaccttaacaaaaaCAACTTTTGGTGTGTGGGGGGTTGTTTTTGGGCGATGCCCGCATACCCCAAGAAAAAGTAATAGTTTTTCCCCTAAAAAATCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                       <  2:72501/131‑1 (MQ=255)
                                               ccgacccaaccaaaaaCAGCTTTTGTGGTGTGAGGGGTTTTTTTTGGGCGATGCCCGCATACCCCAAGAAAAAGTAATAATTTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                    <  2:896170/134‑1 (MQ=255)
                                               ccaaaCCAAACACAACCATCGTTGGTTGTGGGAGGGGCTGCTTTTTGGCGATGCCCGCATAGCCCCAGGAAAATTAATAGTTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                 <  1:201432/147‑1 (MQ=255)
                                               ccGAACCTAACTAAAACACCGTTTGCTGGGTGACGGGGTTTTTTTGGGCGATTTCCGCACAGCCCGAGGAAAAGTAATAAGGTTTACCCGAATAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                      <  1:119968/145‑1 (MQ=255)
                                                  acccaaacaaaaaCAACTTTTGGTGGGTGAGGGGCTTTTTTGTGGGGATGCCCGCATCCCCCGAGAAAAAGTAATAGGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             <  2:125021/141‑1 (MQ=255)
                                                    ccTAACACAAACAACTTTTGCTGTGTGGGGGGTTGTTTTTGGGGGATGCCCGCATACCCCAGGAAAAAGTAATAAGTTTTCCCCGAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                           <  1:916946/151‑1 (MQ=255)
                                                          tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                      >  2:99911/1‑150 (MQ=255)
                                                           aaaaCAATGTTTCTTGTGTGGGGGGCTGTTTTTTGGTGTTGCCCTTTAATCCCCAGGTAAAGTAAAAGGTTTCCCCTTAAAAATTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    <  1:765855/150‑1 (MQ=255)
                                                            aaaCAACTTTTGTTGTGTGAGGGTCTTTTTTTTGGCGATCTCCCCATACCCCCAGGTAAAGTAATAGTTTTTACCTTAATAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    <  2:65104/150‑1 (MQ=255)
                                                               caccttttGTTGTGTGGGGGGCTTTTTTTGGGGGATGTCCCCATCGCCCCAGGAAAAGTAAAATGTTTTCCCCTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                 <  2:399453/145‑1 (MQ=255)
                                                                     ttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                          >  2:872201/1‑151 (MQ=255)
                                                                       gtggtgtgACGGGTTTTTTTTGGGCGTTTCCCGCTTCCCCCCAGAAAAATTAATATTTTTTCCTCTAAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                         <  2:577782/147‑1 (MQ=255)
                                                                          gtgtgCCGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAATAAAAAGGTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       <  1:6009/149‑1 (MQ=255)
                                                                           tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCAGGGATCTTTTGTTGCTCTATAAGCGGGGCAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     >  2:774400/1‑150 (MQ=255)
                                                                                cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                               >  2:595069/1‑151 (MQ=255)
                                                                                        tttttGGCGATGTCCCCACACCCCGAGTAAAAGTAATATGTTTTCCCTAAATAAATCCCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        <  1:583530/150‑1 (MQ=255)
                                                                                        ttttgGGGGATGTCCGCATACCCCCAGTAAAAGAAAAAAGTTTTCCCCTAAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                          <  1:123504/148‑1 (MQ=255)
                                                                                           tgggCGAGGCCCGAAAACCCCCAGAAAAAGTAATAGGTTTTCCCATAAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                     <  1:869698/148‑1 (MQ=255)
                                                                                               cGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                 >  2:90673/1‑150 (MQ=255)
                                                                                               cGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                 >  1:269511/1‑150 (MQ=255)
                                                                                                  tGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  1:514181/1‑148 (MQ=255)
                                                                                                    tCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           >  2:796283/1‑151 (MQ=255)
                                                                                                      cGCATGCCCCCGGGAAAAGTAATAAGTTTTCCCCAAAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                         <  1:914394/151‑1 (MQ=255)
                                                                                                          tAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                      >  2:225010/1‑150 (MQ=255)
                                                                                                            cccccGGTTAAAGAAAAATGTTTTCCCCTAATAATTTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   <  1:213305/150‑1 (MQ=255)
                                                                                                             ccccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   >  2:756685/1‑150 (MQ=255)
                                                                                                             ccccAGGAAAATAAATAATTTTTCCCCTAATAATTCCCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                  <  2:433286/151‑1 (MQ=255)
                                                                                                                  ggTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                >  1:722090/1‑148 (MQ=255)
                                                                                                                       aaataaaaaAGTTTTCCCTGGATAATTCACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                          <  2:828166/148‑1 (MQ=255)
                                                                                                                            ataaGTTTTACCTTGATAAATACCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       <  1:759729/148‑1 (MQ=255)
                                                                                                                              aaGGTTTACCCTGATAAATTACCGCCACCGCCATACCAGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTACATAAAAAGACGCAGAGcg                   >  2:213338/1‑150 (MQ=255)
                                                                                                                                   ttttCCCTAAAAAATACCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                    <  2:907818/145‑1 (MQ=255)
                                                                                                                                   tttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt              >  2:239765/1‑151 (MQ=255)
                                                                                                                                    ttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             >  2:133969/1‑151 (MQ=255)
                                                                                                                                     tttcccttaaTAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                <  2:813928/140‑1 (MQ=255)
                                                                                                                                     tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt            >  2:539597/1‑151 (MQ=255)
                                                                                                                                           tGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc        >  1:158114/1‑150 (MQ=255)
                                                                                                                                              aTAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      >  1:831768/1‑150 (MQ=255)
                                                                                                                                                aaaTTCCCCCCACCGCCATTCCCGTGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      <  1:183323/148‑1 (MQ=255)
                                                                                                                                                 aattcaccgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATa     <  2:499571/144‑1 (MQ=255)
                                                                                                                                                   ttACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  2:112720/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378657‑4378948

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