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→GGGtapmethyl‑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

GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  NC_000913/1969647‑1969942
                                                                                                                                                     |                                                                                                                                                    
ggTGACCGATCTGACAATATCCATCATGGTGGCGGCGGCGTTATTCGCCAGTTTCCAACCCTGCTGGACACGATTGACGCACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGTTACGGCTGGCAAGATTGCGCACTTCACCt                                                                                                                                                   <  1:140728/151‑1 (MQ=255)
    aCCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCa                                                                                                                                               <  1:140117/151‑1 (MQ=255)
         aCTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTAcc                                                                                                                                          >  2:55073/1‑151 (MQ=255)
            gACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGc                                                                                                                                        >  2:173728/1‑150 (MQ=255)
              cAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCaa                                                                                                                                      >  2:478508/1‑150 (MQ=255)
                 tatCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAAcc                                                                                                                                  >  1:397450/1‑151 (MQ=255)
                           ggccgggggggcgTTTTTCCCCAGTTTCAACCCCTGCTGGACACGTTTGACTGCCTCTTCGATCACCCTTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                          <  2:949892/140‑1 (MQ=255)
                           ggTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTg                                                                                                                         >  2:573902/1‑150 (MQ=255)
                            gTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGc                                                                                                                        >  1:765177/1‑150 (MQ=255)
                               gcggcggcgTTATTCACCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGctct                                                                                                                     <  2:213753/150‑1 (MQ=255)
                                cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCtctc                                                                                                                    >  2:283287/1‑150 (MQ=255)
                                 gggtgctttATTCCCCATTTTCGACCCCTGCTGGACTCGATTGACTGACTCTTGGTTCAGCCTTTAGATCTCTTTGGCTGCCTGCGCTCTACGGCGGGCAAGATTGCCCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                                   <  2:963020/143‑1 (MQ=255)
                                    ggcgTTATTCCCCAGTTTCGAACCCGGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTTTTTGGCTGCCGGCGCGCTACGGCTGGCAAGATTGCGCACTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                                <  2:597982/150‑1 (MQ=255)
                                     gcgTTATTCACCAGTTTCGAACCCTGCTGGACACGTTTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcgc                                                                                                              <  2:785706/151‑1 (MQ=255)
                                          aTTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                         >  2:1008589/1‑151 (MQ=255)
                                                 aGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTcc                                                                                                   >  1:33839/1‑150 (MQ=255)
                                                       gAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                            >  2:331582/1‑151 (MQ=255)
                                                         aCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                           >  1:702018/1‑150 (MQ=255)
                                                          cccTGTTGGACACGTTTGACTGCCTCTTCGATCAGCCTTTTGTTCTTTTTTGCTCCCTGCGCGCTCCGGCGGGAAAGATGGCGCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                         <  2:419398/151‑1 (MQ=255)
                                                          cccTGCCGAACTGGATTGCCTGCCTCTTCACTCAGCCTTTTTTTCTTTTTTGCCGCGTGCCCGTTCCGGTGGGCAATCTTGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGc                                                                                           <  1:1058057/149‑1 (MQ=255)
                                                             tgcggGACACGTTTACCTCCCTCTCCGACCGCCCTTTTGTCCTTTTTTGCGCCCTGCGCGCTGCGGGGGGAAAGTTGCCCCATTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                        <  1:464256/149‑1 (MQ=255)
                                                              gctgGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       >  1:120144/1‑149 (MQ=255)
                                                                 ggACACGTTTGCCTGCCTCTTCGATCACCCTTTTGATCTCTTTTGCTGCCTGCGCGCTCCGGCTGGCAAGATGGCCCACTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  2:897243/146‑1 (MQ=255)
                                                                 gaacccgtttGACTGCCTCTCCGATCGCCCTTTTGATCTTTTTTGCTGCCGGCGCGCAAGGGCGGGAAAGTTGCCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  1:171598/138‑1 (MQ=255)
                                                                   acacTATTGACTGCCTCTCCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCGGGCAATTTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTc                                                                                       <  2:713273/144‑1 (MQ=255)
                                                                    cagtaTTGCCTGACTCTCCGATCACCCTTTTGTTCTCTTTTGCTCCCGGCGCGCTCCGGCTGGCAAGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAg                                                                               <  2:849283/147‑1 (MQ=255)
                                                                       gATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAat                                                                            >  1:638368/1‑151 (MQ=255)
                                                                         ttGAGTGCCTCTAGGATCACCCCTTTGCTCTTTTTTGCGCCCTGCGCGCTCCGGCTGCAAAGATTGCGCTCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGa                                                                              <  2:721970/147‑1 (MQ=255)
                                                                                     tCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                                >  1:824316/1‑149 (MQ=255)
                                                                                       gATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                             >  2:223674/1‑150 (MQ=255)
                                                                                        aTCAGCCCTTTGATCTCTTTTGCTGCCCGCGCGCTACGGCGGGCAAGATGGCGCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                              <  2:12820/148‑1 (MQ=37)
                                                                                         tCAGCCTTTTGATCTCTTTCCCTGCCTGCGCCCTAGGCATGGCAAGTTGGCCCACTTCGCCTGCCACTCCCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa                                                                 <  2:413615/144‑1 (MQ=25)
                                                                                         accaccctttggctcttTTTTGCGCCCCGCGCGCTCCGGTGGGCAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGa                                                                 <  2:1104192/132‑1 (MQ=21)
                                                                                           aGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                                        >  2:865426/1‑151 (MQ=255)
                                                                                             cccTTTGGATCTTTTTTCCTGCCTGCGCGCTCGGGCTGGCAAGATTGCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt                                                          <  2:608172/148‑1 (MQ=25)
                                                                                                tttGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCg                                                   >  1:158737/1‑151 (MQ=255)
                                                                                                    gttcttTTTTGCGCCCGGCGCGCTCGGGCTGGAAAGATTCCCCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  2:579209/149‑1 (MQ=25)
                                                                                                     aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                >  2:451890/1‑150 (MQ=255)
                                                                                                        ttttttGCTCCCTGCGCGCTCCGGCTGGAAAGATTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  2:30282/145‑1 (MQ=37)
                                                                                                        tctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGc                                            >  1:923223/1‑151 (MQ=255)
                                                                                                         ctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                           >  1:244336/1‑151 (MQ=255)
                                                                                                          ttttGCTGCCGCCGCGCTGCGGCGGACAAGATGGCGCCTTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                <  1:738255/145‑1 (MQ=25)
                                                                                                           tttGCTGCCTGCGCGCTACGGCTGGCAGGATTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                           <  1:808142/149‑1 (MQ=255)
                                                                                                                  cctgcGCGCTCCGGCTGGCAAGTTTGCGCCCTTCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATg                                    <  2:951642/149‑1 (MQ=255)
                                                                                                                  cctgcGCGCTCCGGCGGCAAAATTTGCCCCCTCCCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGt                                   <  1:832537/150‑1 (MQ=25)
                                                                                                                    tgcGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGc                                >  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)
                                                                                                                                                     |                                                                                                                                                    
GGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTGATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTTCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCT  >  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: