Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,899,821 G→T 7.8% I151I (ATC→ATA)  ygcW ← putative SDR family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,899,8210GT 7.8% 179.6 / 5.1 51I151I (ATC→ATA) ygcWputative SDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base T (2/2);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ATAATAGCCAGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAGCTTACAAATAC  >  NC_000913/2899677‑2899966
                                                                                                                                                |                                                                                                                                                 
ataataGCCAGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGatat                                                                                                                                               <  2:524722/149‑1 (MQ=255)
ataataGCCAGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGatat                                                                                                                                               >  4:417876/1‑149 (MQ=255)
  aataGCCAGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCCTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATta                                                                                                                                             >  4:84560/1‑149 (MQ=255)
         aGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAAttt                                                                                                                                      <  1:618508/149‑1 (MQ=255)
                 tACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCtt                                                                                                                              >  3:1228/1‑149 (MQ=255)
                      tttACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTg                                                                                                                         <  4:340664/149‑1 (MQ=255)
                            tGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGAt                                                                                                                   >  3:179492/1‑149 (MQ=255)
                            tGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTCAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGAt                                                                                                                   >  4:455832/1‑149 (MQ=255)
                             gAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATc                                                                                                                  <  2:558051/149‑1 (MQ=255)
                                  ttatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAt                                                                                                             >  1:284582/1‑149 (MQ=255)
                                   tatattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAtt                                                                                                            >  2:304716/1‑149 (MQ=255)
                                     tattGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAtttt                                                                                                          >  4:263973/1‑149 (MQ=255)
                                      attGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAAttttt                                                                                                         <  1:574281/149‑1 (MQ=255)
                                            cTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCCTGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGcg                                                                                                   >  2:405467/1‑148 (MQ=255)
                                             tAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACATATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCtt                                                                                                  >  3:304778/1‑149 (MQ=255)
                                                 tcatcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATa                                                                                              <  3:455832/149‑1 (MQ=255)
                                                    tcaCAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCt                                                                                           >  4:143666/1‑149 (MQ=255)
                                                              cTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAAt                                                                                 >  1:2797/1‑149 (MQ=255)
                                                               tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGACCAATGAACAGa                                                                                                                                                  <  3:497410/83‑1 (MQ=255)
                                                               tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGa                                                                                                                                                  >  4:497410/1‑83 (MQ=255)
                                                               tttGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATg                                                                                >  1:289185/1‑149 (MQ=255)
                                                                   gTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATgcggc                                                                            <  2:14844/149‑1 (MQ=255)
                                                                          cGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAgg                                                                     >  2:79512/1‑149 (MQ=255)
                                                                          cGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAgg                                                                     <  3:84560/149‑1 (MQ=255)
                                                                               agagCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTcac                                                                <  3:417876/149‑1 (MQ=255)
                                                                                gagCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTcaca                                                               <  3:143666/149‑1 (MQ=255)
                                                                                  gCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACatc                                                             <  2:34463/149‑1 (MQ=255)
                                                                                       tttAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGGATGCGGCGGTCAGGTTCACatcaatca                                                        >  2:66313/1‑149 (MQ=255)
                                                                                       tttAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACatcaatca                                                        >  2:503789/1‑149 (MQ=255)
                                                                                       tttAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACatcaatca                                                        >  2:422390/1‑149 (MQ=255)
                                                                                        ttAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACATATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAatcat                                                       <  4:304778/149‑1 (MQ=255)
                                                                                                 gAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCa                                              >  4:12504/1‑149 (MQ=255)
                                                                                                                cATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGaa                               <  3:216505/149‑1 (MQ=255)
                                                                                                                cATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGaa                               >  2:495030/1‑149 (MQ=255)
                                                                                                                 aTTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAg                              <  2:284582/149‑1 (MQ=255)
                                                                                                                    gTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTcc                           <  1:495030/149‑1 (MQ=255)
                                                                                                                    gTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTcc                           >  3:154494/1‑149 (MQ=255)
                                                                                                                     tCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGt                                                                         <  1:218137/102‑1 (MQ=255)
                                                                                                                     tCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGt                                                                         >  2:218137/1‑102 (MQ=255)
                                                                                                                      ccACCTAAGTAAGAGAACAATGAACATATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAAt                                                                                 <  4:188263/93‑1 (MQ=255)
                                                                                                                      ccACCTAAGTAAGAGAACAATGAACATATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAAt                                                                                 >  3:188263/1‑93 (MQ=255)
                                                                                                                          cTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATgcggcg                                                                           >  3:428306/1‑95 (MQ=255)
                                                                                                                          cTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATgcggcg                                                                           <  4:428306/95‑1 (MQ=255)
                                                                                                                              gtaagAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTAt                 >  4:95966/1‑149 (MQ=255)
                                                                                                                                aagAGAACAATGAACAGATATTAATAATTGTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTc               <  1:624650/149‑1 (MQ=255)
                                                                                                                                  gagaACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAg             <  4:230199/149‑1 (MQ=255)
                                                                                                                                        aaTGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGa                                                   >  1:623935/1‑105 (MQ=255)
                                                                                                                                        aaTGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGa                                                   <  2:623935/105‑1 (MQ=255)
                                                                                                                                        aaTGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAGCTTACa       >  1:291722/1‑149 (MQ=255)
                                                                                                                                            aaCAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAGCTTACAAATa   >  2:222545/1‑149 (MQ=255)
                                                                                                                                             aCAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAGCTTACAAATAc  >  2:590798/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATAATAGCCAGGGGCGATACCATTTACCTGAATATTATATTGACCTAGTTCATCACAATAAGCTTTGGTGAACCCGGCAAGAGCATGTTTAGTGGCAGAATATGCAGGTGACCATTGTCCACCTAAGTAAGAGAACAATGAACAGATATTAATAATTTTGCCGCTTTTCTGCGGGATCATAATTTTTGCAGCTTCATAGCTTAACTCGAATGCGGCGGTCAGGTTCACATCAATCATCGGATCCCAGTCGGCACGACCGAAGTCCAGCACCTTATTCAGCTTACAAATAC  >  NC_000913/2899677‑2899966

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