Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,5970GA44.8% 30.0 / 11.3 29F184F (TTC→TTT) yqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (14/2);  new base A (0/13);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.75e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTT  >  NC_000913/3042449‑3042737
                                                                                                                                                     |                                                                                                                                            
gACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGaa                                                                                                                                            <  1:595589/151‑1 (MQ=255)
 aCGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGaa                                                                                                                                            >  1:164242/1‑150 (MQ=255)
 aCGCTGGCGCGTCTTATCAGGCATACAAAGGCATACCCATTACGCCTGCCCCATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGCTGTATGGAATGCGTTTGC‑AGACGTAGaa                                                                                                                                            >  2:539110/1‑150 (MQ=255)
      ggCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATc                                                                                                                                        <  2:524505/149‑1 (MQ=255)
                 tCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAAt                                                                                                                           >  1:418243/1‑151 (MQ=255)
                               gCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCt                                                                                                             >  1:119832/1‑151 (MQ=255)
                                 aTACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTAc                                                                                                           >  1:113179/1‑151 (MQ=255)
                                               tGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTa                                                                                              >  2:222575/1‑150 (MQ=255)
                                                 cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACg                                                                                           >  1:304999/1‑151 (MQ=255)
                                                                   cGCCAAAAAGTGGCACACACTCCCGCCGAGCACGAAGCCGTGCCAGATGGCAGGTTTGTATGAAATGCTTTTGC‑AGACGTAAAAAATCCCCCCGAGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTa                                                                         <  1:159060/151‑1 (MQ=255)
                                                                       aaaaaGTGCCCCCCACTACCGCCGAGCACAAGGCCGTGCCAGAGGGCTGGTTTTTAGGAAATGCTTTTGC‑AGACGTAAAAAACCCCCCCGGGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                      <  1:170200/148‑1 (MQ=255)
                                                                       aaaaaGGGGCACACCCTCCCGCCGAGCAGGAAGCGGTCCAAGATGGCATGGTGGTATGAAATGGTTTTGA‑AGACGTAAAAAACCCCCCCGGGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTACCGCTGCCCGCCGCGAGCTTAACTg                                                                     <  1:57976/149‑1 (MQ=255)
                                                                                    cacTCCCGCCGAGCACGAACCCGTGCCGGATGGCATGGTTGTATGAAAGGCTTTTGC‑AGACGAAAAAAACCCCCCCGAGGGAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATaa                                                          <  2:554555/149‑1 (MQ=255)
                                                                                       tACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                     >  2:392564/1‑151 (MQ=255)
                                                                                       tACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                     >  2:234565/1‑151 (MQ=255)
                                                                                         ccgcggAGCACGAACCCGGCCCAGTTGGCTGGGTTGTTTGGATGGCTTTGGC‑GGCCGTAAAAATCCCCCCCGGGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                                                      <  2:611473/148‑1 (MQ=255)
                                                                                           gccgAGCACGAAGCGGGCCCAGAGGGCATGGTTGTATGAAAGGCTTTTGA‑AGCGGAAAAAAACCCCCCGGAGGAAATAAACCCCGCCGCCTCCCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccac                                                  <  2:497281/150‑1 (MQ=255)
                                                                                           gccgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccacc                                                 >  2:395544/1‑151 (MQ=255)
                                                                                                 cACGAAGCGGTGCAAGATGGCAGGTTTGTAGGAAAGGCGTTTGC‑AGACGTAAAAAATCCCCCCGGGCAAAAAAACCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGc                                              <  1:621249/148‑1 (MQ=255)
                                                                                                    gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                        >  2:212124/1‑151 (MQ=255)
                                                                                                     aggCCGTGCAAGATGGCATGGTTGTATGAAAGGCTTTTGC‑AGACGTAAAAAACCCCCCCGAGCAAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCg                                             <  2:609643/143‑1 (MQ=255)
                                                                                                                    gCATGGTTGTATGGAATGCGTTTGC‑AGACGTAAAAAATCCCCCCGAGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                        <  1:519025/151‑1 (MQ=255)
                                                                                                                        ggTTGTAGGAAATGC‑TTTGCAAGACTAAAAAAACCCCCCCGGGCGAAAAAACCCCGCCGCCTACCGCCAGTAAGGAAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                                                      <  2:458097/117‑1 (MQ=255)
                                                                                                                            gTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                >  1:48647/1‑151 (MQ=255)
                                                                                                                             tATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                >  1:295734/1‑150 (MQ=255)
                                                                                                                             tATGGAAGGCTTTGGC‑AGAGGAAAAAAATCCCCCCGAGCGAATAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAc                 <  1:402199/149‑1 (MQ=255)
                                                                                                                                ggAATGCGTTTGC‑AGCCGTAAAAAACCCCCCGGAGCAAAAAACCCACGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCa               <  1:395555/148‑1 (MQ=255)
                                                                                                                                     gcttttGA‑AAACGTAAAAAACCCCCCCGGGCAAAAAAACCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACGGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATa         <  2:474691/146‑1 (MQ=255)
                                                                                                                                           gC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtttt  >  1:555186/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTT  >  NC_000913/3042449‑3042737

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