Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 568,752 A→G G71G (GGA→GGG ybcO → DUF1364 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273568,7520AG97.6% 129.8 / ‑3.8 42G71G (GGA→GGGybcODUF1364 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/26);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838
                                                                                       |                                                                                      
tACGAACCCGCCAGCCCTGATTGCCACCATTGCATGTTCCGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGgta                                                                                      <  1:314517/90‑1 (MQ=255)
tACGAAACCGCCAGACCTGATTGCCCCCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGgta                                                                                      <  1:672364/90‑1 (MQ=255)
     aaCCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCaa                                                                                 >  2:317343/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:624829/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:603378/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  1:800475/1‑90 (MQ=255)
         gCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGa                                                                             >  2:44395/1‑90 (MQ=255)
             gACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgc                                                                         <  1:719161/90‑1 (MQ=255)
             gACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgc                                                                         <  1:501007/90‑1 (MQ=255)
               ccTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATgcgc                                                                       <  2:194049/90‑1 (MQ=255)
                  gATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCt                                                                    >  2:390498/1‑90 (MQ=255)
                     tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  2:273767/90‑1 (MQ=255)
                     tGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGa                                                                 <  1:197603/90‑1 (MQ=255)
                            aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          <  1:596733/90‑1 (MQ=255)
                            aTTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATg                                                          >  1:800747/1‑90 (MQ=255)
                                 aTGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCgag                                                     >  1:433843/1‑90 (MQ=255)
                                 aTGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCgag                                                     <  2:450868/90‑1 (MQ=255)
                                      ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                <  2:109639/90‑1 (MQ=255)
                                      ctgcctgcCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAacac                                                <  1:704551/90‑1 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  2:772100/1‑90 (MQ=255)
                                        gcctgccACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAg                                              >  2:51784/1‑90 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  1:772100/90‑1 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  1:603378/90‑1 (MQ=255)
                                              cACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATc                                        <  2:223884/90‑1 (MQ=255)
                                                        tCGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAag                              <  1:624829/90‑1 (MQ=255)
                                                         cGACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAaga                             <  1:404373/90‑1 (MQ=255)
                                                          gACCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAagag                            <  2:53395/90‑1 (MQ=255)
                                                           aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAgg                           <  1:390498/90‑1 (MQ=255)
                                                           aCCGCCGCACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTGATCTGGCTGAAAGAgg                           >  1:695356/1‑90 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       >  1:174435/1‑58 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCg                                                       <  2:174435/58‑1 (MQ=255)
                                                               ccgcACGCATTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGt                       <  1:455147/90‑1 (MQ=255)
                                                                     gcaTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTaa                 <  2:710792/90‑1 (MQ=255)
                                                                     gcaTTTTGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTaa                 <  2:714582/90‑1 (MQ=255)
                                                                        ttttGTTGACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGc              <  2:199451/90‑1 (MQ=255)
                                                                          ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  1:51784/90‑1 (MQ=255)
                                                                          ttgttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGt            <  1:677917/90‑1 (MQ=255)
                                                                             ttgACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGaa         <  2:800475/90‑1 (MQ=255)
                                                                              tgACGCTGGGTATGCAAAAGAAAGGGCGCTGGAAGGTATGGCGAGAACCCCGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAAt        <  2:511519/90‑1 (MQ=255)
                                                                               gACGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATa       >  2:732988/1‑90 (MQ=255)
                                                                                aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATAc      >  1:675328/1‑90 (MQ=255)
                                                                                aCGCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAAAACAGGTTATCTGGCTGAAAGGGGGGGTTATTAAGGCGTGAATAc      >  2:372590/1‑90 (MQ=255)
                                                                                  gCTGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCt    <  2:466870/90‑1 (MQ=255)
                                                                                    tGGGTATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAc  <  1:295259/90‑1 (MQ=255)
                                                                                       |                                                                                      
TACGAAACCGCCAGACCTGATTGCCACCATTGCATGTTCTGCCTGCCACGACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTAC  >  NZ_CP009273/568665‑568838

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