Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,076,301 A→C H22P (CAC→CCC)  sbcB → exodeoxyribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,076,3010AC96.9% 109.9 / ‑3.2 32H22P (CAC→CCC) sbcBexodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (16/15);  total (16/16)
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.82e-01

TGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTAC  >  NZ_CP009273/2076166‑2076436
                                                                                                                                       |                                                                                                                                       
tgGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       <  2:761103/138‑1 (MQ=255)
tgGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       <  2:729981/138‑1 (MQ=255)
tgGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       <  1:473603/138‑1 (MQ=255)
         aGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGa                                                                                                                             <  2:338238/139‑1 (MQ=255)
         aGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGa                                                                                                                             <  1:594894/139‑1 (MQ=255)
                        tGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGtt                                                                                                              >  2:611639/1‑139 (MQ=255)
                        tGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGtt                                                                                                              >  1:441678/1‑139 (MQ=255)
                        tGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGtt                                                                                                              >  2:687354/1‑139 (MQ=255)
                        tGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGtt                                                                                                              <  2:704702/139‑1 (MQ=255)
                          atatTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCg                                                                                                            <  1:659393/139‑1 (MQ=255)
                                               agcagAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       <  1:688712/91‑1 (MQ=255)
                                               agcagAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       >  2:688712/1‑91 (MQ=255)
                                                 cagAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGaa                                                                                     >  2:220659/1‑139 (MQ=255)
                                                                       atgatgAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGa                                                                                                                             <  1:454399/77‑1 (MQ=255)
                                                                       atgatgAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGa                                                                                                                             >  2:454399/1‑77 (MQ=255)
                                                                       atgatgAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCg                                                               <  1:397211/139‑1 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGcc                                                                                                                        >  1:626123/1‑76 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGcc                                                                                                                        <  2:626123/76‑1 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAAc                                                                  <  1:367559/130‑1 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAAc                                                                  >  2:367559/1‑130 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCt                                                         >  2:568709/1‑139 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCt                                                         >  2:382461/1‑139 (MQ=255)
                                                                             aatgaCGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCt                                                         >  2:473904/1‑139 (MQ=255)
                                                                                       aGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGc                                                 >  2:588367/1‑137 (MQ=255)
                                                                                       aGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCaa                                               <  1:170534/139‑1 (MQ=255)
                                                                                                   ccTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       <  1:537988/39‑1 (MQ=255)
                                                                                                   ccTTTTTGTTTCACGATTACGAAACCTTTGGCACGcccc                                                                                                                                       >  2:537988/1‑39 (MQ=255)
                                                                                                   ccTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCtgatg                                    <  2:265664/138‑1 (MQ=255)
                                                                                                       tttGTTTCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTAt                               >  2:623096/1‑139 (MQ=255)
                                                                                                             tCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGcc                                                                                                                        <  1:83405/44‑1 (MQ=255)
                                                                                                             tCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGcc                                                                                                                        >  2:83405/1‑44 (MQ=255)
                                                                                                             tCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCg                                                                     <  1:103162/95‑1 (MQ=255)
                                                                                                             tCACGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCg                                                                     >  2:103162/1‑95 (MQ=255)
                                                                                                                cGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGcaca                                                                                                                 >  2:399632/1‑48 (MQ=255)
                                                                                                                cGATTACGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGcaca                                                                                                                 <  1:399632/48‑1 (MQ=255)
                                                                                                                     aCGAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCccagccag                  <  1:687354/138‑1 (MQ=255)
                                                                                                                       gAAACCTTTGGCACGCCCCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAg               <  2:441678/139‑1 (MQ=255)
                                                                                                                               tGGCACGCCCCACGCGTTAGATCGCCCTGCACAGTTCGCa                                                                                                          >  2:577676/1‑40 (MQ=255)
                                                                                                                               tGGCACGCCCCACGCGTTAGATCGCCCTGCACAGTTCGCa                                                                                                          <  1:577676/40‑1 (MQ=255)
                                                                                                                                    cgcccCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTAc  <  1:623096/139‑1 (MQ=255)
                                                                                                                                    cgcccCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTAc  <  1:473904/139‑1 (MQ=255)
                                                                                                                                    cgcccCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTAc  <  1:382461/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGGCGCTAAAGCTGATTAGCACGGTGATATTTGATACTCTGGCAGACAGCAGAAATAACGGATTTAACCTAATGATGAATGACGGTAAGCAACAATCTACCTTTTTGTTTCACGATTACGAAACCTTTGGCACGCACCCCGCGTTAGATCGCCCTGCACAGTTCGCAGCCATTCGCACCGATAGCGAATTCAATGTCATCGGCGAACCCGAAGTCTTTTACTGCAAGCCCGCTGATGACTATTTACCCCAGCCAGGAGCCGTATTAATTAC  >  NZ_CP009273/2076166‑2076436

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