Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,527,849 A→G 9.1% I302V (ATC→GTC)  yfeH → putative inorganic ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,527,8490AG 9.1% 127.4 / 4.3 44I302V (ATC→GTC) yfeHputative inorganic ion transporter
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (2/2);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CCGACTGAGCTTCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCA  >  NC_000913/2527701‑2527996
                                                                                                                                                    |                                                                                                                                                   
ccGACTGAGCTTCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCGCCTGAGGATTTTCCATCAGa                                                                                                                                                     <  2:195006/149‑1 (MQ=255)
     tGAGCTTCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCa                                                                                                                                                 <  4:60012/148‑1 (MQ=255)
       aGCTTCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAAtt                                                                                                                                              <  1:168339/149‑1 (MQ=255)
          ttCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGa                                                                                                                                            <  4:77794/148‑1 (MQ=255)
           tCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGAt                                                                                                                                           >  3:255992/1‑148 (MQ=255)
           tCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGAt                                                                                                                                           <  3:170006/148‑1 (MQ=255)
                aaGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCt                                                                                                                                      >  3:261164/1‑148 (MQ=255)
                    cAGATGAAATTACTATCGACTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGc                                                                                                                                  >  4:107019/1‑148 (MQ=255)
                         gAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCt                                                                                                                            <  1:289563/149‑1 (MQ=255)
                           aaTTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTgg                                                                                                                          >  1:68535/1‑149 (MQ=255)
                                 tATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGgc                                                                                                                         >  1:279627/1‑144 (MQ=255)
                                 tATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGgc                                                                                                                         <  2:279627/144‑1 (MQ=255)
                                    cGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGa                                                                                                                  >  4:297403/1‑148 (MQ=255)
                                     gTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCCGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATa                                                                                                                <  2:92353/149‑1 (MQ=255)
                                        ttttGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACaa                                                                                                             >  1:323897/1‑149 (MQ=255)
                                         tttGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTTCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACaa                                                                                                             >  4:129462/1‑148 (MQ=255)
                                         tttGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTTCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACaa                                                                                                             <  3:129462/148‑1 (MQ=255)
                                                     aaaaGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGa                                                                                                 >  4:156077/1‑148 (MQ=255)
                                                           gTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGa                                                                                                                  >  2:311066/1‑125 (MQ=255)
                                                           gTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGa                                                                                                                  <  1:311066/125‑1 (MQ=255)
                                                              tGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACa                                                                                        <  4:261164/148‑1 (MQ=255)
                                                              tGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAg                                                                                       <  2:323897/149‑1 (MQ=255)
                                                                                    aaaCATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGGTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAg                                                                                       >  1:264332/1‑127 (MQ=255)
                                                                                    aaaCATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGGTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAg                                                                                       >  1:207844/1‑127 (MQ=255)
                                                                                    aaaCATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGGTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAg                                                                                       <  2:207844/127‑1 (MQ=255)
                                                                                    aaaCATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGGTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAg                                                                                       <  2:264332/127‑1 (MQ=255)
                                                                                          tCTGTTCCCCACATCGGTGGTCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGCCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATaaa                                                            <  4:166427/148‑1 (MQ=255)
                                                                                          tCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATaaa                                                            >  3:178523/1‑148 (MQ=255)
                                                                                            tGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGc                                                          <  4:178523/148‑1 (MQ=255)
                                                                                                  ccACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGc                                                   >  2:61228/1‑149 (MQ=255)
                                                                                                       tCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGAc                                               >  4:73446/1‑148 (MQ=255)
                                                                                                             aTCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAg                                                           <  1:382238/130‑1 (MQ=255)
                                                                                                             aTCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAg                                                           >  2:382238/1‑130 (MQ=255)
                                                                                                              tCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAgg                                       <  1:143639/149‑1 (MQ=255)
                                                                                                                 gTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGc                                    >  2:271352/1‑149 (MQ=255)
                                                                                                                   atgatgGTGCTGCCCCTGATGATTTTCCATCAGAGCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCt                                   <  4:91782/148‑1 (MQ=255)
                                                                                                                     gatgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTgg                                 <  3:156077/148‑1 (MQ=255)
                                                                                                                     gatgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTgg                                 <  4:255992/148‑1 (MQ=255)
                                                                                                                      atgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGa                                <  3:107019/148‑1 (MQ=255)
                                                                                                                       tgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGAcc                              >  1:377357/1‑149 (MQ=255)
                                                                                                                          tgctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGt                           <  2:377357/149‑1 (MQ=255)
                                                                                                                           gctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGtt                          <  2:68535/149‑1 (MQ=255)
                                                                                                                                         ttttCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTc             >  3:105839/1‑148 (MQ=255)
                                                                                                                                                 cAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTaa     >  4:131123/1‑148 (MQ=255)
                                                                                                                                                   gATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATca  >  2:400395/1‑149 (MQ=255)
                                                                                                                                                    aTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATca  >  4:59415/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGACTGAGCTTCAATAAGGCAGATGAAATTACTATCGTCTTTTGTGGTTCGAAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCA  >  NC_000913/2527701‑2527996

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: ‑