Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,162,568 C→T 9.9% intergenic (+515/‑65) PP_2776 → / → PP_2777 homocysteine S‑methyltransferase family protein/acyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,162,5680CT 9.9% 102.5 / 2.0 41intergenic (+515/‑65)PP_2776/PP_2777homocysteine S‑methyltransferase family protein/acyl carrier protein
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base T (2/2);  total (19/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 = 6.38e-01

CCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATAGGG  >  NC_002947/3162429‑3162711
                                                                                                                                           |                                                                                                                                               
ccccTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCt                                                                                                                                        >  2:182329/1‑149 (MQ=255)
 cccTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCtt                                                                                                                                       >  8:2173/1‑149 (MQ=255)
        cAAGCCCCATTGAGCTGCATTAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTgg                                                                                                                                <  7:118796/149‑1 (MQ=255)
                     gCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGTGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTTGCCTGAAAACAGCTTTTTTTGGCCCGATGACTTTa                                                                                                                   <  1:295344/149‑1 (MQ=255)
                      cTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTat                                                                                                                  <  5:111861/149‑1 (MQ=255)
                        gCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGAGTTTattc                                                                                                                >  2:238557/1‑149 (MQ=255)
                           tGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCAGGCACCCGATAGATATGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTattcatt                                                                                                             >  2:40195/1‑149 (MQ=255)
                                       ttCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCttt                                                                                                 >  1:153246/1‑149 (MQ=255)
                                        tCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCttt                                                                                                 <  8:228675/148‑1 (MQ=255)
                                          aGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCTGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGAAGACTTTATTCATTCAATGCGGCTTTTGt                                                                                              <  1:182329/149‑1 (MQ=255)
                                           gggCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTa                                                                                             <  6:158623/149‑1 (MQ=255)
                                             gggcgcTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                                                            <  7:249269/146‑1 (MQ=255)
                                                              acaAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGaat                                                                           >  3:210388/1‑148 (MQ=255)
                                                                 aaaaGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATc                                                                       <  5:151811/149‑1 (MQ=255)
                                                                      cACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTg                                                                  <  3:87851/149‑1 (MQ=255)
                                                                            tGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAAtt                                                             <  7:10850/148‑1 (MQ=255)
                                                                              gCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAgg                                                                                           <  7:116268/116‑1 (MQ=255)
                                                                              gCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAgg                                                                                           >  8:116268/1‑116 (MQ=255)
                                                                               cccGTTGCCTCGCCGCTGAGGGGGAGACTTTAGACTCCCGGCCACCTATAGATCTGGCCTCAAAACATCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGAc                                                         <  4:217676/147‑1 (MQ=255)
                                                                                  gTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGc                                                      <  7:2173/149‑1 (MQ=255)
                                                                                   taaagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAacac                                                                                    >  5:115606/6‑118 (MQ=255)
                                                                                   taaagtctcCCCCTAAGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAacac                                                                                    <  6:115606/113‑1 (MQ=255)
                                                                                        tctcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGtcat                                                 >  5:207852/1‑148 (MQ=255)
                                                                                        tctcCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTGTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTGGGGAAACACTCAATGAATAATCCGTTGGAAGTGGAGAGGGtcat                                                 >  4:55070/1‑148 (MQ=255)
                                                                                                   ggggAGACTTTAGACTACCGGCCCCCGAAAGATATGGCCTCAAAACAGATTTTTTAGCCCCGATCACTTTATTCATTCAATGCGGCTTTTGTATGGAAACACTCAATGAATAATCCGTTTGAATTGGACAGCGACATCAGCAGTAcaca                                     >  6:39099/1‑149 (MQ=255)
                                                                                                     ggAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                                                            <  1:205150/92‑1 (MQ=255)
                                                                                                     ggAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAg                                    >  6:173562/1‑148 (MQ=255)
                                                                                                     ggAGACTTTAGACTCACGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAg                                                                                            >  2:205150/1‑92 (MQ=255)
                                                                                                           tttAGACTCGCTTCCCCCGATAGTTGTGGGCTCAAAACAGCTTTTTTTGGCGCGTTGTCTTTATTCATTCAATGCGGCTTTTGTAGTGAAACTCTCTATGTATAATCCGTTGGATTTGGACATCGTCATCAGCAGTACACAGGATATTc                             <  5:39099/149‑1 (MQ=255)
                                                                                                           tttAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTc                             <  5:173562/149‑1 (MQ=255)
                                                                                                           tttAGACTCCCGGCCCCCGATAGATATGGACTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAAAACTCAATGAATAATCCGTTTGAATTGGACAGCGTCATCAGCAGTACACAGGTTAtta                             >  2:56530/1‑148 (MQ=255)
                                                                                                              aGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGa                                                          >  2:113490/1‑117 (MQ=255)
                                                                                                              aGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGa                                                          <  1:113490/117‑1 (MQ=255)
                                                                                                                  tCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGACACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTgcg                        >  8:225369/1‑147 (MQ=255)
                                                                                                                       gCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGATACACTCAATGAATAATCCGTTTGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTtgc                 >  1:102690/1‑149 (MQ=255)
                                                                                                                         ccacGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTGGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTtgctg               <  3:55070/146‑1 (MQ=255)
                                                                                                                                     tGGCCTTAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTTGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATa     <  1:100797/147‑1 (MQ=255)
                                                                                                                                     tGGCCTTAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATa     >  2:100797/1‑147 (MQ=255)
                                                                                                                                     tGGCCTTAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATa     >  2:104721/1‑147 (MQ=255)
                                                                                                                                       gCCTTAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATa     <  1:104721/145‑1 (MQ=255)
                                                                                                                                       gCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTTCGCAGTCGCTGGTGCTGGATAggg  <  1:56530/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
CCCCTCTCCAAGCCCCATTGAGCTGCATGAATCAAAGGTTTCAGGGCGCGCTGGGATTTCTCACAAAAAGCACTGATGGCACGTTGACTCTCCCCTTAGGGGGAGACTTTAGACTCCCGGCCCCCGATAGATCTGGCCTCAAAACAGCTTTTTTTGGCCCGATGACTTTATTCATTCAATGCGGCTTTTGTAGGGAAACACTCAATGAATAATCCGTTGGAATTGGACAGCGTCATCAGCAGTACACAGGATATTCTTGCGCAGTTGCTGGTGCTGGATAGGG  >  NC_002947/3162429‑3162711

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