Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,799,223 A→G 8.7% intergenic (‑194/+22) PP_2452 ← / ← ansB hypothetical protein/glutaminase‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,799,2230AG 8.7% 124.4 / 3.3 46intergenic (‑194/+22)PP_2452/ansBhypothetical protein/glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base G (2/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTT  >  NC_002947/2799077‑2799364
                                                                                                                                                  |                                                                                                                                             
agtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCg                                                                                                                                            >  1:303252/1‑150 (MQ=255)
 gtaaATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCgg                                                                                                                                           >  1:290260/1‑150 (MQ=255)
      ttCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTg                                                                                                                                      <  2:230401/150‑1 (MQ=255)
           aaCAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgggg                                                                                                                                 <  2:4436/150‑1 (MQ=255)
                 ttGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCa                                                                                                                           >  1:1153355/1‑150 (MQ=255)
                   gTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgggg                                                                                                                                 <  2:850754/142‑1 (MQ=255)
                   gTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAgggg                                                                                                                                 >  1:850754/1‑142 (MQ=255)
                     cAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                       <  1:425044/150‑1 (MQ=255)
                     cAGAAGAAATCCAGCGTTGGCCTGTAACCAATGAGCCGCGTGATAAATTGAAACGATGTCGCCGGTTTATAAGTGCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCa                                                                                                                       <  1:569814/150‑1 (MQ=255)
                        gagaAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACGCAGTa                                                                                                                    >  2:522992/1‑150 (MQ=255)
                           aaaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTc                                                                                                                 <  2:303252/150‑1 (MQ=255)
                            aaTCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTcc                                                                                                                <  1:522992/150‑1 (MQ=255)
                              tCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                              >  2:1011122/1‑150 (MQ=255)
                              tCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTCTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCa                                                                                                              >  2:289952/1‑150 (MQ=255)
                                cAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCaga                                                                                                            >  2:1111400/1‑150 (MQ=255)
                                          cTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTg                                                                                                                                      <  1:181852/114‑1 (MQ=255)
                                          cTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGGGCCGGTGCGCCGGCTCTg                                                                                                                                      >  2:181852/1‑114 (MQ=255)
                                          cTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTg                                                                                                  <  2:718027/150‑1 (MQ=255)
                                                         ccGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTc                                                                                   <  2:1153355/150‑1 (MQ=255)
                                                                gAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTc                                                                            >  2:673129/1‑150 (MQ=255)
                                                                 aaaGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCt                                                                           >  2:556397/1‑150 (MQ=255)
                                                                         aCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCTGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTCGGTCATc                                                                   <  2:397822/150‑1 (MQ=255)
                                                                         aCGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTACAGCTCTTTGCTGTCCTGGGTCTTGGTCATc                                                                   >  2:832050/1‑150 (MQ=255)
                                                                           gATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                                    <  1:243666/83‑1 (MQ=255)
                                                                           gATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAg                                                                                                                                    >  2:243666/1‑83 (MQ=255)
                                                                            aTGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCt                                                                <  2:1135787/150‑1 (MQ=255)
                                                                                         gTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGa                                                   >  1:462545/1‑150 (MQ=255)
                                                                                             ccTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCcg                                               <  2:462545/150‑1 (MQ=255)
                                                                                             ccTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCcg                                               >  1:1129524/1‑150 (MQ=255)
                                                                                              ctctGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCcgc                                              <  1:1011122/150‑1 (MQ=255)
                                                                                              ccctGTGTTTACGCGGTTGAGGGATTGGTAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCCCTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTGATCGCTACCATCGCCAGGATCcgc                                              <  2:290260/148‑1 (MQ=255)
                                                                                                     tCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCttttt                                       <  1:556397/150‑1 (MQ=255)
                                                                                                         cgcgGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCgg                                   <  1:1111400/150‑1 (MQ=255)
                                                                                                              ttGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGAACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCa                              <  1:198503/150‑1 (MQ=255)
                                                                                                               tGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAg                             >  1:962459/1‑150 (MQ=255)
                                                                                                                   gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGt                                                                             <  2:641692/98‑1 (MQ=255)
                                                                                                                   gAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGt                                                                             >  1:641692/1‑98 (MQ=255)
                                                                                                                           aaCGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCAc                 >  1:598393/1‑150 (MQ=255)
                                                                                                                                    cTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc        >  2:568366/1‑150 (MQ=255)
                                                                                                                                    cTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc        <  1:289952/150‑1 (MQ=255)
                                                                                                                                     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTg                                                                                                  <  2:358711/59‑1 (MQ=255)
                                                                                                                                     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTg                                                                                                  >  1:358711/1‑59 (MQ=255)
                                                                                                                                     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTAc                                                              >  2:427038/1‑95 (MQ=255)
                                                                                                                                     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTAc                                                              <  1:427038/95‑1 (MQ=255)
                                                                                                                                     tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCg       >  2:736241/1‑150 (MQ=255)
                                                                                                                                          gCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt  <  1:568366/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
AGTAAATTCGCAACAGCTTGTCAGGAGAAATCCAGCGTTGGCCTGTAAACCATGAGGCCGGTCAGAAAGGTAAACGATGTCGCCGGTTTGTAACCTCTGTGTCTACGCGGTTGAGGAATTGGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTT  >  NC_002947/2799077‑2799364

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