Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,6040AC56.7% 20.5 / 46.8 30V182G (GTG→GGG) yqfAhemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (13/0);  new base C (0/17);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 8.35e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05

TCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATG  >  NC_000913/3042466‑3042753
                                                                                                                                           |                                                                                                                                                      
tCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAAt                                                                                                                                           >  1:418243/1‑151 (MQ=255)
              gCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCt                                                                                                                             >  1:119832/1‑151 (MQ=255)
                aTACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTAc                                                                                                                           >  1:113179/1‑151 (MQ=255)
                              tGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTa                                                                                                              >  2:222575/1‑150 (MQ=255)
                                cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACg                                                                                                           >  1:304999/1‑151 (MQ=255)
                                                  cGCCAAAAAGTGGCACACACTCCCGCCGAGCACGAAGCCGTGCCAGATGGCAGGTTTGTATGAAATGCTTTTGC‑AGACGTAAAAAATCCCCCCG‑AGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTa                                                                                         <  1:159060/151‑1 (MQ=255)
                                                      aaaaaGTGCCCCCCACTACCGCCGAGCACAAGGCCGTGCCAGAGGGCTGGTTTTTAGGAAATGCTTTTGC‑AGACGTAAAAAACCCCCCCG‑GGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACt                                                                                      <  1:170200/148‑1 (MQ=255)
                                                      aaaaaGGGGCACACCCTCCCGCCGAGCAGGAAGCGGTCCAAGATGGCATGGTGGTATGAAATGGTTTTGA‑AGACGTAAAAAACCCCCCCG‑GGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTACCGCTGCCCGCCGCGAGCTTAACTg                                                                                     <  1:57976/149‑1 (MQ=255)
                                                                   cacTCCCGCCGAGCACGAACCCGTGCCGGATGGCATGGTTGTATGAAAGGCTTTTGC‑AGACGAAAAAAACCCCCCCG‑AGGGAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATaa                                                                          <  2:554555/149‑1 (MQ=255)
                                                                      tACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                                     >  2:392564/1‑151 (MQ=255)
                                                                      tACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                                     >  2:234565/1‑151 (MQ=255)
                                                                        ccgcggAGCACGAACCCGGCCCAGTTGGCTGGGTTGTTTGGATGGCTTTGGC‑GGCCGTAAAAATCCCCCCCG‑GGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                                                                      <  2:611473/148‑1 (MQ=255)
                                                                          gccgAGCACGAAGCGGGCCCAGAGGGCATGGTTGTATGAAAGGCTTTTGA‑AGCGGAAAAAAACCCCCCGG‑AGGAAATAAACCCCGCCGCCTCCCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccac                                                                  <  2:497281/150‑1 (MQ=255)
                                                                          gccgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccacc                                                                 >  2:395544/1‑151 (MQ=255)
                                                                                cACGAAGCGGTGCAAGATGGCAGGTTTGTAGGAAAGGCGTTTGC‑AGACGTAAAAAATCCCCCCG‑GGCAAAAAAACCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGc                                                              <  1:621249/148‑1 (MQ=255)
                                                                                   gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                        >  2:212124/1‑151 (MQ=255)
                                                                                    aggCCGTGCAAGATGGCATGGTTGTATGAAAGGCTTTTGC‑AGACGTAAAAAACCCCCCCG‑AGCAAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCg                                                             <  2:609643/143‑1 (MQ=255)
                                                                                                   gCATGGTTGTATGGAATGCGTTTGC‑AGACGTAAAAAATCCCCCCG‑AGCGAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGAt                                        <  1:519025/151‑1 (MQ=255)
                                                                                                       ggTTGTAGGAAATGC‑TTTGCAAGACTAAAAAAACCCCCCCG‑GGCGAAAAAACCCCGCCGCCTACCGCCAGTAAGGAAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                                                                      <  2:458097/117‑1 (MQ=255)
                                                                                                           gTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                                >  1:48647/1‑151 (MQ=255)
                                                                                                            tATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAcc                                >  1:295734/1‑150 (MQ=255)
                                                                                                            tATGGAAGGCTTTGGC‑AGAGGAAAAAAATCCCCCCG‑AGCGAATAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAc                                 <  1:402199/149‑1 (MQ=255)
                                                                                                               ggAATGCGTTTGC‑AGCCGTAAAAAACCCCCCGG‑AGCAAAAAACCCACGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCa                               <  1:395555/148‑1 (MQ=255)
                                                                                                                    gcttttGA‑AAACGTAAAAAACCCCCCCG‑GGCAAAAAAACCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACGGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATa                         <  2:474691/146‑1 (MQ=255)
                                                                                                                          gC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAtttt                  >  1:555186/1‑150 (MQ=255)
                                                                                                                                   aaaaaaCCCCCCCG‑AGCAAAAAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAttt                   <  2:14832/139‑1 (MQ=255)
                                                                                                                                    aaaaaTTCCCCCGAAGCAAAAAACCCAGCCCCCCAACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTg         <  1:660873/150‑1 (MQ=255)
                                                                                                                                     aaaaTCCCCCCG‑AGCAAATAACCCCCCCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATAgg           <  1:517564/147‑1 (MQ=255)
                                                                                                                                     aaaaTCCCCCCG‑AGCAAAAAACCCCCCCCGCCTCCCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgc        <  1:106056/150‑1 (MQ=255)
                                                                                                                                          cACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATg  >  2:590370/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                      
TCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGC‑AGACGTAGAAAATCACCCCG‑AGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCGCGATG  >  NC_000913/3042466‑3042753

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