Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,383 A→G 23.5% intergenic (+189/+40) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3830AG23.5% 137.7 / 40.1 68intergenic (+189/+40)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (9/7);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACT  >  NC_000913/3423249‑3423531
                                                                                                                                      |                                                                                                                                                    
aatttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCtt                                                                                                                                        <  2:399315‑M1/22‑1 (MQ=255)
 atttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTc                                                                                                                                       <  4:457692‑M1/23‑1 (MQ=255)
  tttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTcc                                                                                                                                      <  1:697432‑M1/24‑1 (MQ=255)
       ctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTg                                                                                                                                 <  3:359411‑M1/29‑1 (MQ=255)
         ggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAg                                                                                                                               >  3:457692‑M1/119‑149 (MQ=255)
          gcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                              <  4:641411‑M1/32‑1 (MQ=255)
            ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            >  4:359411‑M1/116‑149 (MQ=255)
            ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt                                                                                                                            >  3:1331794‑M1/116‑149 (MQ=255)
            ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            >  1:399315‑M1/116‑149 (MQ=255)
             gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                                                                                                           >  1:456454‑M1/115‑149 (MQ=255)
             gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                           >  2:697432‑M1/115‑149 (MQ=255)
                 tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt                                                                                                                       >  3:641411‑M1/111‑149 (MQ=255)
                 tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                              >  3:1378499‑M1/111‑142 (MQ=255)
                 tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                              <  4:1378499‑M1/32‑1 (MQ=255)
                            tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         <  2:752766‑M1/37‑1 (MQ=255)
                            tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         >  1:752766‑M1/100‑136 (MQ=255)
                                acctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGc                                                                                                                                                <  4:123090‑M1/14‑1 (MQ=255)
                                acctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGc                                                                                                                                                >  3:123090‑M1/96‑109 (MQ=255)
                                  ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            <  3:327972‑M1/34‑1 (MQ=255)
                                  ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            >  4:327972‑M1/94‑127 (MQ=255)
                                       cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            >  4:1428409‑M1/89‑122 (MQ=255)
                                       cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            <  3:1428409‑M1/34‑1 (MQ=255)
                                        ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                              >  4:429407‑M1/88‑119 (MQ=255)
                                        ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                              <  3:429407‑M1/32‑1 (MQ=255)
                                         catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            >  1:234009‑M1/87‑120 (MQ=255)
                                         catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                            <  2:234009‑M1/34‑1 (MQ=255)
                                            cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                       <  3:1029105/106‑1 (MQ=255)
                                            cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                       >  4:1029105/1‑106 (MQ=40)
                                               gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                             >  4:168154‑M1/81‑113 (MQ=255)
                                               gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                             <  3:168154‑M1/33‑1 (MQ=255)
                                                   tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc                                                                                                                              <  4:178072‑M1/32‑1 (MQ=255)
                                                   tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc                                                                                                                              >  3:178072‑M1/77‑108 (MQ=255)
                                                     agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         <  2:747899‑M1/37‑1 (MQ=255)
                                                     agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         >  1:747899‑M1/75‑111 (MQ=255)
                                                      ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                      <  2:1431621/113‑1 (MQ=255)
                                                      ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                      >  1:1431621/1‑113 (MQ=255)
                                                          ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                       <  3:977342/92‑1 (MQ=255)
                                                          ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                       >  4:977342/1‑92 (MQ=255)
                                                            aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                              >  2:1116491‑M1/68‑99 (MQ=255)
                                                            aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                              <  1:1116491‑M1/32‑1 (MQ=255)
                                                              acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                           >  3:1410685‑M1/66‑100 (MQ=255)
                                                              acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                           <  4:1410685‑M1/35‑1 (MQ=255)
                                                                  gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                          >  4:791340/1‑97 (MQ=255)
                                                                  gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                          <  3:791340/97‑1 (MQ=255)
                                                                      gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                             >  1:865828‑M1/58‑90 (MQ=255)
                                                                      gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                             <  2:865828‑M1/33‑1 (MQ=255)
                                                                       aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   <  3:1036186/83‑1 (MQ=255)
                                                                       aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   >  4:1036186/1‑83 (MQ=255)
                                                                        aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   >  4:260999/1‑82 (MQ=255)
                                                                        aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   <  3:260999/82‑1 (MQ=255)
                                                                                  ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                  >  2:992421/1‑89 (MQ=255)
                                                                                  ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                  <  1:992421/89‑1 (MQ=255)
                                                                                  ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                <  2:253651/107‑1 (MQ=255)
                                                                                  ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                >  1:253651/1‑107 (MQ=255)
                                                                                   gcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                                     >  1:1442283/1‑149 (MQ=255)
                                                                                              cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                           <  1:157606‑M1/35‑1 (MQ=255)
                                                                                              cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                           >  2:157606‑M1/34‑68 (MQ=255)
                                                                                                  tCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                      <  2:1442283/149‑1 (MQ=255)
                                                                                                      agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         <  2:1372437‑M1/37‑1 (MQ=255)
                                                                                                      agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                         >  1:1372437‑M1/26‑62 (MQ=255)
                                                                                                              tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   >  2:734708/1‑44 (MQ=40)
                                                                                                              tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                   <  1:734708/44‑1 (MQ=40)
                                                                                                                        acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                     >  4:1346152/1‑144 (MQ=40)
                                                                                                                        acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                     <  3:1346152/144‑1 (MQ=40)
                                                                                                                                tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                 <  2:539010/76‑1 (MQ=255)
                                                                                                                                tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                 >  1:539010/1‑76 (MQ=255)
                                                                                                                                 aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                           <  4:1303368/129‑1 (MQ=255)
                                                                                                                                 aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                           >  3:1303368/1‑129 (MQ=255)
                                                                                                                                      aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt  >  1:239662/1‑149 (MQ=38)
                                                                                                                                      |                                                                                                                                                    
GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACT  >  NC_000913/3423249‑3423531

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