Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC40.8% 76.7 / 10.3 49G397G (GGT→GGGtapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base C (1/19);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 3.98e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC  >  NC_000913/1969647‑1969941
                                                                                                                                                     |                                                                                                                                                   
ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCt                                                                                                                                                  <  2:287246/151‑1 (MQ=255)
    aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCa                                                                                                                                              <  1:1098446/151‑1 (MQ=255)
             aCAATATCGATCATGGTCGGGGCGGTGTTATTCACCAGTTCCGACCCCTGCGGCACACGATTGACTGCCTCTCCGATCACCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAAGGTTGCGCACTCCACCTGCCACTACCGCaa                                                                                                                                     <  1:456199/151‑1 (MQ=255)
              cAATATCGATCATGGTCGCGGCGGCTTTATTCACCAGTTCCGACCCCTGCTGGACACGATTGACTGCCTCTTCGATCAGCCCTTTGCTCTCTTTTGCTGCCTGCGCGCTGCGGCGGGCAAGATTGCGCCCTTCCCCTGCCACTa                                                                                                                                           <  1:181600/144‑1 (MQ=255)
                      aTCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGcc                                                                                                                            >  2:963057/1‑151 (MQ=255)
                        cATGGTCGCGGCGTCTTTATTCCCCAGTTTCGACCCCGGCTGGACACGTTTGCCTGCCTCTTCGATCCGCCTTTGGCTCTCTTTCGCTCCCTGCGCGCTACGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCAc                                                                                                                               <  1:128849/146‑1 (MQ=255)
                         aTGGTCGCGGCGGCGCTATTCACCAGTTTCGAACCCTGCTGGACACGCTTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGcccc                                                                                                                          >  1:895305/1‑150 (MQ=255)
                            gTCGGGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGACTCTCCGATCAGCCCTTTGCTCTTTTTGCCTCCCGGCGCTCTACGGCGGGCAAGTTGGCGCTCTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                        <  1:417791/149‑1 (MQ=255)
                              cgcgGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACCGACTCTTCGCTCAGCCCTTTGATCTCTTTTGCTCCCTCCGCGCTACGGCCGGCACGATTGCGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGccc                                                                                                                     >  1:441312/1‑148 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                                   >  2:820823/1‑151 (MQ=255)
                                   cggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTAGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                               >  1:264237/1‑151 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                            >  1:157822/1‑149 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                              >  1:611574/1‑147 (MQ=255)
                                           ttCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGCTTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc                                                                                                         >  1:690634/1‑149 (MQ=255)
                                                 aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                  >  2:132259/1‑150 (MQ=255)
                                                   tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCAc                                                                                                >  1:871258/1‑150 (MQ=255)
                                                         aCCCTCTTGGCCACCATTGTCTGCCTTTACTACCCGCCCTTTGGTGTCTTTTCCTTCTTGCGCCCACCGCCGGCAAGGTTTGCGCTTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGTCGCTTCCACTGCCg                                                                                           <  2:446870/149‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGCTTGCGCACTTCACCTGCCACTACAGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCAACTGCCGCATTc                                                                                      >  1:962107/1‑149 (MQ=255)
                                                               ctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                                   >  2:107255/1‑151 (MQ=255)
                                                                 gaccacttttGCCTGCCTCTCCGACCACCCCTTTGATCTCTTTTGCTCCCGGCCCGCAACGGCGGCAATGTTGCCCCACTTCCCCTGCCCCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                      <  1:747209/138‑1 (MQ=255)
                                                                 gaacacttttGACTGCCTCTTCGATCACCCTTTTGATCTTTTTTGCGCCCTGCGCGCTCCGCCGGGAAAGATTGCGCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg                                                                                  <  2:802888/142‑1 (MQ=255)
                                                                        aTTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCGGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                             <  2:1083750/148‑1 (MQ=255)
                                                                         ttGACTGACTTTTCGATCACCCTTTTGATCTTTTTTGCTCCCGGCGCGCTACGGGGGGCAAGATTGCACACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa                                                                            <  2:1000643/148‑1 (MQ=255)
                                                                         ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                                         >  2:59464/1‑151 (MQ=255)
                                                                              tgtctctTCGACCCGCCCTTTGCTCCTTTTTGCCCCCCGGGCTCTACGGCGGGAAAGATGGCGCACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                       <  2:96350/145‑1 (MQ=21)
                                                                              tgcctctTCGACCACCCTTTTTATCTTTTTTGCTCCCTGCGCGCTACGGCTGGCAAGTTTGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                                        <  2:883278/144‑1 (MQ=25)
                                                                                 ctctTCGATCCCCCCTTTGATCTCTTTTTCTGCCGCCGCGCTACGGCGGGAAAGATTGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCtt                                                                  <  1:791046/150‑2 (MQ=25)
                                                                                    ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                               >  1:921393/1‑150 (MQ=255)
                                                                                    tccGATCACCCTTTTATTCTTTTTTGCTGCCTGCGCCCTCCGGCTGCAAAGTTGGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                                    <  2:706688/143‑1 (MQ=21)
                                                                                       gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                            >  1:858138/1‑150 (MQ=255)
                                                                                        aTCACCCCTTTGATCTTTTGTGCTCCCGGCGCGCTACGGCGGGAAAGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                           <  2:848152/150‑1 (MQ=25)
                                                                                          ccgCCTTTTGATCGCTTTGCCGGCCTGCGCGCTGGGGGGGGAAGGTTGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                                        <  1:1024637/149‑1 (MQ=18)
                                                                                                tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                  >  2:434953/1‑151 (MQ=255)
                                                                                                tttGATCTCTTTTGCTCCCGGCGCGCTCCGGCGGGCAAGATTGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                   <  2:634804/150‑1 (MQ=37)
                                                                                                     tcttTTTTGCTCCCTGCGCGCTACGGCTGGAAAGATTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                              <  2:817989/150‑1 (MQ=37)
                                                                                                       ttttttGCTGCCTGCGCGCTCGGGCTGGCAAGTTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                 <  2:292604/143‑1 (MQ=37)
                                                                                                       tctTTTCCTCCCGGCGCGCTCCGGCGGGAAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                              <  2:762698/148‑1 (MQ=25)
                                                                                                        ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                          >  1:1050037/1‑151 (MQ=255)
                                                                                                         ttttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                         >  2:970191/1‑151 (MQ=255)
                                                                                                          tttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                         >  1:470351/1‑150 (MQ=255)
                                                                                                           ttcctgccggcGCGCTCCGGCGGGCAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa                                        <  2:390227/147‑1 (MQ=37)
                                                                                                           ttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGAt                                       >  1:50719/1‑151 (MQ=255)
                                                                                                                     cgcgcTCCGGCTGCCAAGATTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                                <  2:567191/148‑1 (MQ=255)
                                                                                                                          tACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                        >  1:69045/1‑151 (MQ=255)
                                                                                                                            cGGCTGGAAAGTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                        <  2:819186/149‑1 (MQ=25)
                                                                                                                                   cAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt               >  1:860424/1‑151 (MQ=255)
                                                                                                                                             cgcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg        >  2:124251/1‑150 (MQ=255)
                                                                                                                                               cACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtc     >  2:362177/1‑151 (MQ=255)
                                                                                                                                                  ttCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTcc  >  1:634698/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC  >  NC_000913/1969647‑1969941

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: