Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,9040AC56.0% 18.8 / 29.0 25V267G (GTA→GGA) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base A (8/3);  new base C (0/14);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 1.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-06

AATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAG  >  NC_000913/1613761‑1614048
                                                                                                                                               |                                                                                                                                                
aaTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTc                                                                                                                                            >  1:70770/1‑150 (MQ=255)
                                     gTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCAc                                                                                                       >  1:486768/1‑150 (MQ=255)
                                                  cccATTTTCAAGGCGGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGAt                                                                                         <  2:63748/151‑1 (MQ=255)
                                                              cctgctgCAGTTGCCAAAAACCGTCCGTAAATTGCCGAATCATGTCCCTCTCCGAAAAAAAAGCGTTTCGCCGCGCCACATACCTGTACGGTTTTCTGATAAGGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTa                                                                              <  1:271593/149‑1 (MQ=255)
                                                                      agctaaccaaaaaCGATAGGTAACGCCCGAAGAAACTCCTTCGTGGATAATAAACCGTTTTGCGGCTGAAAACCCCTTTCCTGTATTCGGATAACCACCGGCTCCCGCCTCTTCCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAAt                                                                      <  1:625918/141‑1 (MQ=255)
                                                                        cTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAAt                                                                      <  2:66121/148‑1 (MQ=255)
                                                                          ccaacaAACCGTTGGGTAAAGCCCAACGCATTTCCCTTTTCTAAAAAAAAGCTTTTGGCCGCTGCACATCCCTGTCGGTTATTTGGATAACGTCGGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACg                                                                 <  1:451883/148‑1 (MQ=255)
                                                                            caaaaaaCGTTCGGTAAATCCCAAACAATTTCCCTCTTCGATAAAAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTAATAACGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGa                                                               <  1:141466/148‑1 (MQ=255)
                                                                                   cGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAAc                                                        <  2:335059/151‑1 (MQ=255)
                                                                                   cGTTCGGTAAATCCCAAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAAc                                                        <  2:607582/151‑1 (MQ=255)
                                                                                        ggTAAATGCCGAAGAAATTCCCTCTTGGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGcc                                                    <  2:393905/150‑1 (MQ=255)
                                                                                         gTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCa                                                  >  2:426230/1‑151 (MQ=255)
                                                                                              tGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTcc                                             <  2:522366/151‑1 (MQ=255)
                                                                                                   aaccattccccTTTTAGATAAAAAAGCTTTTTGCCGCTGCAAACACCTGTCCGTTATTCTGATAACGGCCGGCTTCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTAc                                        <  2:676288/143‑1 (MQ=255)
                                                                                                   aaGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTa                                         <  2:107505/150‑1 (MQ=255)
                                                                                                   aaGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTa                                         <  2:143769/150‑1 (MQ=255)
                                                                                                   aaGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTa                                         <  2:580341/150‑1 (MQ=255)
                                                                                                   aaGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGAAAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTa                                         <  2:633447/150‑1 (MQ=255)
                                                                                                        aTTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtt                                   >  2:42297/1‑151 (MQ=255)
                                                                                                             ctctTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtttttt                               >  1:37049/1‑150 (MQ=255)
                                                                                                                    ataataAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGc                                                                           <  2:176355/99‑1 (MQ=255)
                                                                                                                    ataataAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGc                                                                           >  1:176354/1‑99 (MQ=255)
                                                                                                                             cGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGCCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc               <  2:620391/150‑1 (MQ=255)
                                                                                                                                         gCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAg  >  1:228479/1‑151 (MQ=255)
                                                                                                                                         gCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAg  >  2:82601/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAG  >  NC_000913/1613761‑1614048

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