Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC50.0% 33.8 / 12.4 34G397G (GGT→GGG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (15/2);  new base C (0/17);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGC  >  NC_000913/1969647‑1969944
                                                                                                                                                    |                                                                                                                                                      
ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCt                                                                                                                                                     <  1:145076/151‑1 (MQ=255)
       gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACt                                                                                                                                               >  1:138388/1‑150 (MQ=255)
                    cGATCATGGTCGCGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGACTCTTCGATCACCCTTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACg                                                                                                                                 <  2:407443/151‑1 (MQ=255)
                            gTCGCGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACCCGATTGCCTGCCTCTTCGATCACCCTTTT‑GATCTCTTTTGCTCCCTGCGCGCTACGGCTGGCAAGTTTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGct                                                                                                                         <  2:54736/151‑1 (MQ=255)
                             tCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctc                                                                                                                        >  2:367225/1‑151 (MQ=255)
                                 gccgccgTTATCCCCAAGTTTCGAACCCCGCTGGACACTATTGACTGACTCTCTGATCGCCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCGGGAAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACg                                                                                                                                 <  1:638437/133‑1 (MQ=255)
                                        ttATTCCCCATTTTCAACCCCGGCTGGACACGATTGCCTACCTCTTCGATCACCCCTGG‑GCCCTTTTTTGCTGCCTGCCCGCTGCGGCGGGCAAGTGTGCGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                               <  2:525033/149‑1 (MQ=255)
                                         tATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc                                                                                                            >  1:119102/1‑151 (MQ=255)
                                         tATTCACCAGTTTCGAACCCTGCTGGACACGAATGACCGCCTCTTCCATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGCTTGGCAAGATTGCGCACTTCATCTGCCACTACAGCAAAACAACGCTCCTGATCTGCGGCCCGTGc                                                                                                            >  1:206420/1‑151 (MQ=255)
                                                    tccGACCCCTGCTGAACACGTTTGACTGACTCTCCGATCAGCCCTTT‑GATCTCTTTTCCTCCCTGCGCCCTACGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACt                                                                                                  <  2:397306/148‑1 (MQ=255)
                                                        aaCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                             >  2:624695/1‑151 (MQ=255)
                                                         aCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                             >  1:351913/1‑150 (MQ=255)
                                                                 gaacacGTTTGACTGCCTCTTCGATCGCCCTTTT‑TATCTCTTTTGCTCCCGGCGCGCGGCGGCTGGCAAGTTTCCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg                                                                                     <  1:805943/148‑1 (MQ=255)
                                                                         ttGACTTCCTCTTCGACCACCCCTTT‑GATCTCTTTTGCTCCCCGCGCGCTCCGGCTGGAAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                                <  2:90831/147‑1 (MQ=255)
                                                                         ttGACTGCCTCTTCTATCCGCCCTTT‑GTTCTTTTTTGCTCCCGGCGCGCTGCGGCGGGAAAGTTGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                                      <  1:350981/141‑1 (MQ=255)
                                                                         ttGACTGCCTCTTCGATCACCCCTTT‑GATCTCTTTTGCTCCCGGCGCCCTCCGGCGGGAAAGTTGGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                              <  2:21843/149‑1 (MQ=25)
                                                                                       gATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                               >  1:493046/1‑150 (MQ=255)
                                                                                        aTCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                              >  2:446299/1‑150 (MQ=255)
                                                                                          acacccTTTGGATCTTTTTCGCTCCCTGGCCGCTCCGGCTGGCAATATGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                                  <  2:743806/142‑1 (MQ=18)
                                                                                            ccctttt‑GCTCTCTTTTGCCCCCGGCCCCCGCCGCCGGGCAAGTTTGCGCCTTCCCCCTGCCACTGCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                                         <  1:639002/147‑1 (MQ=21)
                                                                                             cccTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTcc                                                         >  1:403823/1‑150 (MQ=255)
                                                                                             cccTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCg                                                        >  1:3770/1‑151 (MQ=255)
                                                                                             cccTTT‑ACTCTCTTTTGCTCCCGGCGCGCTACGGCGGCAAAGTTGGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                                    <  2:629743/139‑1 (MQ=21)
                                                                                                ttt‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                     >  2:49273/1‑151 (MQ=255)
                                                                                                    gATCTGTTTGGCCGCCGCGGCGCTCCGCGTGGCAAGATTGCGCTCTCCCCCTGCCACTGCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAGGCAATTCCGt                                                       <  2:164925/146‑1 (MQ=21)
                                                                                                          ttttTCTCCCTGCGCGCTCCGGCTGCAAAGTTGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                                 <  2:799101/146‑1 (MQ=25)
                                                                                                          ttttGCTGCCTGCGCGCTCCGGCGGGAAAGTTTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                               <  2:531172/148‑1 (MQ=37)
                                                                                                          ttttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                            >  2:782821/1‑151 (MQ=255)
                                                                                                                        cccgaCGGCGGCCAATATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAt                              <  2:803606/147‑1 (MQ=25)
                                                                                                                          cggcggcggGAAAGTTTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATaa                                        <  1:651826/136‑1 (MQ=25)
                                                                                                                            aCGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt                          >  1:91387/1‑151 (MQ=255)
                                                                                                                                  gccAAAATGGCCCACTCCCCCTGCCCCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                       <  1:247086/130‑1 (MQ=25)
                                                                                                                                     aaGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt                  >  1:562332/1‑150 (MQ=255)
                                                                                                                                                 ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGATCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTcc     >  1:522154/1‑151 (MQ=255)
                                                                                                                                                    aCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGc  >  2:454682/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGC  >  NC_000913/1969647‑1969944

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