Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,644,276 A→G V10A (GTG→GCG)  ydfU ← Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,644,2760AG98.0% 172.5 / NA 51V10A (GTG→GCG) ydfUQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/25);  minor base C (1/0);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTGAGAGCTTTCTT  >  NC_000913/1644130‑1644416
                                                                                                                                                  |                                                                                                                                            
ccGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAg                                                                                                                                          <  1:724540/151‑1 (MQ=255)
 cGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGcc                                                                                                                                            >  1:592996/1‑148 (MQ=255)
  gCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAg                                                                                                                                          >  2:51399/1‑149 (MQ=255)
     aGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAg                                                                                                                                     <  1:624533/151‑1 (MQ=255)
          cTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCt                                                                                                                                 >  2:40285/1‑150 (MQ=255)
             gCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGt                                                                                                                             <  2:371852/151‑1 (MQ=255)
             gCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGt                                                                                                                             <  2:783762/151‑1 (MQ=255)
              cGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTa                                                                                                                            <  2:8450/151‑1 (MQ=255)
                 aCGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAAt                                                                                                                          >  2:727476/1‑150 (MQ=255)
                      aGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGc                                                                                                                    >  2:111722/1‑151 (MQ=255)
                       ggAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCa                                                                                                                   <  2:616346/151‑1 (MQ=255)
                         aaCGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACg                                                                                                                 <  1:569831/151‑1 (MQ=255)
                          aCGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCTGAACAGGTCTAAGTAATACCCGCACGt                                                                                                                <  2:426762/151‑1 (MQ=255)
                            gACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCCCCAGAACAGGTCTAAGTAATACCCGCACGtt                                                                                                               >  1:215026/1‑150 (MQ=255)
                                 cGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTcc                                                                                                         <  1:718328/151‑1 (MQ=255)
                                         aGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGt                                                                                                 >  2:253206/1‑151 (MQ=255)
                                          gATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGtt                                                                                                <  1:430977/151‑1 (MQ=255)
                                           aTTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTg                                                                                               <  2:302576/151‑1 (MQ=255)
                                            ttACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGc                                                                                              >  2:346296/1‑151 (MQ=255)
                                                  aTGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGaa                                                                                         >  1:219950/1‑150 (MQ=255)
                                                       ttttGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTg                                                                                     <  1:815285/149‑1 (MQ=255)
                                                               ccggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCAc                                                                           <  2:459867/151‑1 (MQ=255)
                                                                 ggttccACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACgt                                                                         >  2:380788/1‑151 (MQ=255)
                                                                       accaGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTggg                                                                   >  1:410644/1‑151 (MQ=255)
                                                                             aCTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCg                                                             >  2:596912/1‑151 (MQ=255)
                                                                              cTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCgg                                                            >  1:246751/1‑151 (MQ=255)
                                                                                 gAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAAt                                                         <  1:462495/151‑1 (MQ=255)
                                                                                   ggATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATaa                                                       <  2:14693/151‑1 (MQ=255)
                                                                                    gATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATa                                                        >  2:843957/1‑149 (MQ=255)
                                                                                     aTTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAgg                                                     >  1:581359/1‑151 (MQ=255)
                                                                                           aaaTATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCt                                               <  2:725541/151‑1 (MQ=255)
                                                                                           aaaTATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCt                                               >  2:686899/1‑151 (MQ=255)
                                                                                           aaaTATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCt                                               <  1:528365/151‑1 (MQ=255)
                                                                                            aaTATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTg                                              <  1:747326/151‑1 (MQ=255)
                                                                                                tCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGAc                                            >  2:272575/1‑149 (MQ=255)
                                                                                                  tGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACag                                          >  1:798083/1‑149 (MQ=255)
                                                                                                   gTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACagag                                        <  1:252421/150‑1 (MQ=255)
                                                                                                   gTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACag                                          >  2:439123/1‑148 (MQ=255)
                                                                                                       ggATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTa                                    <  2:28066/150‑1 (MQ=255)
                                                                                                       ggATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTAt                                   >  2:598716/1‑151 (MQ=255)
                                                                                                         aTTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATcc                                 <  2:449385/151‑1 (MQ=255)
                                                                                                          ttCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATcc                                 >  1:124775/1‑150 (MQ=255)
                                                                                                           tCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATcc                                 >  1:118761/1‑149 (MQ=255)
                                                                                                            cACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAg                               <  2:546964/150‑1 (MQ=255)
                                                                                                            cACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAg                               >  2:238369/1‑150 (MQ=255)
                                                                                                                           ggACCACCAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGAc               >  2:376776/1‑151 (MQ=255)
                                                                                                                             accaccAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTg             >  1:96077/1‑151 (MQ=255)
                                                                                                                             accaccAGCCCAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTg             <  2:841700/151‑1 (MQ=255)
                                                                                                                                     cccAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTGAGAGCttt     <  2:847697/151‑1 (MQ=255)
                                                                                                                                       cAAGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTGAGAGCTTTCt   <  2:820983/151‑1 (MQ=255)
                                                                                                                                        aaGCTCAGGCGCCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTGAGAGCTTTCtt  >  1:408946/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGGCACCAGAACAGGTCTAAGTAATACCCGCACGTTACCTCCAGATCCGTTGCTGGAAAGTGCGGGACGCACGTGGTGGGCGTTCGGAATAAGGCAGCCTGACAGAGATTATCCAGTGCCGATAGTCGAGACTGAGAGCTTTCTT  >  NC_000913/1644130‑1644416

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