Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,328,651 C→T V110I (GTT→ATT)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,328,6510CT98.7% 949.5 / 1.4 234V110I (GTT→ATT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (119/112);  total (121/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 245 total aligned reads displayed.

CCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCAC  >  REL606/3328578‑3328723
                                                                         |                                                                        
ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgc                                                                        >  2:5027332/1‑76 (MQ=255)
ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgc                                                                        >  2:6852410/1‑76 (MQ=255)
ccAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgc                                                                        >  2:3891762/1‑76 (MQ=255)
  aGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgcgc                                                                      >  1:4035717/1‑76 (MQ=255)
  aGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgcgc                                                                      >  2:2761188/1‑76 (MQ=255)
   gcgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAAtgcgcg                                                                     <  2:7307086/76‑1 (MQ=255)
    cgcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCgg                                                                    <  1:1383937/76‑1 (MQ=255)
     gcCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCggg                                                                   >  2:4401214/1‑76 (MQ=255)
      cccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGc                                                                  >  1:3335627/1‑76 (MQ=255)
       ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCt                                                                 <  2:732880/76‑1 (MQ=255)
       ccTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCt                                                                 >  1:5323373/1‑76 (MQ=255)
          gcgcGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTgg                                                              <  1:5752415/76‑1 (MQ=255)
           cgcgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGa                                                             >  2:1422963/1‑76 (MQ=255)
           cgcgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGa                                                             >  1:805012/1‑76 (MQ=255)
           cgcgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGa                                                             <  1:3187723/76‑1 (MQ=255)
             cgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACg                                                           >  2:7089978/1‑76 (MQ=255)
             cgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACg                                                           >  2:7727160/1‑76 (MQ=255)
             cgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACg                                                           <  2:802740/76‑1 (MQ=255)
             cgGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACg                                                           <  1:40801/76‑1 (MQ=255)
              ggATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGc                                                          >  1:879106/1‑76 (MQ=255)
              ggATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGc                                                          >  2:6901493/1‑76 (MQ=255)
               gATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCa                                                         >  1:7224071/1‑76 (MQ=255)
               gATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCa                                                         >  2:3823857/1‑76 (MQ=255)
                aTTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCAt                                                        <  1:7918176/76‑1 (MQ=255)
                aTTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCAt                                                        >  1:6578439/1‑76 (MQ=255)
                aTTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCAt                                                        >  1:376585/1‑76 (MQ=255)
                 ttCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGAATa                                                       >  2:1638499/1‑76 (MQ=255)
                  tCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaa                                                      <  1:3937942/76‑1 (MQ=255)
                   cACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaaa                                                     >  1:1994531/1‑76 (MQ=255)
                   cACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaaa                                                     <  1:1973471/76‑1 (MQ=255)
                   cACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaaa                                                     <  1:2981416/76‑1 (MQ=255)
                   cACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaaa                                                     <  1:681139/76‑1 (MQ=255)
                   cACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATaaa                                                     <  2:5340818/76‑1 (MQ=255)
                    aCGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAg                                                    >  2:1374126/1‑76 (MQ=255)
                       aaTTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTg                                                 >  1:5810986/1‑76 (MQ=255)
                        aTTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGt                                                >  2:3576127/1‑76 (MQ=255)
                          tGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTg                                              <  2:461746/76‑1 (MQ=255)
                              cggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCtgt                                          >  1:6349145/1‑76 (MQ=255)
                              cggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCtgt                                          >  2:3603767/1‑76 (MQ=255)
                               ggcgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCtgtg                                         <  2:5883666/76‑1 (MQ=255)
                                 cgTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCa                                       <  1:260344/76‑1 (MQ=255)
                                  gTTCAACCTGGTTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCAc                                      <  1:1069840/76‑1 (MQ=255)
                                  gTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCAc                                      >  1:4958174/1‑76 (MQ=255)
                                   ttCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACt                                     >  1:7654324/1‑76 (MQ=255)
                                   ttCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACt                                     >  1:4359714/1‑76 (MQ=255)
                                      aaCCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt                                  >  2:6331861/1‑76 (MQ=255)
                                      aaCCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt                                  >  2:2946576/1‑76 (MQ=255)
                                      aaCCTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt                                  >  2:4540675/1‑76 (MQ=255)
                                        ccTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGttt                                <  1:646989/76‑1 (MQ=255)
                                        ccTGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGttt                                <  2:219437/76‑1 (MQ=255)
                                         cTGGGTCACGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtg                               <  1:4281955/76‑1 (MQ=255)
                                          tGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtga                              >  1:2764882/1‑76 (MQ=255)
                                          tGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtga                              <  1:1580486/76‑1 (MQ=255)
                                          tGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtga                              >  2:1854335/1‑76 (MQ=255)
                                          tGGGTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtga                              <  1:757508/76‑1 (MQ=255)
                                            ggTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                            >  2:4899956/1‑76 (MQ=255)
                                            ggTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtgatg                            >  2:3501410/1‑76 (MQ=255)
                                             gTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtgatga                           <  1:526311/76‑1 (MQ=255)
                                             gTCGCGCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTtgatga                           >  2:6789687/1‑76 (MQ=255)
                                               cgcgCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAg                         <  1:4389987/76‑1 (MQ=255)
                                               cgcgCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAg                         <  2:307522/76‑1 (MQ=255)
                                               cgcgCCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAg                         >  2:2315283/1‑76 (MQ=255)
                                                 cgcCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTg                       >  2:7600227/1‑76 (MQ=255)
                                                 cgcCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTg                       <  2:5763037/76‑1 (MQ=255)
                                                 cgcCAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTg                       <  2:7874525/76‑1 (MQ=255)
                                                   ccAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCt                     >  1:128559/1‑76 (MQ=255)
                                                    cAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTg                    <  2:5651880/76‑1 (MQ=255)
                                                    cAACCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTg                    <  2:2071653/76‑1 (MQ=255)
                                                     aaCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTgg                   <  2:5117052/76‑1 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  >  1:7431533/1‑76 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  <  1:6275066/76‑1 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  <  1:5897271/76‑1 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  <  1:464522/76‑1 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  <  1:5116589/76‑1 (MQ=255)
                                                      aCCGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGa                  <  2:313244/76‑1 (MQ=255)
                                                       ccGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                 <  2:7466503/76‑1 (MQ=255)
                                                       ccGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                 <  1:2875543/76‑1 (MQ=255)
                                                       ccGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAg                 <  2:5503885/76‑1 (MQ=255)
                                                        cGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGc                >  2:3982333/1‑76 (MQ=255)
                                                        cGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGc                <  1:2851080/76‑1 (MQ=255)
                                                        cGGCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGc                >  2:5602843/1‑76 (MQ=255)
                                                         ggCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCa               <  1:7153396/76‑1 (MQ=255)
                                                          gCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAt              <  2:1185419/76‑1 (MQ=255)
                                                          gCACACAAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAt              >  2:2664886/1‑76 (MQ=255)
                                                             cacaAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttt           >  2:4708550/1‑76 (MQ=255)
                                                              acaAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAttttt          >  1:1089611/1‑76 (MQ=255)
                                                               caAACCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCAtttttt         <  1:4201426/76‑1 (MQ=255)
                                                                aaaCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTc        >  2:6228160/1‑76 (MQ=255)
                                                                aaaCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTc        >  2:2609260/1‑76 (MQ=255)
                                                                aaaCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTc        <  1:6221182/76‑1 (MQ=255)
                                                                 aaCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCa       <  1:348322/76‑1 (MQ=255)
                                                                  aCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAg      >  2:6117101/1‑76 (MQ=255)
                                                                  aCCAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAg      <  2:7438202/76‑1 (MQ=255)
                                                                   ccAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGt     >  1:609245/1‑76 (MQ=255)
                                                                    cAGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTc    <  1:7400827/76‑1 (MQ=255)
                                                                     aGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCa   >  1:7794609/1‑76 (MQ=255)
                                                                     aGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCa   >  1:5448330/1‑76 (MQ=255)
                                                                     aGAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCa   >  2:5624623/1‑76 (MQ=255)
                                                                      gAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCAc  >  2:1174698/1‑76 (MQ=255)
                                                                      gAATGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCAc  <  2:495301/76‑1 (MQ=255)
                                                                         |                                                                        
CCAGCGCCCTGCGCGGATTCACGAATTGCGCGGCGTTCAACCTGGGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTTTGATGAAGTGCTGGAGCATTTTTTCAGTCAC  >  REL606/3328578‑3328723

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