Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,190,350 C→T 8.0% noncoding (138/161 nt) sroH ← novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,190,3500CT 8.0% 156.3 / 3.4 50noncoding (138/161 nt)sroHnovel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (2/2);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GTAACTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTC  >  NC_000913/4190215‑4190492
                                                                                                                                       |                                                                                                                                              
gTAACTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAAt                                                                                                                                      >  4:75999/1‑146 (MQ=255)
gTAACTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAAt                                                                                                                                      <  3:75999/146‑1 (MQ=255)
 tAACTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTAtc                                                                                                                                  <  3:228853/149‑1 (MQ=255)
    cTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATctct                                                                                                                               <  3:200943/149‑1 (MQ=255)
                ccATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTAtt                                                                                                                         <  4:336016/143‑1 (MQ=255)
                ccATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTAtt                                                                                                                         >  3:336016/1‑143 (MQ=255)
                       tGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTa                                                                                                            >  1:161675/1‑149 (MQ=255)
                                          tCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCt                                                                                         >  2:471108/1‑149 (MQ=255)
                                           cAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCtt                                                                                        <  1:471108/149‑1 (MQ=255)
                                              tGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCaga                                                                                                                                          >  2:384820/1‑96 (MQ=255)
                                              tGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCaga                                                                                                                                          <  1:384820/96‑1 (MQ=255)
                                              tGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCg                                                                                     >  2:429104/1‑149 (MQ=255)
                                                             ttGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGc                                                                      <  3:450228/149‑1 (MQ=255)
                                                             ttGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGc                                                                      >  4:450228/1‑149 (MQ=255)
                                                               gCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCgg                                                                    <  1:429104/149‑1 (MQ=255)
                                                                  tACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTc                                                                                                                        >  3:433242/1‑94 (MQ=255)
                                                                  tACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTc                                                                                                                        <  4:433242/94‑1 (MQ=255)
                                                                        ttATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGa                                                           >  3:74163/1‑149 (MQ=255)
                                                                        ttATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGa                                                           >  1:272281/1‑149 (MQ=255)
                                                                         tATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGaa                                                          <  1:442219/149‑1 (MQ=255)
                                                                          aTTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAAt                                                         <  1:38920/149‑1 (MQ=255)
                                                                                aaaTAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATttcttc                                                   >  1:462588/1‑149 (MQ=255)
                                                                                     agcagcTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTct                                              <  2:161675/149‑1 (MQ=255)
                                                                                        agcTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTctgct                                           <  2:367619/149‑1 (MQ=255)
                                                                                                     cTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACACCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGa                              >  4:227559/1‑149 (MQ=255)
                                                                                                      tCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGAt                             >  2:370631/1‑149 (MQ=255)
                                                                                                      tCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGAt                             <  4:14477/149‑1 (MQ=255)
                                                                                                        aaaaTAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATttc                                                      <  3:117987/122‑1 (MQ=255)
                                                                                                        aaaaTAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATttc                                                      >  4:117987/1‑122 (MQ=255)
                                                                                                             aaCAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACa                      <  2:272281/149‑1 (MQ=255)
                                                                                                              aCAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACAt                     >  4:160033/1‑149 (MQ=255)
                                                                                                               cAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTctgct                                           <  3:462963/126‑1 (MQ=255)
                                                                                                               cAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTctgct                                           >  4:462963/1‑126 (MQ=255)
                                                                                                                gaaaCCCACCTTAAGGTGGGTTTTGCCAGAGAATTATCTCTGGTATTCAGAACGcc                                                                                                                <  2:63071/55‑1 (MQ=255)
                                                                                                                gaaaCCCACCTTAAGGTGGGTTTTGCCAGAGAATTATCTCTGGTATTCAGAACGcc                                                                                                                >  1:63071/2‑56 (MQ=255)
                                                                                                                 aaaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGAcc                                                                                          >  3:226252/1‑77 (MQ=255)
                                                                                                                 aaaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGAcc                                                                                          <  4:226252/77‑1 (MQ=255)
                                                                                                                  aaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTctgc                                            >  1:399965/1‑122 (MQ=255)
                                                                                                                  aaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTctgc                                            <  2:399965/122‑1 (MQ=255)
                                                                                                                  aaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTg                                   >  1:368460/1‑131 (MQ=255)
                                                                                                                  aaCCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTg                                   <  2:368460/131‑1 (MQ=255)
                                                                                                                   aCCCACCTTAAGGTGGGTTTTGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGAc                                                                                           >  2:309106/1‑74 (MQ=255)
                                                                                                                   aCCCACCTTAAGGTGGGTTTTGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGAc                                                                                           <  1:309106/74‑1 (MQ=255)
                                                                                                                     ccACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCt              >  2:232513/1‑149 (MQ=255)
                                                                                                                     ccACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCt              >  2:62783/1‑149 (MQ=255)
                                                                                                                       aCCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTa            >  4:306976/1‑149 (MQ=255)
                                                                                                                          ttAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCt                                                             <  1:50449/97‑1 (MQ=255)
                                                                                                                          ttAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCt                                                             >  2:50449/1‑97 (MQ=255)
                                                                                                                                tGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGtct   >  3:405714/1‑149 (MQ=255)
                                                                                                                                 gggTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGtctc  >  3:202390/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                              
GTAACTAGCGCTTATCCCATTTCTGAAAAACAGAAAGGTAAATCAGTGAAATTTTTCACCATTGCCTACAACTTATTGAGAAATAAGCAGCTTCCTCAGCCCTCAAAATAACAAAACCCACCTTAAGGTGGGTTTCGCCAGAGAATTATCTCTGGTATTCAGAACGCCATTACCGGACTTTGCCTTGACCTTGCGATAATCGCAGGTTGCGGGATGTCTGAATTTCTTCAGTCTGCTGCATCCTGGAAGATGAGAACATGTGTTCTTATTTTCGTCTC  >  NC_000913/4190215‑4190492

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑