Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,035,593 G→A 6.4% intergenic (+350/‑242) yedS → / → hchA pseudogene, outer membrane protein homology; putative outer membrane protein/glyoxalase III and Hsp31 molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,035,5930GA 6.4% 199.6 / 2.9 62intergenic (+350/‑242)yedS/hchApseudogene, outer membrane protein homology; putative outer membrane protein/glyoxalase III and Hsp31 molecular chaperone
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base A (2/2);  total (31/31)
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.91e-01

TCTTTATGAATTAATCCGTAACAAACACGCTTCCCATATAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCG  >  NC_000913/2035449‑2035721
                                                                                                                                                |                                                                                                                                
tCTTTATGAATTAATCCGTAACAAACACGCTTCCCATATAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCa                                                                                                                               >  2:442386/1‑148 (MQ=255)
                     caaacaCGCTTCCCATATAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGt                                                                                                          >  4:242116/1‑148 (MQ=255)
                                       aGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                                                                                  >  4:501917/1‑122 (MQ=255)
                                       aGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                                                                                  <  3:501917/122‑1 (MQ=255)
                                       aGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                                                                                  <  1:41289/122‑1 (MQ=255)
                                       aGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                                                                                  >  2:41289/1‑122 (MQ=255)
                                              tttGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcga                                                                                 <  3:242116/148‑1 (MQ=255)
                                                      aaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcc                                                                                                                             <  1:1088532/96‑1 (MQ=255)
                                                      aaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcc                                                                                                                             >  2:1088532/1‑96 (MQ=255)
                                                          gTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCa                                                                     <  1:442386/148‑1 (MQ=255)
                                                                     aaGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCAt                                                                                                >  1:1074257/1‑110 (MQ=255)
                                                                     aaGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCAt                                                                                                <  2:1074257/110‑1 (MQ=255)
                                                                       gCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                                                                             <  2:673915/111‑1 (MQ=255)
                                                                       gCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                                                                             >  1:673915/1‑111 (MQ=255)
                                                                        cTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAccc                                                                                                                            >  4:928713/1‑79 (MQ=255)
                                                                        cTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAccc                                                                                                                            <  3:928713/79‑1 (MQ=255)
                                                                                 tACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATaa                                                                 >  4:1047327/1‑129 (MQ=255)
                                                                                 tACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATaa                                                                 <  3:1047327/129‑1 (MQ=255)
                                                                                  aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                                                                              <  3:1046768/115‑1 (MQ=255)
                                                                                  aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                                                                              <  3:809835/115‑1 (MQ=255)
                                                                                  aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                                                                              >  4:1046768/1‑115 (MQ=255)
                                                                                  aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                                                                              >  4:809835/1‑115 (MQ=255)
                                                                                     aaTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcc                                                                                                                             <  4:574714/65‑1 (MQ=255)
                                                                                     aaTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAcc                                                                                                                             >  3:574714/1‑65 (MQ=255)
                                                                                     aaTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                                                                             <  3:821700/97‑1 (MQ=255)
                                                                                     aaTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGACAATGGTCTAAAATCATTGa                                                                                             >  4:821700/1‑97 (MQ=255)
                                                                                      aTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATAtt                                         >  3:1016792/1‑148 (MQ=255)
                                                                                        ccGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCa                                       <  4:1016792/148‑1 (MQ=255)
                                                                                             atatAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                                                                             >  3:1032660/1‑89 (MQ=255)
                                                                                             atatAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                                                                             <  4:1032660/89‑1 (MQ=255)
                                                                                                   atatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACctct                                                                                                                   >  2:185809/1‑61 (MQ=255)
                                                                                                   atatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACctct                                                                                                                   <  1:185809/61‑1 (MQ=255)
                                                                                                    tatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGt                                                                           <  3:660461/100‑1 (MQ=255)
                                                                                                    tatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGt                                                                           >  4:660461/1‑100 (MQ=255)
                                                                                                          tcaGCACCACCATATCCCACAGCACTAAATCTCTCCCCACCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGcatcat                     <  1:242350/148‑1 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTc                                        <  3:438106/129‑1 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTc                                        >  4:438106/1‑129 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTACTAATgg                                                                                                           <  2:612215/62‑1 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTAATAATgg                                                                                                           >  1:612215/1‑62 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCACCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGcatcat                     >  4:614671/1‑148 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCACCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGcatcat                     >  2:242350/1‑148 (MQ=255)
                                                                                                          tcaGCACCACAATATTCCACAGCACTAAATCTCTCCCCACCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGcatcat                     <  3:614671/148‑1 (MQ=255)
                                                                                                           caGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAg                                                                <  2:652811/104‑1 (MQ=255)
                                                                                                           caGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAg                                                                >  1:652811/1‑104 (MQ=255)
                                                                                                                ccacAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                                                                                  <  1:184463/81‑1 (MQ=255)
                                                                                                                ccacAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                                                                                  >  2:184463/1‑81 (MQ=255)
                                                                                                                ccacAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATa                                                                  >  4:335417/1‑97 (MQ=255)
                                                                                                                ccacAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATa                                                                  <  3:335417/97‑1 (MQ=255)
                                                                                                                       atTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCt                                                                       <  3:802918/85‑1 (MQ=255)
                                                                                                                       atTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCt                                                                       >  4:802918/1‑85 (MQ=255)
                                                                                                                             caGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATc                                                                                                  >  1:301647/1‑52 (MQ=255)
                                                                                                                             caGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATc                                                                                                  <  2:301647/52‑1 (MQ=255)
                                                                                                                                   tAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAA‑TTGGCGTGAATAATATTCa                                       >  3:284516/1‑104 (MQ=255)
                                                                                                                                   tAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAA‑TTGGCGTGAATAATATTCa                                       <  4:284516/104‑1 (MQ=255)
                                                                                                                                    aaaTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAacac                                   <  2:837925/108‑1 (MQ=255)
                                                                                                                                    aaaTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAacac                                   >  1:837925/1‑108 (MQ=255)
                                                                                                                                     aaTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTg                              <  1:261583/112‑1 (MQ=255)
                                                                                                                                     aaTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTg                              >  2:261583/1‑112 (MQ=255)
                                                                                                                                      aTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAgccg  <  4:63726/139‑1 (MQ=255)
                                                                                                                                      aTCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAgccg  >  3:63726/1‑139 (MQ=255)
                                                                                                                                            ccccGCCACCCCGTACCTCTGACAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTAc               <  3:256863/120‑1 (MQ=255)
                                                                                                                                            ccccGCCACCCCGTACCTCTGACAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTAc               >  4:256863/1‑120 (MQ=255)
                                                                                                                                                |                                                                                                                                
TCTTTATGAATTAATCCGTAACAAACACGCTTCCCATATAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCG  >  NC_000913/2035449‑2035721

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