Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,001,2971.G32.1% 66.7 / 74.4 53coding (953/2259 nt)xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (19/17);  new base G (8/9);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTC  >  NC_000913/3001162‑3001432
                                                                                                                                        |                                                                                                                                       
cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggg                                                                                                                                      >  1:277344/1‑139 (MQ=255)
 aaCCCTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTGACGGGCACTCCATCGCTTCTGCgggg                                                                                                                                     <  2:265429/139‑1 (MQ=255)
           cGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAg                                                                                                                           >  1:1148670/1‑139 (MQ=255)
             ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCg                                                                                                                                        >  1:712140/1‑124 (MQ=255)
             ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCg                                                                                                                                        <  2:712140/124‑1 (MQ=255)
              cACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg                                                                                                                        >  2:603944/1‑139 (MQ=255)
                             gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCt                                                                                                          >  2:1296593/1‑139 (MQ=255)
                                gggCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                       <  1:1296593/139‑1 (MQ=255)
                                     attagggacgAACCGCGACCGAACACTGCCAGCTTATAGTCTGGATGTTCTGTCTAACACCGGCGGTTATGCATCTCACGGGCACTCCAGCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGca                                                                                                 <  2:1430962/130‑2 (MQ=255)
                                      aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                                 >  2:861377/1‑139 (MQ=255)
                                      aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                                 <  1:467725/139‑1 (MQ=255)
                                             tGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                         <  1:359208/139‑1 (MQ=255)
                                                            cATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                      <  2:613365/96‑1 (MQ=255)
                                                            cATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                      >  1:613365/1‑96 (MQ=255)
                                                               tGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGTGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCta                                                                       <  1:798625/139‑1 (MQ=255)
                                                                       tataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                <  2:277344/139‑1 (MQ=255)
                                                                        ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                               >  2:994760/1‑139 (MQ=255)
                                                                        ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCt                                                                                                                      <  2:856216/83‑1 (MQ=255)
                                                                        ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCt                                                                                                                      >  1:856216/1‑83 (MQ=255)
                                                                             cTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                               >  2:759259/1‑117 (MQ=255)
                                                                             cTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                               <  1:759259/117‑1 (MQ=255)
                                                                               ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCg                                                       >  1:66614/1‑139 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCATTCCATCGCTTCTGCg                                                                                                                                        >  1:1259333/1‑57 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCATTCCATCGCTTCTGCg                                                                                                                                        <  2:1259333/57‑1 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                                  >  2:1129054/1‑111 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                                  <  1:1129054/111‑1 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCt                                                                        >  2:1244863/1‑121 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCt                                                                        <  1:1244863/121‑1 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                      >  2:1013757/1‑139 (MQ=255)
                                                                                gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                      >  1:1114215/1‑139 (MQ=255)
                                                                                      cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTg                                                <  1:1013757/139‑1 (MQ=255)
                                                                                        gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                              >  1:1275702/1‑107 (MQ=255)
                                                                                        gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                              <  2:1275702/107‑1 (MQ=255)
                                                                                              cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCt                                                                                                                      >  1:1306185/1‑61 (MQ=255)
                                                                                              cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCt                                                                                                                      <  2:1306185/61‑1 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTGCCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                           <  1:1191339/132‑1 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                           >  2:1191339/1‑132 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                                    >  2:468917/1‑139 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                                    >  2:236041/1‑139 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggg                                                                                                                                     >  1:558568/1‑41 (MQ=255)
                                                                                                   gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggg                                                                                                                                     <  2:558568/41‑1 (MQ=255)
                                                                                                    cgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTa                                   <  1:994760/139‑1 (MQ=255)
                                                                                                        tATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGgc                               <  1:1184006/139‑1 (MQ=255)
                                                                                                              tctcACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                             >  1:435038/1‑102 (MQ=255)
                                                                                                              tctcACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                             <  2:435038/102‑1 (MQ=255)
                                                                                                                 cACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                >  1:578509/1‑96 (MQ=255)
                                                                                                                 cACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                <  2:578509/96‑1 (MQ=255)
                                                                                                                   cGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGt                    >  1:24826/1‑139 (MQ=255)
                                                                                                                         cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAtt              <  2:975682/139‑1 (MQ=255)
                                                                                                                         cTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                     >  1:188132/1‑99 (MQ=255)
                                                                                                                         cTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                     <  2:188132/99‑1 (MQ=255)
                                                                                                                          tCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTg            >  1:35929/1‑139 (MQ=255)
                                                                                                                          tCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTg            >  1:37438/1‑139 (MQ=255)
                                                                                                                               cGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt       <  2:35929/139‑1 (MQ=255)
                                                                                                                                gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt       <  1:236041/139‑1 (MQ=255)
                                                                                                                                   tCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGt   <  2:37438/139‑1 (MQ=255)
                                                                                                                                    cTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                              >  2:1358657/1‑95 (MQ=255)
                                                                                                                                    cTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                              <  1:1358657/95‑1 (MQ=255)
                                                                                                                                    cTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  >  1:1285507/1‑139 (MQ=255)
                                                                                                                                    cTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  <  2:990376/139‑1 (MQ=255)
                                                                                                                                    cTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  >  1:990376/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTC  >  NC_000913/3001162‑3001432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: