Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC57.9% 24.1 / 48.2 38G397G (GGT→GGG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (13/3);  new base C (0/22);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 1.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-05

CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969650‑1969942
                                                                                                                                                 |                                                                                                                                                    
cACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGc                                                                                                                                                 >  2:794275/1‑150 (MQ=255)
       cTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGAAAGATTGCGCACTTCACCTGCCACTAcc                                                                                                                                          <  2:645047/150‑1 (MQ=255)
        tGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGc                                                                                                                                        >  1:812619/1‑151 (MQ=255)
             atatCGATCATGGTCGCGGCGGTTTTATCCACCAGTTTCGACCCCTGCGGGACACGATTGATGGCCTCTTCGATCAGCCCTTT‑GATCTCTTTTGCGGCCTGCGCGCTACGGCTGGAAAGTTTGCGCACTTCCCCTGCCACTACCGCaa                                                                                                                                      <  1:272217/148‑1 (MQ=255)
              tatCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAAcc                                                                                                                                  >  2:937547/1‑151 (MQ=255)
                 cGATCATGGTCGCGGCGGCTTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTCCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCAc                                                                                                                                <  2:256701/150‑1 (MQ=255)
                 cGATCATGGCCGCGGCGGCTTTATCCACCAGTTTCGACCCCTGCTGGACCCTATTGACTGCCTCTCCGTTCACTTTTTT‑GACCTCTTCTCCCCCCACCGCGCTCCCCCGGAAAAAATTCCGCCCTTCCCCTGCCACTACCGCAAAACCACg                                                                                                                               <  1:23973/151‑1 (MQ=255)
                   aTCATGGACGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGcc                                                                                                                             >  2:475702/1‑151 (MQ=255)
                    tCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGccc                                                                                                                            >  1:570009/1‑151 (MQ=255)
                                    gTTTTTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGCCTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTCCCTGCGCGCTCCGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                             <  2:624435/150‑1 (MQ=255)
                                     tttttCCCCATTTTCGACCCCTGCTGGCCACTATTGACGGACTCTTCGATCACCCCTTT‑GATCTTTTTTGCTCCCTGCGCGCTCCGGCTGGCAGGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcgc                                                                                                              <  1:227153/145‑1 (MQ=255)
                                     ttATTCCCCTGTTTCAACCCCGGCTGGACACGATTGACTGACTCTTCGATCACCCCTTT‑GATCTTTTTTGCTCCCGGCGCGCTACGGCTGCAAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  1:518339/141‑1 (MQ=255)
                                          cACCAGTTCCGAACCCTGCTGGACACGTTTGACTGACTCTTCGATCACCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATGCCGCACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCg                                                                                                        <  1:450148/149‑1 (MQ=255)
                                              aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                   >  2:58360/1‑150 (MQ=255)
                                                           gctgGACATGATTGACTGACTCTTCGATCCCCCCTTT‑GATCTCTTTTGCTCCCGGCGCGCTGCGGCGGGAAGGATGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                        <  2:869808/148‑1 (MQ=255)
                                                             gggACACGTTTGACTGACTCTTCGATCAGCCTTTT‑GATCTTTTTTGCTCCCTGCGCGCTCCGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg                                                                                   <  2:944576/150‑1 (MQ=255)
                                                                   cGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                              >  2:853391/1‑150 (MQ=255)
                                                                      ttgccggactCTCCGTCCAGCCCTTG‑GACCTTGTTGGCTGCCGGCGGGCGCGGGCGGGAAAGATTGCCCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa                                                                             <  1:383707/142‑1 (MQ=16)
                                                                               ttttcgtccaccccTTT‑GATCTCTTTTGCTCCCTGCGCTCTACGGCTGCCAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCt                                                                    <  2:828953/137‑1 (MQ=25)
                                                                               ctctcggtccacccTTTGGATCTCTTTCGCTCCCTGCGCGCAGCGGCTGCCAATATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                                     <  1:384697/137‑1 (MQ=21)
                                                                                 ttataccaccctttt‑GATTTTTTTTGCTGCCGGCGCTTTAGGGCTGGAAAGTTTTCGCTCTTCCCCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                                         <  2:107591/129‑1 (MQ=21)
                                                                                 tccGACCAGCCTTTT‑GTTCTCTTTTGCTCCCTGCGCGCTACGGCTGGAAAGATTGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                                     <  2:733470/143‑1 (MQ=25)
                                                                                       cAGCCCTTT‑GATCTCTTTTGATGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCCGGGCCAGAATATTGGTCTGGAAAGCAAt                                                          >  2:790258/1‑150 (MQ=255)
                                                                                        aGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                                        >  2:853351/1‑151 (MQ=255)
                                                                                           ccTTT‑GATCTCTTTTGCTCCCGGCGCGCTCGGGCGGGAAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                                         <  2:86864/147‑1 (MQ=25)
                                                                                             ttt‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                   >  1:429967/1‑151 (MQ=255)
                                                                                                  aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                >  1:701519/1‑150 (MQ=255)
                                                                                                  aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                >  1:337645/1‑150 (MQ=255)
                                                                                                               cttgcGCGCGGCGGCGGCCAAAATTGCCCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                     <  2:120753/146‑1 (MQ=25)
                                                                                                               cctgcGCGCTCCGGCGGCCAAGATGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                  <  2:986278/151‑1 (MQ=255)
                                                                                                               ccggcgcgcTAGGGCTGGAAAGTTGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATa                                       <  1:379255/143‑1 (MQ=25)
                                                                                                               ccggccccttccGGGTGGCAAGTTCGCCCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                     <  2:744418/137‑1 (MQ=25)
                                                                                                                     cgcTCCGGCGGCCAAGTTTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAt                            <  2:71306/151‑1 (MQ=25)
                                                                                                                     cgcTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAt                            >  1:29903/1‑151 (MQ=255)
                                                                                                                        tACGGCTGGAAAGATTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtt                           <  2:687599/149‑1 (MQ=255)
                                                                                                                           ggCTGGAAAGATTGCGCCCTCCCCCTGCCCCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTg                      <  2:910284/151‑1 (MQ=255)
                                                                                                                                  aaGATTGCGCCTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt                <  2:498961/150‑1 (MQ=255)
                                                                                                                                                cACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt  <  2:867327/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969650‑1969942

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