Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,897,230 A→G D309G (GAT→GGT)  manX → PTS mannose transporter subunit IIAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,897,2300AG89.7% 95.8 / 5.6 39D309G (GAT→GGT) manXPTS mannose transporter subunit IIAB
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (15/20);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC  >  NZ_CP009273/1897146‑1897316
                                                                                    |                                                                                      
tttCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGa                                                                                   <  1:66137/90‑1 (MQ=255)
             aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      <  1:338433/90‑1 (MQ=255)
             aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      <  2:551420/90‑1 (MQ=255)
             aaaaGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      >  1:542363/1‑90 (MQ=255)
                  atatCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATc                                                                 <  2:5852/90‑1 (MQ=255)
                  atatCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATc                                                                 >  1:532970/1‑90 (MQ=255)
                          gCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGc                                                         >  1:428004/1‑90 (MQ=255)
                            gTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCaa                                                       >  2:729790/1‑90 (MQ=255)
                                  gAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGa                                                 <  1:265451/90‑1 (MQ=255)
                                    aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  2:430834/90‑1 (MQ=255)
                                    aaCTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATa                                               <  1:97609/90‑1 (MQ=255)
                                        gAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTa                                           <  1:380413/90‑1 (MQ=255)
                                         aaTGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTaa                                          <  2:163741/90‑1 (MQ=255)
                                         aaTGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTaa                                          <  1:839343/90‑1 (MQ=255)
                                              gcgcgGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTAt                                     <  2:330189/90‑1 (MQ=255)
                                                   gTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTgtgt                                <  2:122379/90‑1 (MQ=255)
                                                   gTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTgtgt                                <  2:139102/90‑1 (MQ=255)
                                                           cTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATc                        >  2:267128/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      >  1:202908/1‑42 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAatgat                                                                      <  2:202908/42‑1 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAg                                             >  2:442043/1‑67 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAg                                             <  1:442043/67‑1 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  2:464051/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  2:773180/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  1:89079/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  1:836064/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  1:648614/1‑90 (MQ=255)
                                                             ggAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCAc                      >  1:679042/1‑90 (MQ=255)
                                                                aGTCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCa                   <  1:184755/90‑1 (MQ=255)
                                                                  tCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGt                 >  1:652144/1‑90 (MQ=255)
                                                                  tCCGTAAGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGt                 >  1:457248/1‑90 (MQ=255)
                                                                       aaGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTca            <  1:436535/90‑1 (MQ=255)
                                                                        aGGTTTCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTcac           <  2:648614/90‑1 (MQ=255)
                                                                        aGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTcac           <  1:779055/90‑1 (MQ=255)
                                                                            ttCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACtt       >  2:296295/1‑90 (MQ=255)
                                                                            ttCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACtt       >  1:392319/1‑90 (MQ=255)
                                                                             tCCACCGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTa      <  2:552627/90‑1 (MQ=255)
                                                                                 ccGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTc  <  2:392319/90‑1 (MQ=255)
                                                                                 ccGGTCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTc  <  2:554271/90‑1 (MQ=255)
                                                                                    |                                                                                      
TTTCGGTTGATGAAAAAGATATCGAGGCGTTCAAGAAACTGAATGCGCGCGGTATTGAGCTGGAAGTCCGTAAGGTTTCCACCGATCCGAAACTGAAAATGATGGATCTGATCAGCAAAATCGATAAGTAACGTATTGTGTTGATTATCACTCAGTTTTCACACTTAAGTC  >  NZ_CP009273/1897146‑1897316

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