Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,327,533 G→T R113S (CGT→AGT)  mreC ← cell wall structural complex MreBCD transmembrane component MreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,327,5330GT99.1% 1312.9 / ‑3.5 328R113S (CGT→AGT) mreCcell wall structural complex MreBCD transmembrane component MreC
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (159/165);  total (161/167)
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

Only 100 of 338 total aligned reads displayed.

ATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACT  >  REL606/3327460‑3327605
                                                                         |                                                                        
aTAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCa                                                                        >  2:7096192/1‑76 (MQ=255)
 tAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATACTGACTCAg                                                                       >  2:7518628/1‑76 (MQ=255)
  aacaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGc                                                                      >  1:9167530/1‑76 (MQ=255)
  aacaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGc                                                                      >  2:4340793/1‑76 (MQ=255)
  aacaacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGc                                                                      <  1:7473729/76‑1 (MQ=255)
     aacTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCggg                                                                   <  2:2682110/76‑1 (MQ=255)
        ttGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGaa                                                                <  2:4736382/76‑1 (MQ=255)
        ttGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGaa                                                                <  2:6550654/76‑1 (MQ=255)
         tGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAAc                                                               <  1:7102187/76‑1 (MQ=255)
           aTCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAAccc                                                             <  1:3220947/76‑1 (MQ=255)
           aTCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAAccc                                                             >  1:3073088/1‑76 (MQ=255)
            tCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCca                                                            <  2:1862606/76‑1 (MQ=255)
            tCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCca                                                            >  1:6212559/1‑76 (MQ=255)
            tCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCca                                                            <  2:407015/76‑1 (MQ=255)
             cGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcag                                                           <  2:5125069/76‑1 (MQ=255)
             cGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcag                                                           >  2:8017740/1‑76 (MQ=255)
             cGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcag                                                           >  2:8700725/1‑76 (MQ=255)
              gCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCcagc                                                          >  1:6235672/1‑76 (MQ=255)
               cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcagca                                                         <  1:9261838/76‑1 (MQ=255)
               cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcagca                                                         >  2:3584087/1‑76 (MQ=255)
               cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcagca                                                         >  1:1811269/1‑76 (MQ=255)
               cTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACACAGCGGGGAACCcagca                                                         >  1:4564496/1‑76 (MQ=255)
                tataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcagcag                                                        >  2:10042130/1‑76 (MQ=255)
                tataAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCcagcag                                                        >  1:2101562/1‑76 (MQ=255)
                   aaGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTc                                                     >  1:2519273/1‑76 (MQ=255)
                    aGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTcg                                                    <  2:3703439/76‑1 (MQ=255)
                    aGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTcc                                                    >  2:3327571/1‑75 (MQ=255)
                      gATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTcgcg                                                  <  1:4874774/76‑1 (MQ=255)
                       aTCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCgcgc                                                 >  1:9780308/1‑76 (MQ=255)
                        tCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCgcgcc                                                >  1:667749/1‑75 (MQ=255)
                        tCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCa                                                <  1:10090469/76‑1 (MQ=255)
                        tCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCa                                                <  1:8885249/76‑1 (MQ=255)
                        tCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCa                                                <  1:5283074/76‑1 (MQ=255)
                         cGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAg                                               >  2:5334828/1‑76 (MQ=255)
                         cGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAg                                               <  2:2189614/76‑1 (MQ=255)
                          tttAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGa                                              <  1:2937597/75‑1 (MQ=255)
                          gTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGa                                              >  1:9635207/1‑76 (MQ=255)
                          gTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCGCAGa                                              >  1:2693295/1‑76 (MQ=255)
                          gTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCGCAGa                                              >  2:2006429/1‑76 (MQ=255)
                           ttAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGAc                                             <  2:2927286/76‑1 (MQ=255)
                           ttAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCGCAGAc                                             >  2:4654184/1‑76 (MQ=255)
                            tAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCACCGAcg                                            >  2:6931260/1‑76 (MQ=255)
                             aaCCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGAcgc                                           <  1:5459915/76‑1 (MQ=255)
                             aaCCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCGCCGAcgc                                           >  2:3819790/1‑76 (MQ=255)
                               ccGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGAcgcgc                                         >  2:9497604/1‑76 (MQ=255)
                                cGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGAcgcgcg                                        <  2:418639/76‑1 (MQ=255)
                                 gTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGt                                       >  1:8142064/1‑76 (MQ=255)
                                 gTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCCCGCAGACGCGCGt                                       >  1:1626547/1‑76 (MQ=255)
                                  tGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGAcgcgcggt                                      >  2:9146764/1‑74 (MQ=255)
                                  tGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGtt                                      <  2:9291021/76‑1 (MQ=255)
                                  tGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGtt                                      <  1:4098328/76‑1 (MQ=255)
                                  tGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGtt                                      <  2:64076/76‑1 (MQ=255)
                                    gagaTAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTtct                                    <  2:5878630/76‑1 (MQ=255)
                                     agaTAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTtctc                                   <  1:9068370/76‑1 (MQ=255)
                                      gaTAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGGTTTcc                                  >  2:7400990/1‑76 (MQ=255)
                                       aTAACCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCt                                 <  2:7784810/76‑1 (MQ=255)
                                         aaCCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCCGACCCGCGCTCTCCCGt                               >  1:1053010/1‑76 (MQ=255)
                                         aaCCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGt                               >  1:7583547/1‑76 (MQ=255)
                                         aaCCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGt                               <  2:9032823/76‑1 (MQ=255)
                                         aaCCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCCGt                               >  1:136792/1‑76 (MQ=255)
                                         aaCCTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCTGt                               >  2:3932110/1‑76 (MQ=255)
                                           ccTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCCGCCGCGCGTTCTCCCGttt                             >  1:8613912/1‑76 (MQ=255)
                                           ccTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGttt                             <  2:151151/76‑1 (MQ=255)
                                           ccTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGttt                             <  2:6754423/76‑1 (MQ=255)
                                           ccTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGttt                             <  2:8794379/76‑1 (MQ=255)
                                            cTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTg                            <  1:5613359/76‑1 (MQ=255)
                                            cTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTg                            <  1:501409/76‑1 (MQ=255)
                                            cTGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCCGtttt                            >  2:9385872/1‑75 (MQ=255)
                                             tGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGt                           <  1:4392844/76‑1 (MQ=255)
                                             tGAGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGt                           <  2:1434400/76‑1 (MQ=255)
                                               aGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTAt                         <  2:8192461/76‑1 (MQ=255)
                                               aGTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCTGTTTGTAt                         >  2:7066415/1‑76 (MQ=255)
                                                gTCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTAtt                        <  2:790234/76‑1 (MQ=255)
                                                 tCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTAttt                       >  1:7914675/1‑75 (MQ=255)
                                                 tCACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTg                       <  1:59511/76‑1 (MQ=255)
                                                  cACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGt                      >  1:6419855/1‑76 (MQ=255)
                                                  cACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGt                      <  1:3525460/76‑1 (MQ=255)
                                                  cACCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGt                      <  2:3647111/76‑1 (MQ=255)
                                                   aCCATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTc                     <  2:9748353/76‑1 (MQ=255)
                                                    ccATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTcc                    <  1:9361609/76‑1 (MQ=255)
                                                     cATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCa                   >  1:7149800/1‑76 (MQ=255)
                                                     cATTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCa                   <  1:4008392/76‑1 (MQ=255)
                                                      aTTTTCTGCTCATCCTGACTCAGCGGGGAGCCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCaa                  <  1:3643395/76‑1 (MQ=255)
                                                      aTTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCCGACGCGCGTTCTCCTGTTTGTATTGTCCaa                  >  1:9907295/1‑76 (MQ=255)
                                                      aTTTTCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCaa                  >  2:10142365/1‑76 (MQ=255)
                                                       ttttCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCCCGTTCTCCTGTTTGTATTGTCCCAg                 >  2:1715848/1‑76 (MQ=255)
                                                        tttCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc                <  1:2037624/76‑1 (MQ=255)
                                                        tttCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc                <  1:8290423/76‑1 (MQ=255)
                                                        tttCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGc                >  1:679905/1‑76 (MQ=255)
                                                          tCTGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCAt              <  2:303168/76‑1 (MQ=255)
                                                           cggCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATc             <  1:1148263/74‑1 (MQ=255)
                                                            tGCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATca            <  2:9948403/76‑1 (MQ=255)
                                                             gCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcag           <  2:1358700/76‑1 (MQ=255)
                                                             gCTCATCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcag           <  1:8440926/76‑1 (MQ=255)
                                                                 aTCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATcag           <  1:8660155/72‑1 (MQ=255)
                                                                  tCCTGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGtt      <  1:375605/76‑1 (MQ=255)
                                                                     tGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCAc   <  1:2438179/76‑1 (MQ=255)
                                                                     tGACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGGTCTCCCGTTTTTATTGTCCCAGCATCAGCAGTTTAc   >  2:2702983/1‑76 (MQ=255)
                                                                      gACTCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACt  <  2:560518/76‑1 (MQ=255)
                                                                      gACTCAGCGGGGAACCCAGCAGCTCCCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATAAGCAGTTCACt  >  2:7600241/1‑76 (MQ=255)
                                                                         |                                                                        
ATAACAACTTGATCGCTATAAGGATCGTTAACCGTGGAGATAACCTGAGTCACCATTTTCTGCTCATCCTGACGCAGCGGGGAACCCAGCAGCTCGCGCAGACGCGCGTTCTCCTGTTTGTATTGTCCAAGCATCAGCAGTTCACT  >  REL606/3327460‑3327605

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: