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→GGG) tapmethyl‑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

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC  >  NC_000913/1969647‑1969941
                                                                                                                                                     |                                                                                                                                                   
ggTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCt                                                                                                                                                  <  2:287246/151‑1 (MQ=255)
    aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCa                                                                                                                                              <  1:1098446/151‑1 (MQ=255)
             aCAATATCGATCATGGTCGGGGCGGTGTTATTCACCAGTTCCGACCCCTGCGGCACACGATTGACTGCCTCTCCGATCACCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAA‑GGTTG‑CGCACTCCACCTGCCACTACCGCaa                                                                                                                                     <  1:456199/151‑1 (MQ=255)
              cAATATCGATCATGGTCGCGGCGGCTTTATTCACCAGTTCCGACCCCTGCTGGACACGATTGACTGCCTCTTCGATCAGCCCTTTGCTCTCTTTTGCTGCCTGCGCGCTGCGGCGGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTa                                                                                                                                           <  1:181600/144‑1 (MQ=255)
                      aTCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGcc                                                                                                                            >  2:963057/1‑151 (MQ=255)
                        cATGGTCGCGGCGTCTTTATTCCCCAGTTTCGACCCCGGCTGGACACGTTTGCCTGCCTCTTCGATCCGCCTTTGGCTCTCTTTCGCTCCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCAc                                                                                                                               <  1:128849/146‑1 (MQ=255)
                         aTGGTCGCGGCGGCGCTATTCACCAGTTTCGAACCCTGCTGGACACGCTTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGcccc                                                                                                                          >  1:895305/1‑150 (MQ=255)
                            gTCGGGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACACGTTTGACTGACTCTCCGATCAGCCCTTTGCTCTTTTTGCCTCCCGGCGCTCTACGGCGGGCAA‑GTTGG‑CGCTCTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                        <  1:417791/149‑1 (MQ=255)
                              cgcgGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACCGACTCTTCGCTCAGCCCTTTGATCTCTTTTGCTCCCTCCGCGCTACGGCCGGCAC‑GATTG‑CGCCCTTCACCTGCCACTACCGCAAAACCACGCCCCTGccc                                                                                                                     >  1:441312/1‑148 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                                   >  2:820823/1‑151 (MQ=255)
                                   cggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTAGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                               >  1:264237/1‑151 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                            >  1:157822/1‑149 (MQ=255)
                                        ttATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                              >  1:611574/1‑147 (MQ=255)
                                           ttCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GCTTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGc                                                                                                         >  1:690634/1‑149 (MQ=255)
                                                 aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                  >  2:132259/1‑150 (MQ=255)
                                                   tttCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCAc                                                                                                >  1:871258/1‑150 (MQ=255)
                                                         aCCCTCTTGGCCACCATTGTCTGCCTTTACTACCCGCCCTTTGGTGTCTTTTCCTTCTTGCGC‑CCACCGCCGGCAAGGTTTG‑CGCTTTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGTCGCTTCCACTGCCg                                                                                           <  2:446870/149‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GCTTG‑CGCACTTCACCTGCCACTACAGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCAACTGCCGCATTc                                                                                      >  1:962107/1‑149 (MQ=255)
                                                               ctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAgg                                                                                   >  2:107255/1‑151 (MQ=255)
                                                                 gaccacttttGCCTGCCTCTCCGACCACCCCTTTGATCTCTTTTGCTCCCGGCCCGCAACGGC‑GGCAA‑TGTTGCCCCACTTCCCCTGCCCCTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                      <  1:747209/138‑1 (MQ=255)
                                                                 gaacacttttGACTGCCTCTTCGATCACCCTTTTGATCTTTTTTGCGCCCTGCGCGCTCCGCCGGGAAA‑GATTG‑CGCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAggg                                                                                  <  2:802888/142‑1 (MQ=255)
                                                                        aTTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTTGCTGCCTGCGCGCTACGGCGGGCAA‑GATTG‑CGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                             <  2:1083750/148‑1 (MQ=255)
                                                                         ttGACTGACTTTTCGATCACCCTTTTGATCTTTTTTGCTCCCGGCGCGCTACGGGGGGCAA‑GATTG‑CACACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGaa                                                                            <  2:1000643/148‑1 (MQ=255)
                                                                         ttGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAatat                                                                         >  2:59464/1‑151 (MQ=255)
                                                                              tgtctctTCGACCCGCCCTTTGCTCCTTTTTGCCCCCCGGGCTCTACGGCGGGAAA‑GATGG‑CGCACTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                       <  2:96350/145‑1 (MQ=21)
                                                                              tgcctctTCGACCACCCTTTTTATCTTTTTTGCTCCCTGCGCGCTACGGCTGGCAA‑GTTTG‑CCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATAtt                                                                        <  2:883278/144‑1 (MQ=25)
                                                                                 ctctTCGATCCCCCCTTTGATCTCTTTTTCTGCCGCCGCGCTACGGCGGGAAA‑GATTG‑CGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCtt                                                                  <  1:791046/150‑2 (MQ=25)
                                                                                    ttCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                               >  1:921393/1‑150 (MQ=255)
                                                                                    tccGATCACCCTTTTATTCTTTTTTGCTGCCTGCGCCCTCCGGCTGCAAA‑GTTGG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTc                                                                    <  2:706688/143‑1 (MQ=21)
                                                                                       gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                            >  1:858138/1‑150 (MQ=255)
                                                                                        aTCACCCCTTTGATCTTTTGTGCTCCCGGCGCGCTACGGCGGGAAA‑GTTTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                           <  2:848152/150‑1 (MQ=25)
                                                                                          ccgCCTTTTGATCGCTTTGCCGGCCTGCGCGCTGGGGGGGGAAG‑GTTGG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                                        <  1:1024637/149‑1 (MQ=18)
                                                                                                tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                  >  2:434953/1‑151 (MQ=255)
                                                                                                tttGATCTCTTTTGCTCCCGGCGCGCTCCGGCGGGCAA‑GATTG‑CCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTc                                                   <  2:634804/150‑1 (MQ=37)
                                                                                                     tcttTTTTGCTCCCTGCGCGCTACGGCTGGAAA‑GATTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                              <  2:817989/150‑1 (MQ=37)
                                                                                                       ttttttGCTGCCTGCGCGCTCGGGCTGGCAA‑GTTTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGa                                                 <  2:292604/143‑1 (MQ=37)
                                                                                                       tctTTTCCTCCCGGCGCGCTCCGGCGGGAAA‑GATTG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATaa                                              <  2:762698/148‑1 (MQ=25)
                                                                                                        ctTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                          >  1:1050037/1‑151 (MQ=255)
                                                                                                         ttttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                         >  2:970191/1‑151 (MQ=255)
                                                                                                          tttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                         >  1:470351/1‑150 (MQ=255)
                                                                                                           ttcctgccggcGCGCTCCGGCGGGCAA‑GATTG‑CGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa                                        <  2:390227/147‑1 (MQ=37)
                                                                                                           ttGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGAt                                       >  1:50719/1‑151 (MQ=255)
                                                                                                                     cgcgcTCCGGCTGCCAA‑GATTG‑CGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                                <  2:567191/148‑1 (MQ=255)
                                                                                                                          tACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                        >  1:69045/1‑151 (MQ=255)
                                                                                                                            cGGCTGGAAA‑GTTTG‑CCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttt                        <  2:819186/149‑1 (MQ=25)
                                                                                                                                   cAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCt               >  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)
                                                                                                                                                     |                                                                                                                                                   
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAA‑GATTG‑CGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCC  >  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: ‑