Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,769,486 G→A *274* (TAG→TAA)  yfjQ → CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,769,4860GA100.0% 137.8 / NA 43*274* (TAG→TAA) yfjQCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

AGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGT  >  NC_000913/2769353‑2769617
                                                                                                                                     |                                                                                                                                   
aGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGtt                                                                                                                                <  2:28063/139‑1 (MQ=255)
aGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGtt                                                                                                                                <  1:296974/139‑1 (MQ=255)
          cAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  1:31669/139‑1 (MQ=255)
              gggCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgt                                                                                                                  >  2:958133/1‑139 (MQ=255)
                     aTGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAat                                                                                                           >  1:956233/1‑139 (MQ=255)
                         tAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTaa                                                                                                            <  2:110826/134‑1 (MQ=255)
                         tAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTaa                                                                                                            >  1:110826/1‑134 (MQ=255)
                         tAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCg                                                                                                       <  1:944521/139‑1 (MQ=255)
                                 gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcaccac                                                                                               >  2:316696/1‑139 (MQ=255)
                                 gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcaccac                                                                                               >  2:688801/1‑139 (MQ=255)
                                       cTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTcc                                                                                         <  2:932514/139‑1 (MQ=255)
                                               cgtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTTATATCGGACACCACCTGTCCGCATCGCt                                                                                 <  1:688801/139‑1 (MQ=255)
                                                gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTa                                                                                >  1:800264/1‑139 (MQ=255)
                                                 tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    >  2:1103706/1‑102 (MQ=255)
                                                 tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    <  1:1103706/102‑1 (MQ=255)
                                                     gtgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGc                                                                           <  2:85702/139‑1 (MQ=255)
                                                      tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                >  2:827569/1‑117 (MQ=255)
                                                      tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                <  1:827569/117‑1 (MQ=255)
                                                                     gggATGTTAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCTTCGCTATGTGCTCGTGTGCCTCAATcc                                                           <  2:266376/139‑1 (MQ=255)
                                                                              aGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGacac                                                                                                  >  1:159280/1‑91 (MQ=255)
                                                                              aGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGacac                                                                                                  <  2:159280/91‑1 (MQ=255)
                                                                                cTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                                          <  2:189711/81‑1 (MQ=255)
                                                                                cTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                                          >  1:189711/1‑81 (MQ=255)
                                                                                        tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCtttt                                           <  1:922783/136‑1 (MQ=255)
                                                                                        tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCtttt                                           >  2:922783/1‑136 (MQ=255)
                                                                                          cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    <  2:922662/61‑1 (MQ=255)
                                                                                          cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    >  1:922662/1‑61 (MQ=255)
                                                                                              gtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGt                                  <  1:958133/139‑1 (MQ=255)
                                                                                              gtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGt                                  <  1:1062096/139‑1 (MQ=255)
                                                                                               tgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGtt                                 >  1:84004/1‑139 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                                         <  1:905220/59‑1 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                                         >  2:905220/1‑59 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCGATGTGCTCGTGTGCCTCAATCCCCCGGTTATGGCTTTTAACCCCCGTTACATCTgg                         >  2:512918/1‑139 (MQ=255)
                                                                                                           gAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                                          <  2:437230/102‑1 (MQ=255)
                                                                                                           gAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                                          >  1:437230/1‑102 (MQ=255)
                                                                                                            aaaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                                     <  1:1011413/74‑1 (MQ=255)
                                                                                                            aaaCACTGCTCACGCAACCGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                                     >  2:1011413/1‑74 (MQ=255)
                                                                                                                 cTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTa                                >  1:1056115/1‑122 (MQ=255)
                                                                                                                 cTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTa                                <  2:1056115/122‑1 (MQ=255)
                                                                                                                    cTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                        <  2:33285/111‑1 (MQ=255)
                                                                                                                    cTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                        >  1:33285/1‑111 (MQ=255)
                                                                                                                         gCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATaaaaa       >  1:375206/1‑139 (MQ=255)
                                                                                                                              tGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGt  >  1:322958/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                   
AGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGT  >  NC_000913/2769353‑2769617

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