Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 344,704 A→T 6.4% K177* (AAA→TAA) ‡ yahL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,7040AT 6.4% 224.4 / 5.2 62K177* (AAA→TAA) ‡yahLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCA  >  NC_000913/344557‑344844
                                                                                                                                                   |                                                                                                                                            
gatgatGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaa                                                                                                                                             >  4:692528/1‑149 (MQ=255)
     tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                                                        <  4:327032/149‑1 (MQ=255)
              tGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                               <  3:257197/149‑1 (MQ=255)
              tGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                               <  1:1400496/149‑1 (MQ=255)
              tGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                               <  2:875629/149‑1 (MQ=255)
                   ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAagag                                                                                                                                   <  2:177655/140‑1 (MQ=255)
                   ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAagag                                                                                                                                   >  1:177655/1‑140 (MQ=255)
                     aaaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                               >  2:682705/1‑142 (MQ=255)
                     aaaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTa                                                                                                                               <  1:682705/142‑1 (MQ=255)
                               tGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTgg                                                                                                                      <  1:13980/141‑1 (MQ=255)
                               tGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTgg                                                                                                                      >  2:13980/1‑141 (MQ=255)
                               tGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATgg                                                                                                              <  3:1446952/149‑1 (MQ=255)
                               tGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATgg                                                                                                              <  3:692528/149‑1 (MQ=255)
                                         aCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAg                                                                                                                                           >  1:857836/1‑110 (MQ=255)
                                         aCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAg                                                                                                                                           <  2:857836/110‑1 (MQ=255)
                                           gCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATg                                                                                                  >  4:1096882/1‑149 (MQ=255)
                                           gCTTTACCGGGCTTAAGAAATTAGTACAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATg                                                                                                  >  4:779274/1‑149 (MQ=255)
                                                     gCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                            <  3:1042861/97‑1 (MQ=255)
                                                     gCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                            >  4:1042861/1‑97 (MQ=255)
                                                       ttAAGAAATTAGTACAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                                                      <  3:779274/149‑1 (MQ=255)
                                                         aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtc                                                                                                >  2:1228257/1‑137 (MQ=255)
                                                         aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtc                                                                                                <  1:1228257/137‑1 (MQ=255)
                                                           gAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                      <  1:921199/97‑1 (MQ=255)
                                                           gAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCACTAAAGGGaaa                                                                                                                                      >  2:921199/1‑97 (MQ=255)
                                                                    gCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTcc                                                                                       <  1:326532/135‑1 (MQ=255)
                                                                    gCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTcc                                                                                       >  2:326532/1‑135 (MQ=255)
                                                                       gTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                      <  4:392452/85‑1 (MQ=255)
                                                                       gTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                      >  3:392452/1‑85 (MQ=255)
                                                                        tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtg                                                                                                                         <  3:322651/97‑1 (MQ=255)
                                                                        tttATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtg                                                                                                                         >  4:322651/1‑97 (MQ=255)
                                                                           aTCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                                                                  <  3:1096882/149‑1 (MQ=255)
                                                                                ttCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaaa                                                                                                    >  2:722878/1‑110 (MQ=255)
                                                                                ttCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaaa                                                                                                    <  1:722878/110‑1 (MQ=255)
                                                                                         aaCAGATTGCAAATGTATTCTTCTTGAATACCATACGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgc                                                      <  1:973460/147‑1 (MQ=255)
                                                                                         aaCAGATTGCAAATGTATTCTTCTTGAATACCATACGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTgc                                                      >  2:973460/1‑147 (MQ=255)
                                                                                          aCAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACgcgc                                                                         >  1:1009549/1‑127 (MQ=255)
                                                                                          aCAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACgcgc                                                                         <  2:1009549/127‑1 (MQ=255)
                                                                                                     aTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattat                                        >  2:1364145/1‑149 (MQ=255)
                                                                                                      tGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTAt                                                                                                                              <  2:1196010/62‑1 (MQ=255)
                                                                                                      tGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattattatt                                       <  1:1364145/149‑1 (MQ=255)
                                                                                                      tGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAAGTTTCCCTAAAGGGAAAGAGTGTAt                                                                                                                              >  1:1196010/1‑62 (MQ=255)
                                                                                                         aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                                                     >  2:182897/1‑116 (MQ=255)
                                                                                                         aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                                                     <  1:182897/116‑1 (MQ=255)
                                                                                                              tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGcc                                                                        <  3:785388/108‑1 (MQ=255)
                                                                                                              tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGcc                                                                        >  4:785388/1‑108 (MQ=255)
                                                                                                              tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                                                                  <  1:466342/114‑1 (MQ=255)
                                                                                                              tcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                                                                  >  2:466342/1‑114 (MQ=255)
                                                                                                                ttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgt                                                                                                                                  >  1:993702/1‑48 (MQ=255)
                                                                                                                ttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgt                                                                                                                                  <  2:993702/48‑1 (MQ=255)
                                                                                                                     tACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCgg                                                                                     <  3:26375/88‑1 (MQ=255)
                                                                                                                     tACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCgg                                                                                     >  4:26375/1‑88 (MQ=255)
                                                                                                                            tGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCt                 >  4:1284667/1‑149 (MQ=255)
                                                                                                                                tttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTaaaaaa                           >  4:1159354/1‑135 (MQ=255)
                                                                                                                                tttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGCGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGATACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTaaaaaa                           <  3:1159354/135‑1 (MQ=255)
                                                                                                                                      gagaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                    >  4:634346/1‑72 (MQ=255)
                                                                                                                                      gagaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                    <  3:634346/72‑1 (MQ=255)
                                                                                                                                        gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                                               >  2:1151597/1‑91 (MQ=255)
                                                                                                                                        gaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                                               <  1:1151597/91‑1 (MQ=255)
                                                                                                                                          tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAc                    >  2:721876/3‑132 (MQ=255)
                                                                                                                                          tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAc                    <  1:721876/130‑1 (MQ=255)
                                                                                                                                          tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGtt              >  4:1091517/3‑138 (MQ=255)
                                                                                                                                          tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGtt              <  3:1091517/136‑1 (MQ=255)
                                                                                                                                           ttttCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCa  >  2:981613/1‑149 (MQ=255)
                                                                                                                                                  tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                             >  4:1084933/1‑67 (MQ=255)
                                                                                                                                                  tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                             <  3:1084933/67‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCA  >  NC_000913/344557‑344844

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