Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,026,9410AC44.2% 50.1 / 16.6 43E392D (GAA→GACguaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (18/6);  new base C (0/19);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 2.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTA  >  NC_000913/3026795‑3027084
                                                                                                                                                    |                                                                                                                                               
tccAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAccc                                                                                                                                              >  1:402995/3‑150 (MQ=255)
 aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCAc                                                                                                                                            >  2:304543/1‑151 (MQ=255)
  cAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGGTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGATGCTGATTTCGTGGTGATGGAACCCAcc                                                                                                                                           >  1:113121/1‑151 (MQ=255)
                     ggCTATCTCCTCTCGGCATAGGAAGCTTTTTACCTGGCCCCGCTCGGCGGGGCAAAATCTCTGGGCCTTGACGTTTTGATTGGCAACTTTTTACCGGGAAAAGAGGCTAATTTCGGGGTGATGGACCCCACCGCCACTCCGCTAcagca                                                                                                                          <  2:387545/149‑1 (MQ=255)
                           cGCCTCTCGGCATATGAAGCTTTTTACCTGCCCCCGCTCGGGGGAGCAAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGAGGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                  <  1:139871/151‑1 (MQ=255)
                             cctcgcgcaatagGAGGCTTTTTGCCGGGCCCCGCTCGGCGGGGCAAACTCTCGGGGCTTTGCGGATTGAGTTGCCACTTTTTTCCCTGGCAAAGAGGTGGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgc                                                                                                                     <  2:228621/137‑1 (MQ=255)
                                 tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                             >  1:534105/1‑150 (MQ=255)
                                     catgtaaaGTGTTTTCCCGGGCCACGCCCGGGGGAGCGAATTCTCTGGGCTTTGACTTTTTGATTGGCACCTTTTTCCCGGGCAAAGAGGCTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                                                                                        <  2:505489/145‑1 (MQ=255)
                                           acgtgTTTTCCCTTGCCCCGCTCGTCGGGGAGAAATCTTTGGGGTTAAACTTTTTGATTGGCACCTTTTCCCCTGGCAAAGAGGCTGATTGCGTGGTGATAGACCCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                    <  2:195918/145‑1 (MQ=255)
                                            aagcttttACCCGGCCCCCGCCCGCGGGAGCGAATCCTCTGGGCCTTGCCGATTTATTTGGCATCTTTTTCCCGGGAAATGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                    <  1:432507/145‑1 (MQ=255)
                                                ttttTCCGGCCCGCGCTCGGCGGAGCAAAGCCGCTGGTCCTTGGCGTTTTGAGTGGAACTTTTTGCCCTGCCAAAGAGGCTTTTTCCGTGGGGAGGATCCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCttt                                                                                               <  2:611449/149‑1 (MQ=255)
                                                 tttCCCTGGCCCCGCTGGGCGGAGCAAAATCTCTGGGCCTTGACGTTTTGATTGGCATTTTTTTCCCTGGAAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                             <  1:521316/150‑1 (MQ=255)
                                                     ccTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGAc                                                                                        >  1:154831/1‑151 (MQ=255)
                                                           ccACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                   >  1:53991/1‑151 (MQ=255)
                                                            cccGCTGGGCGGAGCGAAATCTCTGGGCCTTGATGATTTGATTGGCAACTTTTCACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTg                                                                                  <  2:189942/149‑1 (MQ=255)
                                                               gCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc                                                                               >  2:381020/1‑151 (MQ=255)
                                                                cTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                              >  1:350369/1‑151 (MQ=255)
                                                                    gcggAGCGAAATCTCTGGGCTTTGACTATTTGATGGGCATCTTTTTCCCTGGCAAAGGGGCTGATTTCGTGGTGATGAACCCCGCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCg                                                                              <  2:414483/147‑1 (MQ=255)
                                                                      ggAGCGAAATCTCGGGCCCTTGACGATTTGATGGGCACCTTTTTCCCGGGAAAAGGGGTTGTTTTGGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatg                                                                         <  1:113426/150‑1 (MQ=255)
                                                                           gAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc                                                                    >  2:228461/1‑150 (MQ=255)
                                                                            aattctcGGGCCCTTGACGATTTGATGGACAACTTTTCCCCTGCCAAAGGGGTTGGTTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga                                                                     <  2:179775/145‑1 (MQ=255)
                                                                             attctcGGGGCCTTGACGTTTTGATTGGCAATTTTTTCCCTGGAAAAGGGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt                                                                 <  1:82866/149‑1 (MQ=255)
                                                                                ctcgGGGCCTTGACGATTTGATTGGAAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTACGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGtt                                                                 <  2:443278/148‑1 (MQ=255)
                                                                                 tcggggCCTTGACGTTTTGTTTGGCAACTTTTCACCCGGAAAGGACCCTAGTTCGGTGGTGTAGGACCCCTCCGCCACTCCGCTACAGCAGCTGCGATATTACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGAGGACGtt                                                                 <  2:163929/144‑1 (MQ=255)
                                                                                   tGGGCCTTGACGATTTGATTGGCAACTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGAGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTggg                                                              <  1:164147/148‑1 (MQ=255)
                                                                                   tGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCg                                                            >  1:112727/1‑150 (MQ=255)
                                                                                    gggCCTTGACGTTTTGATTGGCATTTTTTTCCCTGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTg                                                                <  1:511371/145‑1 (MQ=255)
                                                                                    gggCCTTGACGATTTGATGGGCACTTTTTTCCCGGGAAAAGAGGTTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGc                                                             <  1:220355/148‑1 (MQ=255)
                                                                                     ggCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATg                                                         >  1:592131/1‑151 (MQ=255)
                                                                                        cTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGAcc                                                      >  1:77803/1‑151 (MQ=255)
                                                                                          tGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGt                                                    >  2:594446/1‑151 (MQ=255)
                                                                                                             tttttCCCGGGAAGAGGGGCTGTTTTCGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCAc                                      <  2:237325/146‑1 (MQ=255)
                                                                                                              ttttCCCTGGCAAAGAGGCTGATTTCGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTAc                                  <  1:69048/149‑1 (MQ=255)
                                                                                                              ttttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGt                                >  1:612195/1‑151 (MQ=255)
                                                                                                                ttACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTg                              >  1:363731/1‑151 (MQ=255)
                                                                                                                   tccTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGAt                            >  2:218197/2‑151 (MQ=255)
                                                                                                                             ggggCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTgg                  <  1:195796/149‑1 (MQ=255)
                                                                                                                                  tGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTAc             >  1:568012/1‑151 (MQ=255)
                                                                                                                                    aTTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg            >  1:421987/1‑150 (MQ=255)
                                                                                                                                     tttCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACg            <  1:113377/149‑1 (MQ=255)
                                                                                                                                          tggtgATGAAACCCTCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCa      <  2:519456/150‑1 (MQ=255)
                                                                                                                                           ggtgATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAAc    <  1:606713/151‑1 (MQ=255)
                                                                                                                                              gagggAACCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTGCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGTGCGATGACCGTTCGATCTAACGCACCTACGTTGATGGTCGTCTGGTGCACGAACGCaagtg  <  2:113094/147‑4 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTA  >  NC_000913/3026795‑3027084

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