Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,555,722 C→T intergenic (‑53/+104) maeA ← / ← sra malate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,7220CT98.1% 178.2 / ‑4.4 53intergenic (‑53/+104)maeA/sramalate dehydrogenase, NAD‑requiring/30S ribosomal subunit protein S22
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/31);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTT  >  NC_000913/1555585‑1555856
                                                                                                                                         |                                                                                                                                      
cTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGttt                                                                                                                                      <  2:30224/140‑1 (MQ=255)
   tATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTAc                                                                                                                                   >  1:956201/1‑140 (MQ=255)
             cGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCg                                                                                                                         >  1:772774/1‑140 (MQ=255)
              ggAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGt                                                                                                                        <  1:42521/140‑1 (MQ=255)
                aaaTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCg                                                                                                                      <  2:1145345/140‑1 (MQ=255)
                aaaTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCg                                                                                                                      <  1:193096/140‑1 (MQ=255)
                aaaTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCg                                                                                                                      <  1:345246/140‑1 (MQ=255)
                  aTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGc                                                                                                                    <  1:199617/140‑1 (MQ=255)
                  aTTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTCCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGc                                                                                                                    <  2:363367/140‑1 (MQ=255)
                     ccAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTc                                                                                                                 <  2:986550/140‑1 (MQ=255)
                     ccAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTc                                                                                                                 <  2:902740/140‑1 (MQ=255)
                     ccAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTc                                                                                                                 <  2:154117/140‑1 (MQ=255)
                     ccAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTc                                                                                                                 <  1:1095672/140‑1 (MQ=255)
                     ccAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTc                                                                                                                 <  1:1505607/140‑1 (MQ=255)
                              cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGta                                                                                                            >  1:958900/1‑136 (MQ=255)
                              cAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGta                                                                                                            <  2:958900/136‑1 (MQ=255)
                                     aGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATa                                                                                                 <  2:32095/140‑1 (MQ=255)
                                      gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaa                                                                                                >  2:896059/1‑140 (MQ=255)
                                      gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaa                                                                                                >  2:181106/1‑140 (MQ=255)
                                      gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaa                                                                                                >  1:682070/1‑140 (MQ=255)
                                      gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaa                                                                                                >  1:748045/1‑140 (MQ=255)
                                      gCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaa                                                                                                >  2:44118/1‑140 (MQ=255)
                                          aaGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaac                                                                                            <  1:26176/140‑1 (MQ=255)
                                           aGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATaaaaaca                                                                                           >  2:1449945/1‑140 (MQ=255)
                                                 ataAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAaaaacaaaaaca                                                                                     >  2:580688/1‑140 (MQ=255)
                                                     aGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCAt                                                                                 <  1:44118/140‑1 (MQ=255)
                                                        gAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                                              <  1:787189/140‑1 (MQ=255)
                                                        gAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgc                                                                              <  1:1449945/140‑1 (MQ=255)
                                                         aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgcg                                                                             >  1:1479676/1‑140 (MQ=255)
                                                         aaCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgcg                                                                             >  2:317950/1‑140 (MQ=255)
                                                          aCGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTgcgc                                                                            >  1:1302911/1‑140 (MQ=255)
                                                           cGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCa                                                                           <  2:523743/140‑1 (MQ=255)
                                                             cTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAgg                                                                         <  2:650225/140‑1 (MQ=255)
                                                               gTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCa                                                                       <  1:580688/140‑1 (MQ=255)
                                                                          tttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg                                                                    <  1:1468098/132‑1 (MQ=255)
                                                                          tttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg                                                                    >  2:1468098/1‑132 (MQ=255)
                                                                                          tcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTg                                            <  1:181106/140‑1 (MQ=255)
                                                                                                ttttttGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGtt                                      <  1:1220451/140‑1 (MQ=255)
                                                                                                        atatCCATCCCTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTa                              >  1:654366/1‑140 (MQ=255)
                                                                                                                cccTGGGGGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAg                      <  2:1453001/140‑1 (MQ=255)
                                                                                                                      ggggCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGc                >  2:1417181/1‑140 (MQ=255)
                                                                                                                      ggggCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGc                >  1:172890/1‑140 (MQ=255)
                                                                                                                         gCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGtt             <  2:1029551/140‑1 (MQ=255)
                                                                                                                          cTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGttt            >  2:264035/1‑140 (MQ=255)
                                                                                                                          cTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGttt            >  2:315666/1‑140 (MQ=255)
                                                                                                                          cTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGttt            >  1:896700/1‑140 (MQ=255)
                                                                                                                          cTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGttt            >  1:114983/1‑140 (MQ=255)
                                                                                                                               aTCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGc       <  1:150087/140‑1 (MQ=255)
                                                                                                                                  gTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCgg    <  1:1417181/140‑1 (MQ=255)
                                                                                                                                  gTCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCgg    <  2:748045/140‑1 (MQ=255)
                                                                                                                                   tCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCGATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGt   <  1:1223187/140‑1 (MQ=255)
                                                                                                                                   tCTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGt   >  2:1482451/1‑140 (MQ=255)
                                                                                                                                    cTTTGTTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt  <  2:1302911/140‑1 (MQ=255)
                                                                                                                                    cTTTGTTTTACCACCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt  <  2:956201/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CTTTATTCAACAACGGAAATTCCAGCAGTACAGGGCCAGCGTAAGGGATATAAAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTT  >  NC_000913/1555585‑1555856

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