Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 20,770 (A)7→8 intergenic (‑262/+45) insA‑1 ← / ← rpsT IS1 protein InsA/30S ribosomal subunit protein S20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091320,7631.A96.2% 177.6 / ‑0.5 52intergenic (‑255/+52)insA‑1/rpsTIS1 protein InsA/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (21/29);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAG  >  NC_000913/20624‑20894
                                                                                                                                            |                                                                                                                                    
tatTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCAt                                                                                                                                       >  2:210632/1‑140 (MQ=255)
         tCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaaaaac                                                                                                                            >  1:534755/1‑138 (MQ=255)
          cAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAacc                                                                                                                           >  1:340807/1‑140 (MQ=255)
             aaTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGc                                                                                                                         >  2:88016/1‑140 (MQ=255)
                  aTCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTgcg                                                                                                                    <  1:172836/140‑1 (MQ=255)
                  aTCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTgcg                                                                                                                    <  1:240394/140‑1 (MQ=255)
                  aTCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAAACCGCTTgcg                                                                                                                    >  1:2250/1‑140 (MQ=255)
                       ccGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGc                                                                                                               >  2:198113/1‑140 (MQ=255)
                           tGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCtttt                                                                                                           >  1:448392/1‑140 (MQ=255)
                                cTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTcaca                                                                                                      <  1:593307/140‑1 (MQ=255)
                                                aCGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTg                                                                                      >  1:766986/1‑140 (MQ=255)
                                                aCGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTg                                                                                      <  1:621636/140‑1 (MQ=255)
                                                        ttCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTaa                                                                              >  1:395760/1‑140 (MQ=255)
                                                          cAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGc                                                                            <  1:198809/140‑1 (MQ=255)
                                                            gATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCa                                                                          <  2:774690/140‑1 (MQ=255)
                                                            gATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCa                                                                          <  1:770451/140‑1 (MQ=255)
                                                            gATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCa                                                                          <  2:487299/140‑1 (MQ=255)
                                                            gATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCa                                                                          <  1:621045/140‑1 (MQ=255)
                                                            gATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCa                                                                          <  2:245194/140‑1 (MQ=255)
                                                                  aCAACGGCATGTCTCTTTCGGCTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAACTTCAGCAAATTGGCGATTAAGCCCGAttgt                                                                    >  2:436284/1‑140 (MQ=255)
                                                                    aaCGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTttgttg                                                                  >  2:632504/1‑140 (MQ=255)
                                                                     aCGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGa                                                                 <  1:198113/140‑1 (MQ=255)
                                                                      cGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGAt                                                                >  2:84451/1‑140 (MQ=255)
                                                                         cATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCtg                                                             >  2:39167/1‑140 (MQ=255)
                                                                         cATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCtg                                                             >  1:242689/1‑140 (MQ=255)
                                                                           tGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGa                                                                 >  2:343605/1‑134 (MQ=255)
                                                                           tGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGa                                                                 <  1:343605/134‑1 (MQ=255)
                                                                           tGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCtgtg                                                           >  1:117502/1‑140 (MQ=255)
                                                                           tGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCtgtg                                                           <  1:323091/140‑1 (MQ=255)
                                                                              ctctTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAg                                                        <  1:39167/140‑1 (MQ=255)
                                                                              ctctTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAg                                                        <  2:448392/140‑1 (MQ=255)
                                                                                    tGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCTGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGt                                                  <  2:554142/140‑1 (MQ=255)
                                                                                    tGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGt                                                  <  1:436103/140‑1 (MQ=255)
                                                                                     gACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGtt                                                 <  1:210632/140‑1 (MQ=255)
                                                                                      aCTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTa                                                <  2:82200/140‑1 (MQ=255)
                                                                                      aCTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTa                                                <  2:2250/140‑1 (MQ=255)
                                                                                       cTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAg                                               >  1:449976/1‑140 (MQ=255)
                                                                                          tCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAGACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCt                                            >  2:354543/1‑140 (MQ=255)
                                                                                           cTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCtt                                           >  1:810857/1‑140 (MQ=255)
                                                                                                   cGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTg                                   <  2:659056/140‑1 (MQ=255)
                                                                                                    ggCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGc                                  >  2:224249/1‑140 (MQ=255)
                                                                                                     gCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCa                                 >  1:710404/1‑140 (MQ=255)
                                                                                                         tgtgtTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCtt                             <  1:632504/140‑1 (MQ=255)
                                                                                                                 ctcACGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTtg                     <  2:710404/140‑1 (MQ=255)
                                                                                                                     cGCATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGa                 >  2:14268/1‑140 (MQ=255)
                                                                                                                       cATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCTAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATc               <  2:2135/140‑1 (MQ=255)
                                                                                                                       cATCACAAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATc               <  2:340807/140‑1 (MQ=255)
                                                                                                                          cacaAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGa            >  2:375241/1‑140 (MQ=255)
                                                                                                                            caAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGAcc          >  2:750420/1‑140 (MQ=255)
                                                                                                                            caAAAGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGAcc          >  2:15697/1‑140 (MQ=255)
                                                                                                                               aaGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCttt       <  2:395760/140‑1 (MQ=255)
                                                                                                                               aaGCAGCAGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCttt       <  1:84451/140‑1 (MQ=255)
                                                                                                                                 gcagcaGGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTag     <  1:522584/140‑1 (MQ=255)
                                                                                                                                   agcagGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTagca   <  2:242689/140‑1 (MQ=255)
                                                                                                                                    gcagGCATAAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTagcag  >  1:223396/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                    
TATTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCCAGTTTGTTGATCTGTGCAGTCAGGTTAGCCTTATGACGTGCAGCTTTGTTTTTGTGGATCAGACCTTTAGCAG  >  NC_000913/20624‑20894

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: