Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,909,799 C→T intergenic (‑8/+89) proQ ← / ← msrC RNA chaperone ProQ/L‑methionine (R)‑S‑oxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,909,7990CT77.8% 57.0 / 13.1 27intergenic (‑8/+89)proQ/msrCRNA chaperone ProQ/L‑methionine (R)‑S‑oxide reductase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (10/11);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

GTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCG  >  NZ_CP009273/1909723‑1909883
                                                                            |                                                                                    
gTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACgg                                                                         >  2:459301/1‑90 (MQ=255)
gTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACgg                                                                         >  2:307907/1‑90 (MQ=255)
  ggggAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGAc                                                                       >  2:578427/1‑90 (MQ=255)
    ggAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGt                                                                     >  2:677839/1‑90 (MQ=255)
             tCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCcaacaa                                                            >  1:702587/1‑90 (MQ=255)
                 ccAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAgcgc                                                        <  1:335905/90‑1 (MQ=255)
                 ccAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAgcgc                                                        <  2:117666/90‑1 (MQ=255)
                    gAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAgg                                                     >  1:521733/1‑90 (MQ=255)
                         gcgATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGa                                                >  2:674224/1‑90 (MQ=255)
                          cgATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGaa                                               >  1:581513/1‑90 (MQ=255)
                             ttACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg                                            <  1:677839/90‑1 (MQ=255)
                               aCTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGc                                          <  1:462767/90‑1 (MQ=255)
                                  tcttTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCTAACAAGCGCAGGCATGAACAGGCGTc                                       <  1:518082/90‑1 (MQ=255)
                                          gCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCATGCATGAACAGGCGTCATTATAAt                               >  2:186694/1‑90 (MQ=255)
                                            tATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGa                             >  2:450364/1‑90 (MQ=255)
                                             aTTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGAc                            >  2:26891/1‑90 (MQ=255)
                                             aTTCAACTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGAc                            >  1:107524/1‑90 (MQ=255)
                                             aTTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGAc                            >  2:677512/1‑90 (MQ=255)
                                                   cTTAGGTTGATTTTCCATGAAATTTTCTGATTACAACGGACTTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATc                      <  2:309212/90‑1 (MQ=255)
                                                       ggTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTa                  <  2:415361/90‑1 (MQ=255)
                                                       ggTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTa                  <  1:45232/90‑1 (MQ=255)
                                                       ggTTGATTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTa                  <  2:87167/90‑1 (MQ=255)
                                                       ggTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTa                  <  1:585586/90‑1 (MQ=255)
                                                            aTTTTCCATGAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGc             >  2:317009/1‑90 (MQ=255)
                                                                   aTGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTa      <  2:581513/90‑1 (MQ=255)
                                                                     gAAATTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATc    <  1:674224/90‑1 (MQ=255)
                                                                       aaTTTTCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCg  <  1:578427/90‑1 (MQ=255)
                                                                            |                                                                                    
GTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCG  >  NZ_CP009273/1909723‑1909883

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