Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,530,076 (A)7→8 coding (155/483 nt) ydcD → putative immunity protein for RhsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,0691.A91.5% 192.5 / 10.8 59coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base A (27/27);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGG  >  NC_000913/1529975‑1530158
                                                                                                 |                                                                                          
aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaa                                                                                      >  2:909908/1‑100 (MQ=255)
 tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaa                                                                                     <  2:604559/100‑1 (MQ=255)
   tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaa                                                                                   >  2:129766/1‑99 (MQ=255)
       gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                               >  1:107037/1‑95 (MQ=255)
       gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                               >  1:1134913/1‑95 (MQ=255)
         ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAtattat                                                                              >  2:489147/1‑100 (MQ=255)
          tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                             >  1:261837/1‑100 (MQ=255)
           aCCTTGGGTTTGATAACAATTCAAATTCTGACATCTTATTTGATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                            >  2:711178/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                           >  2:27231/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                           >  2:348055/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                           >  1:498075/1‑100 (MQ=255)
             cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                          <  2:828314/100‑1 (MQ=255)
              ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTtata                                                                         >  1:1187119/1‑100 (MQ=255)
                     tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataat                                                                  >  1:898340/1‑100 (MQ=255)
                      gATAAAAAATCAAATTCTGACATCATATATTAAAAATAAAAAACACAAAAGAGCGTTGCGATAAAAAAAAAAAAAAAAAAATATTATTTATAAAAtaatt                                                                 <  2:713062/100‑1 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                >  2:442098/1‑100 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                <  2:656723/100‑1 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTACGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                <  1:434939/100‑1 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                               >  1:529897/1‑100 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                               >  1:1098519/1‑100 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                               >  2:1034654/1‑100 (MQ=255)
                          aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATg                                                             <  2:444648/100‑1 (MQ=255)
                             aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                          <  1:754784/100‑1 (MQ=255)
                              ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                         >  2:356203/1‑100 (MQ=255)
                                     tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                  <  1:361656/100‑1 (MQ=255)
                                      cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                  <  2:975056/99‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:742044/100‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:728424/100‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:1117828/100‑1 (MQ=255)
                                         acataaaaatttATGGACATAAAAAACCAAAGAGCGTTGAAAAATAACTTTAAAAAAAAAAATATTATTTATAAAATAATAAAAAAAAAAAAAAtt                                                  <  2:1168042/87‑1 (MQ=255)
                                         aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                              <  1:129768/100‑1 (MQ=255)
                                           aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                            >  1:232559/1‑100 (MQ=255)
                                           aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                            >  2:240842/1‑100 (MQ=255)
                                           aTCGTATTTTATAGACATAAAACACCAAAGAGCATAGAAATATATCTTTCTAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                   <  2:18277/93‑1 (MQ=255)
                                            tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                           >  2:686740/1‑100 (MQ=255)
                                              gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                         >  2:645124/1‑100 (MQ=255)
                                                aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                      <  1:766796/100‑1 (MQ=255)
                                                  tttAAGAACATAAAACACCAAAGAGCATTGAGATATATCTTTATGAAAAAAAATATTATTTATAAACTAATTAATGCCCCCCAAACTAATAGAGGGCATg                                     <  2:261834/100‑1 (MQ=255)
                                                    tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                   <  2:529893/100‑1 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               >  1:1261555/1‑100 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               <  1:1107179/100‑1 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               >  2:599040/1‑100 (MQ=255)
                                                         aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                              >  1:752161/1‑100 (MQ=255)
                                                         aCATAAAACACCAAAGAGCGTTGAGATATAACATTAAAAAAAAAAATATTATTTATAAAATAAATAATAAAAAAAAAATTAGTAGAGGGAATGGTCa                                 <  2:31517/97‑1 (MQ=255)
                                                          cATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                             <  1:240844/100‑1 (MQ=255)
                                                           aTAAAACACCAAAGAGCGTTGAGATATAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                            <  2:664382/100‑1 (MQ=255)
                                                            tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                           >  1:535412/1‑100 (MQ=255)
                                                             aaaaCACCAAAGAGCGTTGAGACATAAATTTCTGAAAGAAAATATTACTTATAAAATAATTAATAAAAAAAAAATTAGTAGAGGGAATGGTCATTTtata                          <  2:879383/100‑1 (MQ=255)
                                                              aaaCACCAAAGAGCGTTGAGATATATCCTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                         <  2:370619/100‑1 (MQ=255)
                                                               aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  1:450002/100‑1 (MQ=255)
                                                               aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  1:57091/100‑1 (MQ=255)
                                                               aaCACCAAAAAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  2:870249/100‑1 (MQ=255)
                                                                acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATGTTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                       <  2:675678/100‑1 (MQ=255)
                                                                acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                       >  1:597793/1‑100 (MQ=255)
                                                                  acCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                     <  2:129306/100‑1 (MQ=255)
                                                                    cAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                   >  2:375150/1‑100 (MQ=255)
                                                                         agCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt              >  1:838457/1‑100 (MQ=255)
                                                                           cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa            <  1:405551/100‑1 (MQ=255)
                                                                             ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAACATAATTCATCACCAAAAAATTAGTTTTGGGAAAGTTAATTATATAAATATCATGCTtaata          >  2:643172/1‑100 (MQ=255)
                                                                             ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata          >  2:739815/1‑100 (MQ=255)
                                                                             ttGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa         >  1:974528/1‑100 (MQ=255)
                                                                                   tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa    <  1:1049957/100‑1 (MQ=255)
                                                                                    atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg   <  1:686744/100‑1 (MQ=255)
                                                                                     tatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAgg  >  1:527729/1‑100 (MQ=255)
                                                                                                 |                                                                                          
ATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGG  >  NC_000913/1529975‑1530158

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