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

ATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑AAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGG  >  NC_000913/1529975‑1530158
                                                                                                 |                                                                                          
aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑aaa‑aa                                                                                      >  2:909908/1‑100 (MQ=255)
 tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑aaa‑aaa                                                                                     <  2:604559/100‑1 (MQ=255)
   tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑aaa‑aaaaa                                                                                   >  2:129766/1‑99 (MQ=255)
       gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑aaa‑aaaaatatt                                                                               >  1:107037/1‑95 (MQ=255)
       gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑aaa‑aaaaatatt                                                                               >  1:1134913/1‑95 (MQ=255)
         ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAAtattat                                                                              >  2:489147/1‑100 (MQ=255)
          tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATattatt                                                                             >  1:261837/1‑100 (MQ=255)
           aCCTTGGGTTTGATAACAATTCAAATTCTGACATCTTATTTGATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATAttattt                                                                            >  2:711178/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATAttattta                                                                           >  2:27231/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATAttattta                                                                           >  2:348055/1‑100 (MQ=255)
            ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATAttattta                                                                           >  1:498075/1‑100 (MQ=255)
             cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATAttatttat                                                                          <  2:828314/100‑1 (MQ=255)
              ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTtata                                                                         >  1:1187119/1‑100 (MQ=255)
                     tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAaataat                                                                  >  1:898340/1‑100 (MQ=255)
                      gATAAAAAATCAAATTCTGACATCATATATTAAAAATAAAAAACACAAAAGAGCGTTGCGATAAAAAAAAAAA‑A‑AAA‑AAAATATTATTTATAAAAtaatt                                                                 <  2:713062/100‑1 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaatta                                                                >  2:442098/1‑100 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaatta                                                                <  2:656723/100‑1 (MQ=255)
                       aTAACAATTCAAATTCTGACATCGTATTTTACGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaatta                                                                <  1:434939/100‑1 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaattaa                                                               >  1:529897/1‑100 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaattaa                                                               >  1:1098519/1‑100 (MQ=255)
                        tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAAtaattaa                                                               >  2:1034654/1‑100 (MQ=255)
                          aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATg                                                             <  2:444648/100‑1 (MQ=255)
                             aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGAcc                                                          <  1:754784/100‑1 (MQ=255)
                              ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCa                                                         >  2:356203/1‑100 (MQ=255)
                                     tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                  <  1:361656/100‑1 (MQ=255)
                                      cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                  <  2:975056/99‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:742044/100‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:728424/100‑1 (MQ=255)
                                       tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                <  1:1117828/100‑1 (MQ=255)
                                         acataaaaatttATGGACATAAAAAACCAAAGAGCGTTGAAAAATAACTTTA‑A‑AAAAA‑AAAATATTATTTATAAAATAATAAAAAAAAAAAAAAtt                                                  <  2:1168042/87‑1 (MQ=255)
                                         aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                              <  1:129768/100‑1 (MQ=255)
                                           aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                            >  1:232559/1‑100 (MQ=255)
                                           aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                            >  2:240842/1‑100 (MQ=255)
                                           aTCGTATTTTATAGACATAAAACACCAAAGAGCATAGAAATATATCTTTC‑TAA‑AAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                   <  2:18277/93‑1 (MQ=255)
                                            tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                           >  2:686740/1‑100 (MQ=255)
                                              gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                         >  2:645124/1‑100 (MQ=255)
                                                aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑AAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                      <  1:766796/100‑1 (MQ=255)
                                                  tttAAGAACATAAAACACCAAAGAGCATTGAGATATATCTTTA‑T‑GAAAA‑AAAATATTATTTATAAACTAATTAATGCCCCCCAAACTAATAGAGGGCATg                                     <  2:261834/100‑1 (MQ=255)
                                                    tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                   <  2:529893/100‑1 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               >  1:1261555/1‑100 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               <  1:1107179/100‑1 (MQ=255)
                                                        gACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                               >  2:599040/1‑100 (MQ=255)
                                                         aCATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                              >  1:752161/1‑100 (MQ=255)
                                                         aCATAAAACACCAAAGAGCGTTGAGATATAACATTA‑A‑AAAAA‑AAAATATTATTTATAAAATAAATAATAAAAAAAAAATTAGTAGAGGGAATGGTCa                                 <  2:31517/97‑1 (MQ=255)
                                                          cATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                             <  1:240844/100‑1 (MQ=255)
                                                           aTAAAACACCAAAGAGCGTTGAGATATAAAAAAA‑A‑AAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                            <  2:664382/100‑1 (MQ=255)
                                                            tAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                           >  1:535412/1‑100 (MQ=255)
                                                             aaaaCACCAAAGAGCGTTGAGACATAAATTTC‑T‑G‑AAAGAAAATATTACTTATAAAATAATTAATAAAAAAAAAATTAGTAGAGGGAATGGTCATTTtata                          <  2:879383/100‑1 (MQ=255)
                                                              aaaCACCAAAGAGCGTTGAGATATATCCTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                         <  2:370619/100‑1 (MQ=255)
                                                               aaCACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  1:450002/100‑1 (MQ=255)
                                                               aaCACCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  1:57091/100‑1 (MQ=255)
                                                               aaCACCAAAAAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                        <  2:870249/100‑1 (MQ=255)
                                                                acacCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATGTTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                       <  2:675678/100‑1 (MQ=255)
                                                                acacCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                       >  1:597793/1‑100 (MQ=255)
                                                                  acCAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                     <  2:129306/100‑1 (MQ=255)
                                                                    cAAAGAGCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                   >  2:375150/1‑100 (MQ=255)
                                                                         agCGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt              >  1:838457/1‑100 (MQ=255)
                                                                           cGTTGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa            <  1:405551/100‑1 (MQ=255)
                                                                             ttGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAACATAATTCATCACCAAAAAATTAGTTTTGGGAAAGTTAATTATATAAATATCATGCTtaata          >  2:643172/1‑100 (MQ=255)
                                                                             ttGAGATATATCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata          >  2:739815/1‑100 (MQ=255)
                                                                             ttGAGATATATCTTTC‑T‑G‑AAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa         >  1:974528/1‑100 (MQ=255)
                                                                                   tatatCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa    <  1:1049957/100‑1 (MQ=255)
                                                                                    atatCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg   <  1:686744/100‑1 (MQ=255)
                                                                                     tatCTTTC‑T‑GAAAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAgg  >  1:527729/1‑100 (MQ=255)
                                                                                                 |                                                                                          
ATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTC‑T‑G‑AAA‑AAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGG  >  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: ‑