Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,161,946 A→G A161A (GCA→GCG)  iceT → putative citrate/iron‑citrate/zinc‑citrate efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,161,9460AG96.6% 200.9 / ‑1.8 59A161A (GCA→GCG) iceTputative citrate/iron‑citrate/zinc‑citrate efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/35);  total (23/36)
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.18e-01

GAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGG  >  NC_000913/2161851‑2162042
                                                                                               |                                                                                                
gAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGcgtcgt                                                                                              >  1:913125/1‑100 (MQ=255)
gAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGcgtcgt                                                                                              >  2:736565/1‑100 (MQ=255)
   cAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGc                                                                                           <  2:1019011/100‑1 (MQ=255)
    aaTATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCa                                                                                          >  2:443828/1‑100 (MQ=255)
    aaTATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCa                                                                                          <  1:337832/100‑1 (MQ=255)
     atatatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCAc                                                                                         >  1:767097/1‑100 (MQ=255)
       atatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTg                                                                                       <  1:1905453/100‑1 (MQ=255)
        tatGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTgg                                                                                      >  1:1713429/1‑100 (MQ=255)
         atGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGa                                                                                     >  2:1292025/1‑100 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCGGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCtttt                                                                               <  1:1905818/99‑1 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCttttt                                                                              <  2:255985/100‑1 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCttttt                                                                              <  1:1817136/100‑1 (MQ=255)
                cgATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCtttt                                                                               <  1:976251/99‑1 (MQ=255)
                cgATGACCTTTGCCACGTGGCCCGGTCAGGTCGGTCGGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCt                                                                                  <  1:1223116/96‑1 (MQ=255)
                 gATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTg                                                                             >  2:933007/1‑100 (MQ=255)
                    gACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATc                                                                          <  1:176123/100‑1 (MQ=255)
                    gACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATc                                                                          <  2:1210345/100‑1 (MQ=255)
                    gACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATc                                                                          <  2:1142126/100‑1 (MQ=255)
                           gTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTc                                                                   >  2:1915152/1‑100 (MQ=255)
                            tCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTcc                                                                  >  2:1250268/1‑100 (MQ=255)
                            tCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTcc                                                                  <  1:1853319/100‑1 (MQ=255)
                             cACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCg                                                                 >  1:876153/1‑100 (MQ=255)
                              aCGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCgg                                                                <  1:93770/100‑1 (MQ=255)
                               cGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGt                                                               >  1:1392348/1‑100 (MQ=255)
                                   aCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTgggg                                                           >  1:1551011/1‑100 (MQ=255)
                                   aCCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTgggg                                                           <  2:447431/100‑1 (MQ=255)
                                    cccGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGa                                                          >  1:432504/1‑100 (MQ=255)
                                     ccGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGAt                                                         >  1:1174806/1‑100 (MQ=255)
                                            ggtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGt                                                  >  2:1435197/1‑100 (MQ=255)
                                             gtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTg                                                 <  1:1015469/100‑1 (MQ=255)
                                             gtcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTg                                                 <  1:443826/100‑1 (MQ=255)
                                              tcggtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGc                                                >  2:926183/1‑100 (MQ=255)
                                                 gtcCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGc                                                <  1:644149/97‑1 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCg                                           <  1:1012895/100‑1 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCg                                           <  2:860648/100‑1 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCg                                           >  2:1702819/1‑100 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCg                                           <  1:1810506/100‑1 (MQ=255)
                                                   ccGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTATTTGATCAACATTCCGGTGGGGATCATCGGTGCGATCg                                           <  1:1199378/100‑1 (MQ=255)
                                                    cGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGc                                          <  2:1296122/100‑1 (MQ=255)
                                                     gctgctCGGTCCGGCGCTCGGCGGGCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGcc                                         <  1:795973/100‑1 (MQ=255)
                                                        gctCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCcaca                                      <  1:1292018/100‑1 (MQ=255)
                                                             gTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCt                                 <  2:17689/100‑1 (MQ=255)
                                                              tCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTg                                <  2:252336/100‑1 (MQ=255)
                                                               ccGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGt                               <  2:890808/100‑1 (MQ=255)
                                                               ccGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGt                               <  2:1186183/100‑1 (MQ=255)
                                                                cGGCGCTCGGCGGTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGtt                              <  1:201149/100‑1 (MQ=255)
                                                                       cgatcttCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCg                       <  2:848484/94‑1 (MQ=255)
                                                                        ggcggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATTCCGa                      <  1:1620999/100‑1 (MQ=255)
                                                                         gcggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGaa                     <  2:913126/100‑1 (MQ=255)
                                                                          cggTCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTTATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAAc                    <  2:685483/100‑1 (MQ=255)
                                                                            ccgatctgGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTa                  <  1:581110/95‑1 (MQ=255)
                                                                             tCTGCTGGTGGAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGagatcg                 >  1:848483/1‑95 (MQ=255)
                                                                                 ctgGTGGAGTACGCGGCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTaca                <  1:695254/97‑1 (MQ=255)
                                                                                  tggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCAt            >  2:1901673/1‑100 (MQ=255)
                                                                                   ggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTaagatcgg           >  2:581112/1‑93 (MQ=255)
                                                                                   ggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCAt            >  1:1555403/1‑99 (MQ=255)
                                                                                   ggtggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATg           >  1:720014/1‑100 (MQ=255)
                                                                                     tggAGTACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCa         <  1:1915139/100‑1 (MQ=255)
                                                                                         gTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAcg     >  1:1458459/1‑100 (MQ=255)
                                                                                          tACGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAcgc    >  2:1726172/1‑100 (MQ=255)
                                                                                           aCGCGTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAcgcg   >  1:1791187/1‑100 (MQ=255)
                                                                                            cgcgTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGAcgcg   <  1:1184384/99‑1 (MQ=255)
                                                                                            cgcgTCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACgcgg  >  2:645602/1‑100 (MQ=255)
                                                                                               |                                                                                                
GAGCAATATATGGCGGCGATGACCTTTGTCACGTTACCCGGTCAGGTCGGTCCGCTGCTCGGTCCGGCGCTCGGCGGTCTGCTGGTGGAGTACGCATCGTGGCACTGGATCTTTTTGATCAACATTCCGGTGGGGATTATCGGTGCGATCGCCACATTGCTGTTAATGCCGAACTACACCATGCAGACGCGG  >  NC_000913/2161851‑2162042

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑