Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,584 A→G 9.5% intergenic (‑431/+15) PP_0022 ← / ← PP_0023 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294727,5840AG 9.5% 113.7 / 3.1 42intergenic (‑431/+15)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGC  >  NC_002947/27449‑27729
                                                                                                                                       |                                                                                                                                                 
gCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttttt                                                                                                                                      <  1:132255/149‑1 (MQ=255)
  aCTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                    >  1:168908/1‑149 (MQ=255)
    tGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCtt                                                                                                                                   >  8:163759/1‑148 (MQ=255)
          gATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTcccc                                                                                                                            >  4:66889/1‑149 (MQ=255)
            tAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCt                                                                                                                          >  2:137490/1‑149 (MQ=255)
            tAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCt                                                                                                                          >  3:165726/1‑149 (MQ=255)
               ccGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCTCGGGAGCTGCttttt                                                                                                                                       <  6:46026/133‑1 (MQ=255)
               ccGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCTCGGGAGCTGCttttt                                                                                                                                       >  5:46026/1‑133 (MQ=255)
                           tAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTgg                                                                                                            <  7:163759/148‑1 (MQ=255)
                                   gAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGt                                                                                                  >  5:57901/1‑149 (MQ=255)
                                      tCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCTCGGGAGCTGCttttttt                                                                                                                                     <  7:172045/112‑2 (MQ=255)
                                      tCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCTCGGGAGCGGCttttttt                                                                                                                                     >  8:172045/1‑111 (MQ=255)
                                                 aagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCtt                                                                                     <  8:177580/149‑1 (MQ=255)
                                                            tCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTACCGAGAGCTGCTTTTTTATTAATCCCCCTTTTTGTGCTTGTGGGGTTTAAGGTcc                                                                                                >  6:185645/1‑127 (MQ=255)
                                                              acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTTTGCTTGTGGGGTTTCAGGTCCGGATGGt                                                                                         >  2:197801/1‑132 (MQ=255)
                                                              acgacgCGGTCCGATTGATCAGCGAGACGGACGTGAGGGTTTGCACACAGCCCCAAAAAAAACGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                        <  5:263447/149‑1 (MQ=255)
                                                                       tCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCa                                                                           <  6:234415/137‑1 (MQ=255)
                                                                       tCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCa                                                                           >  5:234415/1‑137 (MQ=255)
                                                                          gATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc                                                            <  1:137490/149‑1 (MQ=255)
                                                                                       gaCGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAg                                               <  4:165726/149‑1 (MQ=255)
                                                                                         gcgTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGGCTCGCTGCCTTTTGTCCGGAGTg                                             <  2:19788/147‑1 (MQ=255)
                                                                                              gTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATc                                        <  2:179131/149‑1 (MQ=255)
                                                                                               tGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                  <  8:155449/138‑1 (MQ=255)
                                                                                               tGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                  >  7:155449/1‑138 (MQ=255)
                                                                                               tGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCt                                       >  1:210448/1‑149 (MQ=255)
                                                                                                   ggTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGa                                                >  8:23100/1‑136 (MQ=255)
                                                                                                   ggTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGa                                                <  7:23100/136‑1 (MQ=255)
                                                                                                    gTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCAtt                                  <  2:210448/149‑1 (MQ=255)
                                                                                                    gTTTGCACACAGCCCC‑AAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTg                                 <  6:57901/149‑1 (MQ=255)
                                                                                                                 cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGttt                                                                                                       <  2:254381/67‑1 (MQ=255)
                                                                                                                 cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGttt                                                                                                       >  1:254381/1‑67 (MQ=255)
                                                                                                                  ccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCt                                                                                >  8:123868/1‑89 (MQ=255)
                                                                                                                   cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGt                                                                      >  2:79622/1‑98 (MQ=255)
                                                                                                                   cAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGt                                                                      <  1:79622/98‑1 (MQ=255)
                                                                                                                    aaaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTTAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCttt                        >  3:239172/1‑143 (MQ=255)
                                                                                                                    aaaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTTAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCttt                        <  4:239172/143‑1 (MQ=255)
                                                                                                                    aaaaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                  >  4:235945/1‑149 (MQ=255)
                                                                                                                            agcagcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGAt          <  3:66889/149‑1 (MQ=255)
                                                                                                                                  tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCAt    >  6:43673/1‑149 (MQ=255)
                                                                                                                                   cccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCtt                                                         <  7:164509/95‑1 (MQ=255)
                                                                                                                                   cccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCtt                                                         >  8:164509/1‑95 (MQ=255)
                                                                                                                                     cGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGc  >  6:222800/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
GCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGC  >  NC_002947/27449‑27729

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: ‑