Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9010TA52.6% 7.6 / 25.0 19Y44F (TAT→TTT) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base A (1/9);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 5.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-03

CAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGC  >  NC_000913/4336751‑4337051
                                                                                                                                                       |                                                                                                                                                       
cagTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                         >  2:214729/1‑151 (MQ=255)
 agTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATa                                                                                                                                                       <  1:357750/151‑1 (MQ=255)
  gTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATa                                                                                                                                                       <  1:11946/150‑1 (MQ=255)
                  tGAAAACCGAGAAACGGGCAAAAGCAGCGGCGGGCGTAAGCGTAAGAACCGCCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCa                                                                                                                                       <  1:152028/150‑1 (MQ=255)
                                       aggcaggGGGGGGGGTAAGCGTAAGAACCACCACGACTTGGGTCGGGGAACGACATTTGGGCGTATACCATCGAGGGCCCCAGTGCACCGATAATCGTCTCCACCCATCCAAAAAGGGCAGTCCCGCCAGTAGAGGCCAGGTTTGCAGGTa                                                                                                                 <  1:303611/149‑1 (MQ=255)
                                              ggcgggcgTAAGCGTAAGAACCCCCGGAACTTGGGTCGAGGAACGCCATTGTGGCGTATACCATCGAGAACCCCAGTGCAGCGAAAATCCGCACCACCCATCAAAAAAGGGCAGCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaa                                                                                                          <  1:25187/151‑1 (MQ=255)
                                                   gcgaaagcgtaaGACCCCCCAGGACTTGGGTCGAGGACCGACTTTTTGGCGTATCCCATCGAGACCCCCAGTGCCCCGAAAATCGCCCCCAACCATCCAAAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAacac                                                                                                     <  2:176731/151‑1 (MQ=255)
                                                         gcttaaaAACCCCCAGGTTTGGGGCCGGGAAACGCCATTTTGGCGTAACCCAACGAGAGCCCCACGGCCCGGACAATGGTCACCAACCACCCAAAAATGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa                                                                                                 <  2:137347/146‑1 (MQ=255)
                                                                 accaccAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                                                                                        >  2:154634/1‑150 (MQ=255)
                                                                            tggggACGAGAACCGACATTTTGGGGTATCCCATGGAGAGCCCCAGTGCACCGTTAATCGCCCCCACCCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGaca                                                                             <  2:202324/148‑1 (MQ=255)
                                                                                   gAGGACCGACATTTGGGCGTATCCCATCGAGACCCCCAGTGCCCCAATAATCGTCCCCACCCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCa                                                                       <  2:82781/149‑1 (MQ=255)
                                                                                   gAGGAACGACTGTTGGGCATATCCAGCGGAGGGCCCCAGTCCCCCAATAAGCGTCAGCAACCATCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACGCCATCAg                                                                      <  1:81844/150‑1 (MQ=255)
                                                                                          gACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                              >  2:132885/1‑151 (MQ=255)
                                                                                                          ccATCAAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAGAAAAACACCAGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGAATTAAGCCCACttt                                              >  1:1385/1‑151 (MQ=255)
                                                                                                               cgagccccccAGTCCCCAGATAATCGTCCCCACCCATCCAAAAATGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACttt                                              <  1:377751/141‑1 (MQ=255)
                                                                                                                                cGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCGGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATcc                              >  1:168516/1‑146 (MQ=255)
                                                                                                                                  aTAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGAATCAGAATCCGAAGAc                        >  1:305996/1‑150 (MQ=255)
                                                                                                                                    aaccGTCACCAACCATCAATAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACa                       <  2:67684/146‑1 (MQ=255)
                                                                                                                                    aaTCGCCACCACCCACCCAAAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGa                         <  2:259562/147‑1 (MQ=255)
                                                                                                                                         tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                 >  2:33303/1‑150 (MQ=255)
                                                                                                                                          tccccgccacccataAAAAGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtg                  <  1:31712/143‑2 (MQ=255)
                                                                                                                                                ccatccaaAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc            <  2:310278/148‑1 (MQ=255)
                                                                                                                                                       tAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGc  >  2:172407/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGC  >  NC_000913/4336751‑4337051

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