Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,595,0440TA48.6% 26.5 / 19.2 35Y146F (TAC→TTC) livMbranched‑chain amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (14/4);  new base A (0/17);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCACC  >  NC_000913/3594900‑3595181
                                                                                                                                                   |                                                                                                                                         
aGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAAc                                                                                                                                     >  1:181635/1‑151 (MQ=255)
          cgcgCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTAccc                                                                                                                           >  2:39172/1‑151 (MQ=255)
                     agcaCCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGCGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCccagcaccagc                                                                                                                 >  2:321308/1‑150 (MQ=255)
                        ccccgAAAACCGGGGGGGGAGCCCGCCCCCCGCGCCATTAAACCCACAAACCGGCGGGAGGTCCAAGAGGCCCAGGCCGTAATAGTGATTGAGCAGCACAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACcagcag                                                                                                               <  1:141182/146‑1 (MQ=255)
                        aCCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGAc                                                                                                             >  2:463678/1‑151 (MQ=255)
                            gAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGAcc                                                                                                            >  2:318385/1‑148 (MQ=255)
                             aaaaCCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGAGATTGAGCAGCGCAAAAGTGTAAGCGCCGCTGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGa                                                                                                         >  1:211034/1‑150 (MQ=255)
                                   gAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGAcc                                                                                                   >  2:173775/1‑150 (MQ=255)
                                      aagaaacccggccgccgtTGCCATAAACCCGAAATTGGCAAGGAAGGTCAAAAACCCCAATCCGAAAGAGTGATTGAGCAGCGAAAAGGTGTAGGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACcaac                                                                                                <  2:401515/140‑1 (MQ=255)
                                        gaacgCCCGCCGCCGCTCCCATAAACCAACAAACGGGCGGGCGGGTCAAAAAGCCCAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACcaacaa                                                                                              <  2:293719/146‑1 (MQ=255)
                                               ccccgccgcTGCCATTAACCCAGAAACCGGCGGGAGGGTCAAAAGGCCCAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTc                                                                                      <  2:43260/149‑1 (MQ=255)
                                                      gcTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCg                                                                                >  1:198333/1‑150 (MQ=255)
                                                             ttAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCg                                                                          >  2:132910/1‑149 (MQ=255)
                                                                       aaccGGCAGGCAGGTCCAGACCCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGAt                                                                     <  1:463559/141‑1 (MQ=255)
                                                                            ggcgggcgGGCCCAGAAGCCAAAGCCGTAATGGTGATTGAGCGGCGAAAAAGTGTAAGCGCCGATGGCGAAAAACCCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGt                                                                <  1:490545/145‑1 (MQ=255)
                                                                                ggccaggTCCAGACCCCCAACCCGAAATAGGGATGGCGCAGGGCAAAATTGTAACCCCGAAGGGAGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATgg                                                       <  1:173715/148‑1 (MQ=255)
                                                                                aggcaggTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATgg                                                       >  2:127182/1‑151 (MQ=255)
                                                                                 gccgggTCAAGACGCCCAGGCGGTAATAGTGATGGAGCAGCGCAAAAGTGTAGGCGCGGATGGGATAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATg                                                        <  2:77653/145‑1 (MQ=255)
                                                                                   caggTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGTCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTc                                                     <  1:360543/150‑1 (MQ=255)
                                                                                     gggcaagaaaGCCAAAGCTGAAAAAGTGTTTGAGCAGCGCAAAAGTGAAAGCGCGGATGGAGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAg                                                   <  1:282618/143‑1 (MQ=255)
                                                                                     ggcccAGACCCCCAAGCCGTAATGTTGATAGAGAAGCGCAAAGGGGTAAGCGCCGTGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTc                                                     <  1:313654/145‑1 (MQ=255)
                                                                                                cccAAGCCGTAATAGTGATTGAGCAGCGAAAAAGTGAAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTgg                                              <  2:54918/145‑1 (MQ=255)
                                                                                                 caaaGCCGTAATAGTGATTGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATc                                       <  2:228846/149‑1 (MQ=255)
                                                                                                   aaGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATcc                                      >  2:351003/1‑150 (MQ=255)
                                                                                                      ccGTAAGAGAGATTGGGCGGCGCAAAAGTGGAAGCGCCGGGGGCAAAAAAACCTCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCg                                  <  2:179596/151‑1 (MQ=255)
                                                                                                              gTGATGGAGCAGCGAAAAAGGGAAAGCGCCGAGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGt                          <  1:268915/151‑1 (MQ=255)
                                                                                                                  ttGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGaa                       >  1:181072/1‑150 (MQ=255)
                                                                                                                  ttGAGCAGCGAAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTg                         <  2:181695/148‑1 (MQ=255)
                                                                                                                     accatcgaaaaaGTAAAACCGCCGATGGCATAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGaa                       <  1:48336/139‑1 (MQ=255)
                                                                                                                      gcagcGAAAAAGTGTAAGCGCCGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAAcc                    <  2:181132/149‑1 (MQ=255)
                                                                                                                          cgcAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAAc              >  1:399762/1‑151 (MQ=255)
                                                                                                                             aaaaGTGTAGGCGCGGGGGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACgg            <  1:319298/150‑1 (MQ=255)
                                                                                                                               aaGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGgcc          >  2:185093/1‑150 (MQ=255)
                                                                                                                                  tggaaGCGCCGATGGCGAAAAACCCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGgccac        <  1:533604/146‑1 (MQ=255)
                                                                                                                                          ccGATGGCGAAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCAcc  <  2:474054/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
AGATAGTCACCGCGCAAACGCAGCACCGGAAAACCGAGCAGGAAGCCCGCCGCCGCTGCCATTAATCCAGCAATCGGCAGGCAGGTCCAGAAGCCCAAGCCGTAATAGTGATTGAGCAGCGCAAAAGTGTAAGCGCCGATGGCGTAAAAACCGCCGTACCCCAGCACCAGCAGACCAGAAAGACCAACAACCACGTTCAGCCCGAGACCGAGGATAATGTAGATCATGGTCAGGGTGGCAATATCCACCGTCCCGCGTGAAACCATAAACGGCCACGCCACC  >  NC_000913/3594900‑3595181

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: