Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,012,766 A→G T256A (ACA→GCA)  pqiA → intermembrane transport protein PqiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,012,7660AG100.0% 93.1 / NA 27T256A (ACA→GCA) pqiAintermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/14);  total (13/14)

AGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAA  >  NC_000913/1012636‑1012897
                                                                                                                                  |                                                                                                                                   
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTcgcg                                                                                                                             >  2:201126/1‑139 (MQ=255)
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTcgcg                                                                                                                             <  1:201126/139‑1 (MQ=255)
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCgcgc                                                                                                                            <  1:168889/140‑1 (MQ=255)
aGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCgcgc                                                                                                                            <  1:50992/140‑1 (MQ=255)
   cGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCgctgc                                                                                                                         >  2:327782/1‑140 (MQ=255)
                cGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCcat                                                                                                            <  1:76980/140‑1 (MQ=255)
                 gTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCcatc                                                                                                           <  2:192696/140‑1 (MQ=255)
                 gTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCcatc                                                                                                           <  1:37498/140‑1 (MQ=255)
                          tGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTAt                                                                                                  >  1:249200/1‑140 (MQ=255)
                           gCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATc                                                                                                 >  2:188164/1‑140 (MQ=255)
                             gTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCtt                                                                                               <  1:123342/140‑1 (MQ=255)
                                            gtACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGc                                                                                          <  1:323675/130‑1 (MQ=255)
                                            gtACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGc                                                                                          >  2:323675/1‑130 (MQ=255)
                                                  cAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGcc                                                                              >  2:106554/1‑136 (MQ=255)
                                                  cAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGcc                                                                              <  1:106554/136‑1 (MQ=255)
                                                    aTCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCcatcat                                                                        >  2:281367/1‑140 (MQ=255)
                                                         tcctGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGa                                                                   <  2:354444/140‑1 (MQ=255)
                                                                    aaCCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATttatta                                                        <  1:327782/140‑1 (MQ=255)
                                                                                         gtaCCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTcgacg                                   >  1:184790/1‑140 (MQ=255)
                                                                                             cAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATtc                               >  1:311052/1‑140 (MQ=255)
                                                                                                     aCGTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTggg                       >  2:259835/1‑140 (MQ=255)
                                                                                                       gTTCGGCGTAGAAACAGCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGt                     <  2:311052/140‑1 (MQ=255)
                                                                                                                       gCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCa                   >  2:141176/1‑126 (MQ=255)
                                                                                                                       gCCTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCa                   <  1:141176/126‑1 (MQ=255)
                                                                                                                         cTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGa   >  2:249691/1‑140 (MQ=255)
                                                                                                                         cTGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGa   >  2:268326/1‑140 (MQ=255)
                                                                                                                          tGCAGTGGGCACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGaa  <  1:281367/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
AGGCGTCACGGGGATACGTCAGGGGCTGCGTTCGTGCTCCTGTTGTACGGCAATCCTTCCTGCTGATGAACCCGTGTGCCCGCGTTGTAGTACCAAAGGGTACGTTCGGCGTAGAAACAGCCTGCAGTGGACACTCGCGCTGCTTGTAACGTCCATCATGCTGTATCTTCCGGCTAATATTTTGCCCATCATGGTGACGGATTTATTAGGCTCGAAGATGCCGTCGACGATTCTCGCTGGGGTCATTCTGTTATGGAGCGAA  >  NC_000913/1012636‑1012897

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