Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,806,042 (C)6→7 coding (42/759 nt) ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,0361.C94.2% 174.3 / 2.5 52coding (48/759 nt)ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (25/24);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGG  >  NC_000913/1805905‑1806169
                                                                                                                                    |                                                                                                                                       
cTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                <  1:1082359/139‑2 (MQ=255)
cTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                <  1:135941/139‑2 (MQ=255)
cTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                >  2:1082359/1‑138 (MQ=255)
 ttttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAGACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                <  1:186847/138‑2 (MQ=255)
 ttttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccccccg                                                                                                                               <  1:564048/139‑3 (MQ=255)
 ttttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                >  2:186847/1‑137 (MQ=255)
  tttCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccccccgc                                                                                                                              <  1:668142/139‑4 (MQ=255)
             cTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGAt                                                                                                                   <  2:1216998/139‑1 (MQ=255)
                aTTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATgg                                                                                                                 >  2:639813/1‑139 (MQ=255)
                aTTGTGCAAGATATCCTGCATTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATgg                                                                                                                 >  2:17065/1‑139 (MQ=255)
                    tgCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAgg                                                                                                             >  1:1274486/1‑139 (MQ=255)
                    tgCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAgg                                                                                                             >  1:552773/1‑139 (MQ=255)
                    tgCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAgg                                                                                                             >  1:493315/1‑139 (MQ=255)
                      cAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  1:259885/139‑1 (MQ=255)
                      cAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           >  1:686842/1‑139 (MQ=255)
                      cAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  2:676387/139‑1 (MQ=255)
                              ccTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           >  1:933542/1‑131 (MQ=255)
                              ccTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  2:933542/131‑1 (MQ=255)
                              ccTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTc                                                                                                   >  1:1198914/1‑139 (MQ=255)
                                tGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCaa                                                                                                 <  2:686842/139‑1 (MQ=255)
                                   gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATg                                                                                                                                          <  2:1230832/97‑1 (MQ=255)
                                   gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATg                                                                                                                                          >  1:1230832/1‑97 (MQ=255)
                                   gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAg                                                                                              >  2:91261/1‑139 (MQ=255)
                                    ttCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                                                             >  1:518764/1‑139 (MQ=255)
                                      cAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCtt                                                                                           <  2:1071788/139‑1 (MQ=255)
                                                tCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCaa                                                                                 <  2:552773/139‑1 (MQ=255)
                                                 cGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACg                                                                                                                     >  1:4661/1‑102 (MQ=255)
                                                 cGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACg                                                                                                                     <  2:4661/102‑1 (MQ=255)
                                                  gAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAaga                                                                               <  1:407192/139‑1 (MQ=255)
                                                   aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  2:310841/110‑1 (MQ=255)
                                                   aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           >  1:310841/1‑110 (MQ=255)
                                                   aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAagag                                                                              >  1:960360/1‑139 (MQ=255)
                                                          ttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                                                       <  2:1198914/139‑1 (MQ=255)
                                                          ttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                                                       <  2:365834/139‑1 (MQ=255)
                                                            ttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaa                                                                     <  2:1108281/139‑1 (MQ=255)
                                                               gcgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTAcaacaa                                                                  <  1:91261/139‑1 (MQ=255)
                                                                 gcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACaaaa                                                                >  1:354656/1‑139 (MQ=255)
                                                                      cAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAaga                                                           <  1:550851/139‑1 (MQ=255)
                                                                              tcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcg                                                   <  2:1274486/139‑1 (MQ=255)
                                                                               ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                                                       <  1:1210382/118‑1 (MQ=255)
                                                                               ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                                                       >  2:1210382/1‑118 (MQ=255)
                                                                               ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcgc                                                  >  2:801514/1‑139 (MQ=255)
                                                                                   gACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCtt                                                                                           >  1:322055/1‑94 (MQ=255)
                                                                                   gACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCtt                                                                                           <  2:322055/94‑1 (MQ=255)
                                                                                   gACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTc                                              >  1:546413/1‑139 (MQ=255)
                                                                                     cAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                <  2:402395/54‑2 (MQ=255)
                                                                                     cAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccccccc                                                                                                                                >  1:402395/1‑53 (MQ=255)
                                                                                        tAAAAACGGTATCATCGGCAGCTGCATTCAGTGACGCAAACAAGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                                                             <  1:928322/87‑1 (MQ=255)
                                                                                        tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  2:1088828/73‑1 (MQ=255)
                                                                                        tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           >  1:1088828/1‑73 (MQ=255)
                                                                                        tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                                                             >  2:928322/1‑87 (MQ=255)
                                                                                        tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCa                                                                                         <  1:899876/91‑1 (MQ=255)
                                                                                        tAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCa                                                                                         >  2:899876/1‑91 (MQ=255)
                                                                                                        ggCAGCTGCATTCAGGGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATc                         <  1:724084/139‑1 (MQ=255)
                                                                                                                 aTTCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAg                <  2:354656/139‑1 (MQ=255)
                                                                                                                 aTTCAGTGACGCAAACATGCCGCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAg                <  2:420480/139‑1 (MQ=255)
                                                                                                                  ttCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           <  2:424152/47‑1 (MQ=255)
                                                                                                                  ttCAGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGaa                                                                                                           >  1:424152/1‑47 (MQ=255)
                                                                                                                  ctCAGTGACGCAAACATGCTCCCCCGCCAACATAACGATGGCAGGAACTGGCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGc               <  1:1172978/138‑1 (MQ=255)
                                                                                                                     aGTGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTc            >  1:1053134/1‑139 (MQ=255)
                                                                                                                       tGACGCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTccc          >  1:604433/1‑139 (MQ=255)
                                                                                                                           gCAAACATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGa      <  2:731415/139‑1 (MQ=255)
                                                                                                                               aCATGCCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAgg  >  1:1038042/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGG  >  NC_000913/1805905‑1806169

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: