Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,969,7940AC45.1% 68.1 / 19.6 52G397G (GGT→GGG) tapmethyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (24/4);  new base C (0/23);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.27e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942
                                                                                                                                                     |                                                                                                                                                    
ggTGACCGATCTGACAATATCCATCATGGTGGCGGCGGCGTTATTCGCCAGTTTCCAACCCTGCTGGACACGATTGACGCACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGTTACGGCTGGCAAGA‑TTGCGCACTTCACCt                                                                                                                                                   <  1:140728/151‑1 (MQ=255)
    aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCa                                                                                                                                               <  1:140117/151‑1 (MQ=255)
         aCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTAcc                                                                                                                                          >  2:55073/1‑151 (MQ=255)
            gACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGc                                                                                                                                        >  2:173728/1‑150 (MQ=255)
              cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCaa                                                                                                                                      >  2:478508/1‑150 (MQ=255)
                 tatCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAAcc                                                                                                                                  >  1:397450/1‑151 (MQ=255)
                           ggccgggggggcgTTTTTCCCCAGTTTCAACCCCTGCTGGACACGTTTGACTGCCTCTTCGATCACCCTTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:949892/140‑1 (MQ=255)
                           ggTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                         >  2:573902/1‑150 (MQ=255)
                            gTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGc                                                                                                                        >  1:765177/1‑150 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  2:213753/150‑1 (MQ=255)
                                cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                                    >  2:283287/1‑150 (MQ=255)
                                 gggtgctttATTCCCCATTTTCGACCCCTGCTGGACTCGATTGACTGACTCTTGGTTCAGCCTTTA‑GATCTCTTTGGCTGCCTGCGCTCTACGGCGGGCAAGA‑TTGCCCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                                   <  2:963020/143‑1 (MQ=255)
                                    ggcgTTATTCCCCAGTTTCGAACCCGGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTTTTTGGCTGCCGGCGCGCTACGGCTGGCAAGA‑TTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                                <  2:597982/150‑1 (MQ=255)
                                     gcgTTATTCACCAGTTTCGAACCCTGCTGGACACGTTTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcgc                                                                                                              <  2:785706/151‑1 (MQ=255)
                                          aTTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                         >  2:1008589/1‑151 (MQ=255)
                                                 aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                   >  1:33839/1‑150 (MQ=255)
                                                       gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                            >  2:331582/1‑151 (MQ=255)
                                                         aCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                           >  1:702018/1‑150 (MQ=255)
                                                          cccTGTTGGACACGTTTGACTGCCTCTTCGATCAGCCTTTT‑GTTCTTTTTTGCTCCCTGCGCGCTCCGGCGGGAAAGA‑TGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                         <  2:419398/151‑1 (MQ=255)
                                                          cccTGCCGAACTGGATTGCCTGCCTCTTCACTCAGCCTTTT‑TTTCTTTTTTGCCGCGTGCCCGTTCCGGTGGGCAATC‑TTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                           <  1:1058057/149‑1 (MQ=255)
                                                             tgcggGACACGTTTACCTCCCTCTCCGACCGCCCTTTT‑GTCCTTTTTTGCGCCCTGCGCGCTGCGGGGGGAAAGT‑TGCCCCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                        <  1:464256/149‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       >  1:120144/1‑149 (MQ=255)
                                                                 ggACACGTTTGCCTGCCTCTTCGATCACCCTTTT‑GATCTCTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGA‑TGGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  2:897243/146‑1 (MQ=255)
                                                                 gaacccgtttGACTGCCTCTCCGATCGCCCTTTT‑GATCTTTTTTGCTGCCGGCGCGCAAGGGCGGGAAAGT‑TGCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  1:171598/138‑1 (MQ=255)
                                                                   acacTATTGACTGCCTCTCCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCGGGCAATT‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  2:713273/144‑1 (MQ=255)
                                                                    cagtaTTGCCTGACTCTCCGATCACCCTTTT‑GTTCTCTTTTGCTCCCGGCGCGCTCCGGCTGGCAAGA‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg                                                                               <  2:849283/147‑1 (MQ=255)
                                                                       gATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                            >  1:638368/1‑151 (MQ=255)
                                                                         ttGAGTGCCTCTAGGATCACCCCTTT‑GCTCTTTTTTGCGCCCTGCGCGCTCCGGCTGCAAAGA‑TTGCGCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                              <  2:721970/147‑1 (MQ=255)
                                                                                     tCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                                >  1:824316/1‑149 (MQ=255)
                                                                                       gATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                             >  2:223674/1‑150 (MQ=255)
                                                                                        aTCAGCCCTTT‑GATCTCTTTTGCTGCCCGCGCGCTACGGCGGGCAAGA‑TGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              <  2:12820/148‑1 (MQ=37)
                                                                                         tCAGCCTTTT‑GATCTCTTTCCCTGCCTGCGCCCTAGGCATGGCAAGT‑TGGCCCACTTCGCCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa                                                                 <  2:413615/144‑1 (MQ=25)
                                                                                         accaccctttggctcttTTTTGCGCCCCGCGCGCTCCGGTGGGCAAGA‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa                                                                 <  2:1104192/132‑1 (MQ=21)
                                                                                           aGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                                        >  2:865426/1‑151 (MQ=255)
                                                                                             cccTTTGGATCTTTTTTCCTGCCTGCGCGCTCGGGCTGGCAAGA‑TTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt                                                          <  2:608172/148‑1 (MQ=25)
                                                                                                ttt‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                   >  1:158737/1‑151 (MQ=255)
                                                                                                    gttcttTTTTGCGCCCGGCGCGCTCGGGCTGGAAAGA‑TTCCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  2:579209/149‑1 (MQ=25)
                                                                                                     aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                >  2:451890/1‑150 (MQ=255)
                                                                                                        ttttttGCTCCCTGCGCGCTCCGGCTGGAAAGA‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  2:30282/145‑1 (MQ=37)
                                                                                                        tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                                            >  1:923223/1‑151 (MQ=255)
                                                                                                         ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                           >  1:244336/1‑151 (MQ=255)
                                                                                                          ttttGCTGCCGCCGCGCTGCGGCGGACAAGA‑TGGCGCCTTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  1:738255/145‑1 (MQ=25)
                                                                                                           tttGCTGCCTGCGCGCTACGGCTGGCAGGA‑TTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                           <  1:808142/149‑1 (MQ=255)
                                                                                                                  cctgcGCGCTCCGGCTGGCAAGT‑TTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                                    <  2:951642/149‑1 (MQ=255)
                                                                                                                  cctgcGCGCTCCGGC‑GGCAAAATTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                                   <  1:832537/150‑1 (MQ=25)
                                                                                                                    tgcGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                                >  1:1103486/1‑151 (MQ=255)
                                                                                                                             ccggcggacaaattggcgcACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAttttt                        <  1:895588/136‑1 (MQ=25)
                                                                                                                                          ttGCCCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtggt             <  1:9241/149‑1 (MQ=255)
                                                                                                                                          tgccccccTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtgg              <  1:429730/141‑1 (MQ=255)
                                                                                                                                            gcgcACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCg         >  1:109339/1‑151 (MQ=255)
                                                                                                                                                aCTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAtct     >  1:282980/1‑151 (MQ=255)
                                                                                                                                                   tCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCt  >  1:764358/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTT‑GATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGA‑TTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942

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