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.A95.0% 204.9 / 3.9 60coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (27/30);  total (28/32)
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.57e-01

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                     |                                                                                             
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                <  2:195696/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                <  1:895278/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                >  1:1570005/1‑100 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑a                                                                                              >  1:109254/1‑100 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑a                                                                                              >  2:436167/1‑100 (MQ=255)
      tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaa                                                                                         <  1:670324/100‑1 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaa                                                                                       >  2:551630/1‑99 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaat                                                                                      <  2:1214670/100‑3 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aaaaaaaatat                                                                                    >  1:642053/1‑96 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta                                                                               >  1:1447697/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta                                                                               <  1:599202/100‑1 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATAttattta                                                                               >  2:571954/1‑99 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTtata                                                                             >  2:1266457/1‑100 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAaat                                                                         >  2:1178070/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAAAAAAAA‑AAAAAAAAATATTATTTATAAaataat                                                                      <  2:109254/99‑1 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACAAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaatt                                                                     <  2:657005/100‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:275575/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  1:1393542/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAAtaattaat                                                                  >  2:1454916/1‑100 (MQ=255)
                               aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATg                                                                 >  1:337504/1‑100 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                              <  1:1286379/100‑1 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                              >  2:1176911/1‑99 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                          <  1:1170897/100‑1 (MQ=255)
                                        aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        <  1:948447/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     >  1:277106/1‑100 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     <  1:1034004/100‑1 (MQ=255)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATAAAAAAAAAA‑AAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                  <  1:649740/100‑1 (MQ=255)
                                              aCAACGTATATTATGGAAATAAAAAACAAAAGAGGGTTGAGATATAAATTTCA‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   <  2:448403/99‑1 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                 <  2:1272649/99‑1 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  1:549300/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:642051/100‑1 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  1:1316875/1‑100 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             <  1:976689/100‑1 (MQ=255)
                                                    tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa                                            <  1:436168/100‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          >  1:1539552/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          <  1:381531/100‑1 (MQ=255)
                                                       tttATGGACATAAAACACCAATGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTCTAAAATGATTAATGACCAACAAACTAGTatagtgcat                                          >  2:32462/1‑94 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                       <  1:1266462/100‑1 (MQ=255)
                                                         tATGGACATAAAAAACTAAAGAGCGTTGAAATATAACTATCTCA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:210216/99‑1 (MQ=255)
                                                          aaagaaaaaaaaCACCAAAGAGAGTTAAGATAAAACAATCT‑AAAAAAAAATATTATTTATAAAAAAAAAAAAAAAAAAAAAATTAGTAGAGGGAATGGt                                       <  2:1403027/91‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  1:1178075/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  2:713586/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:1384231/1‑100 (MQ=255)
                                                             gACATAAAACACAAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:606844/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:20668/100‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  1:194820/100‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 >  1:827712/1‑100 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                >  1:449618/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  1:571955/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGTGGTAATGGTCATTTTATaa                             >  2:1354480/1‑100 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:1345082/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:699974/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATTACCAAAAAATTAGTTGAGGGAATGGTCTTTTTATAAg                            >  2:1228566/1‑99 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           >  1:933306/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          >  1:800046/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        <  1:1186149/100‑1 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta                 >  1:427192/1‑100 (MQ=255)
                                                                                   tGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             >  2:344450/1‑100 (MQ=255)
                                                                                   tGAGATATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  2:1569995/98‑1 (MQ=255)
                                                                                   tGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTaataat            <  1:518329/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:186189/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  2:1098650/1‑100 (MQ=255)
                                                                                      gaTATATCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  1:998740/1‑100 (MQ=255)
                                                                                        tatatCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  2:1467584/100‑1 (MQ=255)
                                                                                        tatatCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  2:967604/1‑100 (MQ=255)
                                                                                            tCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    <  2:449617/100‑1 (MQ=255)
                                                                                            tCTTTCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    <  1:617474/100‑1 (MQ=255)
                                                                                              tttCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  2:1384223/100‑1 (MQ=255)
                                                                                               ttCT‑GAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATGAGTAGAGGGAATGGTCATTTTATAAGGATAATGGTTAATAATTACAGGACGc  >  2:1260512/1‑99 (MQ=255)
                                                                                                     |                                                                                             
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162

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