Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 220,255 Δ1 bp 100% coding (674/816 nt) metQ ← DL‑methionine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913220,2520G.52.9% 100.6 / 150.6 68coding (677/816 nt)metQDL‑methionine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (14/18);  new base . (16/20);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

GCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCG  >  NC_000913/220153‑220348
                                                                                                   |                                                                                                
gctgctTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTtacgtacg                                                                                                  >  2:14180/1‑100 (MQ=255)
 ctgctTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACgg                                                                                                 <  1:1925521/100‑1 (MQ=255)
  tgctTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACggg                                                                                                >  2:196587/1‑100 (MQ=255)
  tgctTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACggg                                                                                                <  2:1635574/100‑1 (MQ=255)
   gctTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACgggg                                                                                               <  2:599611/100‑1 (MQ=255)
       cGTAAACTTCGTCAGACTGATAAGCCTTGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACtggactc                                                                                            >  1:370531/1‑95 (MQ=255)
       cGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGActct                                                                                          <  2:1911/100‑1 (MQ=255)
        gTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGActct                                                                                          >  2:166447/1‑100 (MQ=255)
          aaaCTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCttt                                                                                        >  2:1213180/1‑100 (MQ=255)
          aaaCTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTa                                                                                       <  2:970421/100‑1 (MQ=255)
          aaaCTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTa                                                                                       <  1:1495919/100‑1 (MQ=255)
            aCTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATc                                                                                     <  2:2393044/100‑1 (MQ=255)
             cTTCGTCAGACTGTTAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATc                                                                                     >  2:1760133/1‑100 (MQ=255)
             cTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATc                                                                                     <  2:399883/100‑1 (MQ=255)
             cTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCt                                                                                    >  2:877754/1‑100 (MQ=255)
                cGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTc                                                                                  >  2:358265/1‑100 (MQ=255)
                cGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTca                                                                                 >  1:1222081/1‑100 (MQ=255)
                cGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTca                                                                                 <  2:2266406/100‑1 (MQ=255)
                 gTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGTCTCTTTATCttga                                                                                 <  2:1756964/100‑3 (MQ=255)
                 gTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTcaa                                                                                >  1:2033739/1‑100 (MQ=255)
                 cttccgatcttAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTca                                                                                 <  2:1053740/90‑1 (MQ=255)
                  tCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTcaa                                                                                <  2:1555111/100‑1 (MQ=255)
                   cAGACTGATAAGCCTGTACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTcaaca                                                                              <  1:669764/100‑1 (MQ=255)
                      acctaTAAGCCTGGACGAATTTCTTCTCGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTGCGTACGGGACTCTTTATCTTcaaca                                                                              <  1:1498295/93‑1 (MQ=255)
                          aTAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATAc                                                                       >  2:1457237/1‑100 (MQ=255)
                           tAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTcaagatcggaag                                                                       >  1:1053737/1‑91 (MQ=255)
                             aGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCg                                                                     >  2:823100/1‑100 (MQ=255)
                              gCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGt                                                                    <  2:400319/100‑1 (MQ=255)
                               ccTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCt                                                                  >  1:2287930/1‑100 (MQ=255)
                               ccTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCt                                                                  >  1:963526/1‑100 (MQ=255)
                                cTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCtt                                                                 >  1:478930/1‑100 (MQ=255)
                                   gACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCg                                                              <  2:2397790/100‑1 (MQ=255)
                                     cGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGc                                                             <  2:1158641/100‑1 (MQ=255)
                                     cGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGcc                                                            >  1:1842349/1‑100 (MQ=255)
                                     cGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGcc                                                            >  1:1821969/1‑100 (MQ=255)
                                     cGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGcc                                                            >  1:1220471/1‑100 (MQ=255)
                                       aaTTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGAACAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCgg                                                          >  1:1434738/1‑100 (MQ=255)
                                         tttCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGCCTTTCGCCGGAg                                                        <  2:1699155/100‑1 (MQ=255)
                                          ttcttcACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAg                                                        >  2:1233663/1‑100 (MQ=255)
                                            cttcACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAAGTTTACGTACGGCGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTc                                                      <  1:1196189/100‑1 (MQ=255)
                                             ttcACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAg                                                    >  1:1053969/1‑100 (MQ=255)
                                               cACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAgg                                                   <  2:635853/100‑1 (MQ=255)
                                                aCGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGc                                                  >  1:1208001/1‑100 (MQ=255)
                                                aCGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGcc                                                 <  1:2017914/100‑1 (MQ=255)
                                                 cGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGcc                                                 <  2:1001774/100‑1 (MQ=255)
                                                  gTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGcc                                                 <  1:783633/98‑1 (MQ=255)
                                                  gTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCaa                                               >  2:1925267/1‑100 (MQ=255)
                                                     ctcGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATc                                             <  1:2249193/100‑1 (MQ=255)
                                                     ctcGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATct                                            <  2:1809143/100‑1 (MQ=255)
                                                      tcGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATct                                            >  1:1779966/1‑100 (MQ=255)
                                                       cGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctg                                           <  2:1895766/100‑1 (MQ=255)
                                                        ggCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggc                                         <  2:700447/100‑1 (MQ=255)
                                                          cGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATctggctn                                       <  2:529098/100‑2 (MQ=255)
                                                            tCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCtggctggc                                     >  2:566590/1‑100 (MQ=255)
                                                             cTTTGTTATTTTCAATCGTCACGATCATATTTACGTACGAGACTCTTTATTTTCAACAAAGTTACCGTCTTTCGCCGCAGTCAGGCCAATctggctggta                                    <  1:1357996/100‑3 (MQ=255)
                                                               ttGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCAt                                   >  2:1594228/1‑100 (MQ=255)
                                                                 gTTATCTTCACGCGTGACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAg                                 <  2:2555574/100‑1 (MQ=255)
                                                                  ttATCTTCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggt                               >  1:297442/1‑100 (MQ=255)
                                                                   tATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggt                               >  1:2449739/1‑100 (MQ=255)
                                                                   tATCTTCACGCGTCACGATCAGGTTTACGTACGGGGAATCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAagt                               >  1:109522/1‑100 (MQ=255)
                                                                      cTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAggtggt                            >  2:1259303/1‑100 (MQ=255)
                                                                       ttCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAt                          <  1:1287417/100‑1 (MQ=255)
                                                                        tCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAt                          <  1:1045030/100‑1 (MQ=255)
                                                                        tCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAt                          <  1:1732355/100‑1 (MQ=255)
                                                                        tCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAtt                         >  2:54277/1‑100 (MQ=255)
                                                                        tCACGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTAtt                         <  2:1281469/100‑1 (MQ=255)
                                                                          aCGCGTCACGATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGa                       <  2:852734/100‑1 (MQ=255)
                                                                             cgTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATa                     <  2:1276205/100‑1 (MQ=255)
                                                                               tCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAAc                   <  1:1259299/100‑1 (MQ=255)
                                                                                cACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACt                  >  1:1134365/1‑100 (MQ=255)
                                                                                   gATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGcc               <  2:398230/100‑1 (MQ=255)
                                                                                   gATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGcc               <  1:897719/100‑1 (MQ=255)
                                                                                   gATCAGGTTTACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCTGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCa              <  1:694901/100‑1 (MQ=255)
                                                                                          tttACGTACGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGAt       <  1:1925257/100‑1 (MQ=255)
                                                                                            tacgtacgGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGAttt     >  2:1637089/1‑100 (MQ=255)
                                                                                              cgtacgGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTgc   >  2:1717442/1‑100 (MQ=255)
                                                                                               gtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTgc   >  1:415111/1‑100 (MQ=255)
                                                                                               gtacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCACAGGCGGTATTGATAACTGCCAGAGCGCTTTgc   >  1:56795/1‑100 (MQ=255)
                                                                                                tacgGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTgcg  >  2:771751/1‑100 (MQ=255)
                                                                                                   |                                                                                                
GCTGCTTCGTAAACTTCGTCAGACTGATAAGCCTGGACGAATTTCTTCACGTTCTCGGCGTCTTTGTTATCTTCACGCGTCACGATCAGGTTTACGTACGGGGACTCTTTATCTTCAACAAAGATACCGTCTTTCGCCGGAGTCAGGCCAATCTGGCTGGCATAGGTGGTATTGATAACTGCCAGAGCGATTTGCG  >  NC_000913/220153‑220348

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