Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,511 G→C 13.5% intergenic (+35/‑117) PP_3010 → / → PP_3011 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,5110GC13.5% 129.8 / 13.6 52intergenic (+35/‑117)PP_3010/PP_3011hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/18);  new base C (4/3);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

GTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCG  >  NC_002947/3402366‑3402648
                                                                                                                                                   |                                                                                                                                         
gTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGagag                                                                                                                                         <  8:191673/148‑1 (MQ=255)
gTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGc                                                                                                                                        >  3:309927/1‑149 (MQ=255)
     gTTGCGGGGCACGGGAAGCGAGCGCCACTTCCTCGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCACCCCCCACACCAATAGCCCGCTCTCGACGGCGGGc                                                                                                                                    >  1:39507/1‑149 (MQ=255)
     gTTGCGGGGCACGGGAAGCGAGCGCCACTTCAGCGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCAAACCAATAGCCCGCTCTCGACGGCGGGc                                                                                                                                    >  1:45261/1‑149 (MQ=255)
       tGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttt                                                                                                                                 >  7:361684/1‑149 (MQ=255)
         cGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCtattg                                                                                                                               <  2:45261/149‑2 (MQ=255)
         cGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCtattg                                                                                                                               <  2:39507/149‑2 (MQ=255)
         cGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCttttt                                                                                                                               <  3:349558/149‑1 (MQ=255)
            ggCACGGGAAGCGAGCTCCACTTCACGGCCTATCAACGTGATGCCTTGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgt                                                                                                                             <  5:232128/148‑1 (MQ=255)
            ggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCACAACACCAATAGCCAGCCGTCGAGAGCGGGCTATTtgt                                                                                                                             >  6:232128/1‑148 (MQ=255)
                 gggAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGcc                                                                                                                        >  6:108562/1‑148 (MQ=255)
                              cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                          <  3:212035/149‑1 (MQ=255)
                                        ggCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGtt                                                                                                   >  7:289741/1‑148 (MQ=255)
                                          ccTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGACCGGCAACGTTCa                                                                                                 >  3:261710/1‑148 (MQ=255)
                                              tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGGGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaa                                                                                            <  6:75965/149‑1 (MQ=255)
                                                aaCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAg                                                                                          >  6:341425/1‑149 (MQ=255)
                                                 aCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAgg                                                                                         >  1:51989/1‑149 (MQ=255)
                                                    tGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgt                                                                                                                           <  2:310428/112‑1 (MQ=255)
                                                    tGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgt                                                                                                                           >  1:310428/1‑112 (MQ=255)
                                                     gATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCgg                                                                                     >  6:92264/1‑149 (MQ=255)
                                                             gTACGACCATTGCTGGTTGCTTGATGCTGCACTGATGTGTTCAATAGCGTAATCATCCATCCCCCACACCAATAGGCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGcc                                                                             <  4:261710/149‑1 (MQ=255)
                                                                   ccATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                          >  5:152909/1‑114 (MQ=255)
                                                                   ccATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                          <  6:152909/114‑1 (MQ=255)
                                                                         cAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGa                                                                 <  2:51989/149‑1 (MQ=255)
                                                                            gTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTGATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGag                                                                >  8:57914/1‑147 (MQ=255)
                                                                                      gctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTg                                                                  >  4:110636/1‑135 (MQ=255)
                                                                                      gctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTg                                                                  <  3:110636/135‑1 (MQ=255)
                                                                                       ctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGtt                                                         >  6:61052/1‑143 (MQ=255)
                                                                                       ctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGtt                                                         <  5:61052/143‑1 (MQ=255)
                                                                                        tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTg                                                   >  6:206283/1‑148 (MQ=255)
                                                                                          cACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCAt                                                >  2:234182/1‑149 (MQ=255)
                                                                                          cACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCAt                                                >  5:229762/1‑149 (MQ=255)
                                                                                               aTGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGc                                            <  8:289741/148‑1 (MQ=255)
                                                                                                    ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                                      >  4:224187/1‑149 (MQ=255)
                                                                                                    ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                                      >  7:133838/1‑149 (MQ=255)
                                                                                                    ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                                      >  8:27318/1‑149 (MQ=255)
                                                                                                             gTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTaa                             <  6:229762/149‑1 (MQ=255)
                                                                                                             gTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTaa                             >  4:48246/1‑149 (MQ=255)
                                                                                                                 tcatcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATAtt                         >  5:146274/1‑149 (MQ=255)
                                                                                                                          ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAg                                                                               >  1:344947/1‑86 (MQ=255)
                                                                                                                          ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAg                                                                               <  2:344947/86‑1 (MQ=255)
                                                                                                                           cccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt               <  7:27318/149‑1 (MQ=255)
                                                                                                                            gcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg              >  2:357465/2‑149 (MQ=38)
                                                                                                                            ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg              >  3:97743/1‑149 (MQ=255)
                                                                                                                                acCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg          <  5:92264/149‑1 (MQ=255)
                                                                                                                                 caaaaaGCCCGCTCTCGACGGCGGGCTATTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                     <  6:151686/132‑1 (MQ=39)
                                                                                                                                 caaaaaGCCCGCTCTCGACGGCGGGCTATTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                     >  5:151686/6‑137 (MQ=39)
                                                                                                                                  cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgga        <  4:340759/149‑1 (MQ=255)
                                                                                                                                        cccGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg              <  7:70471/137‑1 (MQ=255)
                                                                                                                                        cccGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtg              >  8:70471/1‑137 (MQ=255)
                                                                                                                                          cgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCg  >  8:49683/1‑147 (MQ=255)
                                                                                                                                            ccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg          <  7:165822/137‑1 (MQ=255)
                                                                                                                                            ccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg          >  8:165822/1‑137 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCG  >  NC_002947/3402366‑3402648

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