Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,7480CA51.7% 20.5 / 44.6 29E55* (GAG→TAG) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (12/2);  new base A (0/15);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.75e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-04

CTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATG  >  NC_000913/3177600‑3177882
                                                                                                                                                    |                                                                                                                                      
cTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCa                                                                                                                                       >  1:1129801/1‑150 (MQ=255)
      cccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAttt                                                                                                                                >  1:1166195/1‑151 (MQ=255)
                                ccccAGTAGCCAAGGGGTTTCGCGGGTGTCGTACGGGGCAATCGGAATCTGCTCAAGCAGCCCAACTTCATCACGCCCGGGGTCAAAGGGACGCAAGTCTGGGGGGTGCCCGTGCCAAAAAATTTCCCGCCGCACCTCATGACTCATTTGc                                                                                                      <  1:4018973/149‑1 (MQ=255)
                                  ccAGTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTAAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGccc                                                                                                    <  2:1002036/151‑1 (MQ=255)
                                  ccAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTg                                                                                                       <  2:2688549/148‑1 (MQ=255)
                                       aGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGTTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa                                                                                               <  2:3867547/151‑1 (MQ=255)
                                              gggTTTGGCTGGTTTGGTACGCGGCAATCCGAATCTGCTCAATCACCACAACTCCATCCCGCACTGGGTCAAAGGGTAGCAATACTGCGGCTTGCCCCCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAtt                                                                                            <  1:555798/147‑2 (MQ=255)
                                                 tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                     >  1:1820046/1‑151 (MQ=255)
                                                  ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAACCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCACAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                     <  1:394821/150‑1 (MQ=255)
                                                                aCGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                        >  2:1400150/1‑149 (MQ=255)
                                                                 cgcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                     <  2:2351685/151‑1 (MQ=255)
                                                                    ggCAACCCGAATCTGCACAATCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCTTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                   <  1:2240252/150‑1 (MQ=255)
                                                                          ccGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                             >  2:2539396/1‑150 (MQ=255)
                                                                                cggcccAATCAGCACAACTCCATCCCCCACGGGGTCAAAGGTTAACAAGACTGCGGGGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTACCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                        <  1:1732056/144‑1 (MQ=255)
                                                                                  gCTCAATCAGCACATCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAg                                                      <  1:2618841/149‑1 (MQ=255)
                                                                                     cAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                 >  2:2617521/1‑151 (MQ=255)
                                                                                     cAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                 >  2:118243/1‑151 (MQ=255)
                                                                                       aTCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                               >  1:2052625/1‑151 (MQ=255)
                                                                                            caaaaCTCCATCCCGCACTGGGTCAAGGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                             <  1:3672697/145‑1 (MQ=255)
                                                                                               aaCTTCACCCCGCACTGGGTCAAGGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt                                          <  1:2793426/148‑1 (MQ=255)
                                                                                                              tGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                        >  1:1619098/1‑151 (MQ=255)
                                                                                                               gggTCAAAGGGTAGCAAGACGGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                            <  1:2263303/146‑1 (MQ=255)
                                                                                                               gagTCAAAGGAAAGAGACACCGCGGCGCGACCGAGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt                       <  1:3701034/149‑1 (MQ=255)
                                                                                                                  tCAAAGGGAAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCTCGATAAAGAGTTTCTCGTGCAATAAttt                       >  2:3264856/1‑148 (MQ=255)
                                                                                                                    aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                   >  2:3440933/1‑150 (MQ=255)
                                                                                                                    aaaGGGTAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                  <  1:1819871/151‑1 (MQ=255)
                                                                                                                              aaaaCTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGc        >  2:707125/4‑151 (MQ=255)
                                                                                                                                gACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc       <  1:768565/150‑1 (MQ=255)
                                                                                                                                      ggcgTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATg  <  1:4061459/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
CTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATG  >  NC_000913/3177600‑3177882

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