Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,005,083 T→C intergenic (+51/‑139) unknown → / → unknown Possible protease sohB (EC 3.4.21.‑)/Outer membrane protein A precursor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,005,0830TC96.0% 60.7 / NA 25intergenic (+51/‑139)unknown/unknownPossible protease sohB (EC 3.4.21.‑)/Outer membrane protein A precursor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (12/12);  total (13/12)
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

TTAGTGGATGAGATTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAATTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGTAACTCTTGTTGCTTGTGGTAACTTAAGTAAAGTATCTAAAG  >  1/1004850‑1005303
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            
ttaGTGGATGAGATTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATc                                                                                                                                                                                                             <  1:875439/251‑1 (MQ=255)
       aTGAGATTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTg                                                                                                                                                                                                      >  1:461655/1‑251 (MQ=255)
           gaTTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAAt                                                                                                                                                                                                  >  1:421770/1‑251 (MQ=255)
            aTTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAATTCGGATTTACTTTAATCCGTTTTGTAAtt                                                                                                                                                                                                 >  1:574510/1‑251 (MQ=255)
                      gCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGa                                                                                                                                                                                       >  1:297748/1‑251 (MQ=255)
                         atgatTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGc                                                                                                                                                                                     <  1:222096/250‑1 (MQ=255)
                                                           ggTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATAc                                                                                                                                                                                          <  1:570482/211‑1 (MQ=255)
                                                                                 gTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATc                                                                                                                                    <  1:611922/243‑1 (MQ=255)
                                                                                 gTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGtt                                                                                                                             >  1:64488/1‑250 (MQ=255)
                                                                                  tGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTc                                                                                                                            <  1:288813/250‑1 (MQ=255)
                                                                                           cGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGAt                                                                                                                    >  1:561012/1‑249 (MQ=255)
                                                                                                   aaaaaaGTCGGTAAACAGGCGGAAAAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAAAAAAATGAAAGAACATTGATATAACAACCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTc                                                                                                          <  1:737566/250‑1 (MQ=255)
                                                                                                                                                  aCGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAAc                                                                                      <  1:663293/224‑1 (MQ=255)
                                                                                                                                                                                gATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGCAACTCTTGTTGCtt                             >  1:872490/1‑251 (MQ=255)
                                                                                                                                                                                        aTCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTc                                                                                                                            >  1:824891/1‑148 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCtttt                                                                                                                                                          >  1:696123/1‑100 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTc                                                                                                               <  1:880380/143‑1 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAAc                                                                                      >  1:77781/1‑168 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGCAACTCTTGTTGCTTGTGGTAACTTa                  <  1:866680/236‑1 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGCAACTCTTGTTGCTTGTGGTAACTTAAGTAAAGTATCTaaa   >  1:691544/1‑251 (MQ=255)
                                                                                                                                                                                                          cTATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGCAACTCTTGTTGCTTGTGGTAACTTAAGTAAAGTATCTaaa   >  1:780165/1‑251 (MQ=255)
                                                                                                                                                                                                           tATTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGCAACTCTTGTTGCTTGTGGTAACTTAAGTAAAGTATCTAaag  <  1:114305/251‑1 (MQ=255)
                                                                                                                                                                                                            aTTTATACAGATTGATTAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCttat                                                              <  1:99612/190‑1 (MQ=255)
                                                                                                                                                                                                                         gattAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTc                                                                                                          >  1:158878/1‑133 (MQ=255)
                                                                                                                                                                                                                         gattAAAAATAGAAAACTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAAt                                                                                                <  1:655606/143‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                            
TTAGTGGATGAGATTATGACCAGCGATGATTTGCTATTACAGGCGATGAAAGAGAAATCGGTTATCGGCGTAAAATATGCGGTGAAGAAATCGTTATTACAAAAAGTCGGTAAACAGGCGGAAGAAAGTGCCGATAATATTGCTTTACGTTGGCTTAAGAAAAATGAAAGAACATTGATATAACATCCTATTTTTTTCGATTCTATTTATACAGATTGATTAAAAATAGAAAATTCGGATTTACTTTAATCCGTTTTGTAATTACCATACGGAGCTGAAATTGAAGTTGCACTATTGCTTTTTGAGCCTTCTTTTTATTATATCTATGGTTCCTTTGGATTGTTCTACTCACCTAAAAATAAGGTAAAACAGATGAAATTATCACATTTCTTATTATCTGCAATCGCAGCGGTAACTCTTGTTGCTTGTGGTAACTTAAGTAAAGTATCTAAAG  >  1/1004850‑1005303

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