Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,172,259 G→A intergenic (‑48/‑271) ygiZ ← / → mdaB inner membrane protein/NADPH quinone reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,172,2590GA100.0% 129.6 / NA 38intergenic (‑48/‑271)ygiZ/mdaBinner membrane protein/NADPH quinone reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/19);  total (19/19)

ACCAATAAGTAAGCATTATATATAGCAATGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACGTTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTAT  >  NC_000913/3172148‑3172369
                                                                                                               |                                                                                                              
aCCAATAAGTAAGCATTATATATAGCAATGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCa                                                                                     <  2:744608/139‑1 (MQ=255)
          aagCATTATATATAGCAATGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATAtt                                                                           >  2:741678/1‑139 (MQ=255)
                         cAATGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCTTTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCaaaaa                                                            >  2:742222/1‑139 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCTATTTATAATTATCAACATGACACATTGATACTTTCATCt                                                                                                 >  1:3098/1‑99 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCt                                                                                                 <  2:3098/99‑1 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTAc                                                                                            <  1:611039/104‑1 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTAc                                                                                            >  2:611039/1‑104 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAAc                                                         >  1:207695/1‑139 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAAc                                                         <  1:332204/139‑1 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAAc                                                         >  2:362824/1‑139 (MQ=255)
                            tGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAAc                                                         >  1:720633/1‑139 (MQ=255)
                              cTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCg                                                       >  1:43505/1‑139 (MQ=255)
                                   aaTATTGTTTTAATTTTCTTCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCa                                                  <  2:43505/139‑1 (MQ=255)
                                    atatTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCTTTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCaa                                                 <  1:742222/139‑1 (MQ=255)
                                             tAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTaaa                                        <  2:207695/139‑1 (MQ=255)
                                                        tttAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTg                             <  2:145790/139‑1 (MQ=255)
                                                                      gTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACgt                                             <  2:197186/109‑1 (MQ=255)
                                                                      gTACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACgt                                             >  1:197186/1‑109 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGt                              >  1:504652/1‑123 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGt                              <  2:504652/123‑1 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATg              >  2:256465/1‑139 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATg              >  2:56016/1‑139 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATg              >  1:171623/1‑139 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATg              >  1:772496/1‑139 (MQ=255)
                                                                       tACGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATg              >  1:187298/1‑139 (MQ=255)
                                                                        aCGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGt                              <  2:345763/122‑1 (MQ=255)
                                                                        aCGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGt                              >  1:345763/1‑122 (MQ=255)
                                                                         cGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAg            <  1:633821/139‑1 (MQ=255)
                                                                         cGCTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAg            >  2:585159/1‑139 (MQ=255)
                                                                           cTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGa          <  1:56016/139‑1 (MQ=255)
                                                                           cTCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGa          >  1:211926/1‑139 (MQ=255)
                                                                            tCATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGaa         <  2:171623/139‑1 (MQ=255)
                                                                             cATAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAAt        <  1:741678/139‑1 (MQ=255)
                                                                              aTAAAAATCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATg       <  2:211926/139‑1 (MQ=255)
                                                                                 aaaaTCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCt    <  2:187298/139‑1 (MQ=255)
                                                                                   aaTCCATTTATAATTATCAACATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTAt  <  2:772496/139‑1 (MQ=255)
                                                                                                     aaCATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGGATTATGAGGAATGGCTAt  <  2:434861/121‑1 (MQ=255)
                                                                                                     aaCATGACACATTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTAt  >  1:434861/1‑121 (MQ=255)
                                                                                                               |                                                                                                              
ACCAATAAGTAAGCATTATATATAGCAATGCTTCAAATATTGTTTTAATTTTCTGCTTTAACATATATCAGTACGCTCATAAAAATCCATTTATAATTATCAACATGACACGTTGATACTTTCATCTTTTACTTATTCAAAGCCATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTAT  >  NC_000913/3172148‑3172369

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