Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,640,376 T→C G70G (GGA→GGG)  essQ ← Qin prophage; putative S lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,640,3760TC100.0% 149.1 / NA 41G70G (GGA→GGG) essQQin prophage; putative S lysis protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/21);  total (20/21)

GGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCA  >  NC_000913/1640243‑1640512
                                                                                                                                     |                                                                                                                                        
ggCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTActcccctctc                                                                                                                                    >  2:581598/1‑140 (MQ=255)
 gCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCg                                                                                                                                   <  1:163268/140‑1 (MQ=255)
  cGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGc                                                                                                                                  <  2:197973/140‑1 (MQ=255)
  cGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGc                                                                                                                                  <  1:742004/140‑1 (MQ=255)
     cATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGc                                                                                                                                  >  1:417929/1‑137 (MQ=255)
     cATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGc                                                                                                                                  <  2:417929/137‑1 (MQ=255)
     cATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGc                                                                                                                               >  1:287601/1‑140 (MQ=255)
     cATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGc                                                                                                                               >  2:39429/1‑140 (MQ=255)
            gTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGc                                                                                                                        <  1:667456/140‑1 (MQ=255)
                aaaaTTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGt                                                                                                                    <  1:581598/140‑1 (MQ=255)
                aaaaTTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGt                                                                                                                    <  1:5688/140‑1 (MQ=255)
                 aaaTTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGtc                                                                                                                   <  2:287601/140‑1 (MQ=255)
                         ccGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTaa                                                                                                           <  1:503906/140‑1 (MQ=255)
                                             cTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCC‑TTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCa                                                                                      <  2:1025690/140‑1 (MQ=255)
                                                tAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTGGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGa                                                                                    >  2:946948/1‑140 (MQ=255)
                                                        cAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTcag                                                                            >  2:993290/1‑140 (MQ=255)
                                                          gTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCa                                                                                      >  2:756366/1‑128 (MQ=255)
                                                          gTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCa                                                                                      <  1:756366/128‑1 (MQ=255)
                                                           ttACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTcagcag                                                                         >  2:294460/1‑140 (MQ=255)
                                                             aCCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGc                                                                       >  1:373954/1‑140 (MQ=255)
                                                              ccGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGcc                                                                      >  1:586815/1‑140 (MQ=255)
                                                                            tcAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTAc                                                        <  1:946948/140‑1 (MQ=255)
                                                                                aaTTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAg                                                    <  1:993290/140‑1 (MQ=255)
                                                                                     cGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCacac                                               <  2:373954/140‑1 (MQ=255)
                                                                                          tCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTAc                                                        >  1:590417/1‑126 (MQ=255)
                                                                                          tCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTAc                                                        <  2:590417/126‑1 (MQ=255)
                                                                                           ccTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTg                                         <  1:275587/140‑1 (MQ=255)
                                                                                             tttCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCa                                       <  2:423734/140‑1 (MQ=255)
                                                                                                      cAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGcc                                    >  1:819974/1‑134 (MQ=255)
                                                                                                      cAGTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGcc                                    <  2:819974/134‑1 (MQ=255)
                                                                                                        gTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGa                            >  2:201164/1‑140 (MQ=255)
                                                                                                        gTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGa                            >  1:718317/1‑140 (MQ=255)
                                                                                                        gTTCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGa                            >  1:57502/1‑140 (MQ=255)
                                                                                                         ttCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCa                                  >  2:386201/1‑133 (MQ=255)
                                                                                                         ttCATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCa                                  <  1:386201/133‑1 (MQ=255)
                                                                                                           cATAGTCTTGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACgg                         >  1:63451/1‑140 (MQ=255)
                                                                                                                  ttGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGcc                                    >  2:557446/1‑122 (MQ=255)
                                                                                                                  ttGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGcc                                    <  1:557446/122‑1 (MQ=255)
                                                                                                                  ttGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACt                  >  2:997015/1‑140 (MQ=255)
                                                                                                                  ttGAGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACt                  >  2:841216/1‑140 (MQ=255)
                                                                                                                     aGTCATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTa               >  1:1021533/1‑140 (MQ=255)
                                                                                                                        cATTGGATTACTCCCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCg            <  2:41317/140‑1 (MQ=255)
                                                                                                                                  ctcccctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCa  <  1:997015/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                        
GGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTACCGTAACTTTATTCTCAAAATTACGGACTCCTTTCACAACCAGTTCATAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCA  >  NC_000913/1640243‑1640512

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: ‑