Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,631,564 A→T 7.8% intergenic (+1338/‑97) PP_4094 → / → PP_4096 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,631,5640AT 7.8% 151.6 / 2.7 51intergenic (+1338/‑97)PP_4094/PP_4096hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base T (2/2);  total (26/25)
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.63e-01

AGGGGTAGAGGGACATTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAG  >  NC_002947/4631423‑4631710
                                                                                                                                             |                                                                                                                                                  
aGGGGTAGAGGGACATTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCaa                                                                                                                                             <  8:206953/149‑1 (MQ=255)
    gTAGAGGGACATTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGAc                                                                                                                                         <  6:144181/149‑1 (MQ=255)
          ggACATTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAAcc                                                                                                                                   <  1:147633/149‑1 (MQ=255)
              aTTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTg                                                                                                                                 >  8:80298/1‑147 (MQ=255)
                    tgAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCaa                                                                                                                         <  6:332484/149‑1 (MQ=255)
                      aaTTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAg                                                                                                                        <  4:277345/148‑1 (MQ=255)
                        tttACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGc                                                                                                                       >  5:253068/1‑147 (MQ=255)
                             ttCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGcct                                                                                                                 >  2:250090/1‑146 (MQ=255)
                             ttCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGt                                                                                                                 >  2:123587/1‑148 (MQ=255)
                                cAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATa                                                                                                             <  3:384924/149‑1 (MQ=255)
                                cAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATa                                                                                                             <  3:154969/149‑1 (MQ=255)
                                         gcggcgCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATcc                                                                                                    <  8:275037/149‑1 (MQ=255)
                                              gcCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCaa                                                                                               >  7:274365/1‑149 (MQ=255)
                                                 aCTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTc                                                                                                        >  7:283482/1‑137 (MQ=255)
                                                 aCTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTc                                                                                            >  7:361322/1‑149 (MQ=255)
                                                   ttCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCgg                                                                                          >  7:243190/1‑149 (MQ=255)
                                                      tATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCAt                                                                                       >  3:194053/1‑149 (MQ=255)
                                                       aTCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGACCCTGAAGAGCAAGCCAGCGTCATAGCCTc                                                                                                        >  7:307401/1‑131 (MQ=255)
                                                       aTCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTc                                                                                                        <  8:307401/131‑1 (MQ=255)
                                                        tCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAg                                                                                                  <  6:354323/136‑1 (MQ=255)
                                                        tCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCa                                                                                     <  2:293655/149‑1 (MQ=255)
                                                        tCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCa                                                                                     <  8:361322/149‑1 (MQ=255)
                                                        tCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAACAGCAAGCCAGCGTCATAGCCTCATCCAg                                                                                                  >  5:354323/1‑136 (MQ=255)
                                                           aaaGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTCGCACTGGCAAGGAAGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGCCATCAGTg                                                                                  >  2:196954/1‑149 (MQ=255)
                                                                    aaCCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTc                                                                             >  3:150876/1‑145 (MQ=255)
                                                                    aaCCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTc                                                                             <  4:150876/145‑1 (MQ=255)
                                                                     aCCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATgc                                                                        <  6:212173/149‑1 (MQ=255)
                                                                               gTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCg                                                               >  2:196333/1‑148 (MQ=255)
                                                                               gTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCg                                                               >  2:89236/1‑148 (MQ=255)
                                                                                   aTCGTTGAAACAAAGGAGTAGCCGTTGTTCTGCACCGTGCCACCCGTTCATATCTGGCTCTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGc                                                                                                           <  1:242344/100‑1 (MQ=255)
                                                                                   aTCGTTGAAACAAAGGAGTAGCCGTTGTTCTGCACCGTGCCACCCGTTCATATCTGGCTCTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGc                                                                                                           >  2:236353/1‑100 (MQ=255)
                                                                                   aTCGTTGAAACAAAGGAGTAGCCGTTGTTCTGCACCGTGCCACCCGTTCATATCTGGCTCTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGc                                                                                                           <  1:236353/100‑1 (MQ=255)
                                                                                   aTCGTTGAAACAAAGGAGTAGCCGTTGTTCTGCACCGTGCCACCCGTTCATATCTGGCTCTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGc                                                                                                           >  2:242344/1‑100 (MQ=255)
                                                                                      gTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACt                                                        >  8:378308/1‑148 (MQ=255)
                                                                                       ttGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCa                                                      <  4:194053/149‑1 (MQ=255)
                                                                                         gAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAg                                                     <  4:295235/148‑1 (MQ=255)
                                                                                          aaaCAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGa                                                   >  4:163437/1‑149 (MQ=255)
                                                                                           aaCAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGAt                                                  >  5:149349/1‑149 (MQ=255)
                                                                                           aaCAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGAt                                                  >  7:119708/1‑149 (MQ=255)
                                                                                              aaaGGAGTAGCCGTTGTTCCACTCCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATgg                                                >  1:234033/1‑148 (MQ=255)
                                                                                                       gCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGgctc                                       >  6:266503/1‑146 (MQ=255)
                                                                                                              ttCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGttt                               <  7:80298/149‑1 (MQ=255)
                                                                                                               tCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGttt                               <  8:243190/148‑1 (MQ=255)
                                                                                                                  gCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCAt                                                            >  4:393806/1‑116 (MQ=255)
                                                                                                                  gCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCAt                                                            <  3:393806/116‑1 (MQ=255)
                                                                                                                  gCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAg                           <  1:196333/149‑1 (MQ=255)
                                                                                                                   cACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGc                          >  5:306180/1‑149 (MQ=255)
                                                                                                                      cGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGt                       <  6:149349/149‑1 (MQ=255)
                                                                                                                              ccGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTg               >  4:70382/1‑149 (MQ=255)
                                                                                                                              ccGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTg               <  2:234033/149‑1 (MQ=255)
                                                                                                                                    atatCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTGCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTTGAAtt          <  1:250090/148‑1 (MQ=255)
                                                                                                                                            cactggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTGCGATCTGGAATTGCCAGGAg  >  2:100071/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
AGGGGTAGAGGGACATTGTGTGAATTTACTTCCAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAG  >  NC_002947/4631423‑4631710

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