Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 699,722 A→G 51.8% intergenic (+115/‑129) PP_t06 → / → PP_23SD tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,7220AG51.8% ‑6.5 / 69.8 54intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (12/14);  new base G (15/13);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

CGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT  >  NC_002947/699582‑699803
                                                                                                                                              |                                                                                   
cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAacac                                                                             <  1:18350/149‑1 (MQ=37)
 ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACAcc                                                                            >  4:86037/1‑149 (MQ=37)
   ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                                                                          >  8:96913/1‑149 (MQ=37)
   ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg                                                                          >  1:106193/1‑149 (MQ=37)
          ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCAct                                                                 >  4:337158/1‑149 (MQ=37)
          ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCAc                                                                  <  3:306060/148‑1 (MQ=37)
           cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCActg                                                                <  6:348878/149‑1 (MQ=37)
            gCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAac                                                                               >  4:294088/1‑135 (MQ=37)
            gCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAac                                                                               <  3:294088/135‑1 (MQ=37)
                  cTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCtggtgg                                                         <  4:77162/149‑1 (MQ=37)
                   tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                          <  6:195876/133‑1 (MQ=37)
                   tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                          >  5:195876/1‑133 (MQ=37)
                   tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGa                                                        >  4:233290/1‑149 (MQ=37)
                   tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGa                                                        >  5:159697/1‑149 (MQ=37)
                     caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGAt                                                       >  6:375926/1‑148 (MQ=37)
                         acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGc                                                  >  1:205546/1‑149 (MQ=37)
                          cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                  <  2:368544/118‑1 (MQ=39)
                          cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                  >  1:368544/1‑118 (MQ=39)
                          cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGc                                                  >  1:58810/1‑148 (MQ=37)
                             tCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTaa                                               <  7:110308/148‑1 (MQ=37)
                                ccGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGt                                            >  8:58467/1‑148 (MQ=37)
                                   tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCtggt                                                           >  6:240593/1‑130 (MQ=38)
                                   tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCtggt                                                           <  5:240593/130‑1 (MQ=38)
                                       cAGTTACTTTTCATAAATTATAAATTAAAATTCATCTCTAATTTTTATTTTTGGCTTTTTTCATATCGTTCTTTAATATATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                                     <  1:117146/149‑1 (MQ=11)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  5:345017/149‑1 (MQ=37)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGAATCAGGCTAAGGTAAAATTtgt                                   <  7:58467/149‑1 (MQ=255)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAAttttt                                   >  3:313823/1‑147 (MQ=37)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  1:148087/149‑1 (MQ=37)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  3:57299/149‑1 (MQ=37)
                                        aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                   <  8:264864/149‑1 (MQ=37)
                                          ttACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgtg                                  >  3:341610/1‑148 (MQ=25)
                                           tACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                 <  4:341610/148‑1 (MQ=25)
                                            aCTTGTCAGACCTTATAAAATAACATTCAGCTCGACTGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                <  8:306475/148‑1 (MQ=11)
                                            aCTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                <  4:350726/148‑1 (MQ=21)
                                            aCTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                >  3:350726/1‑148 (MQ=21)
                                              ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAgg                                                   <  5:88400/127‑1 (MQ=35)
                                              ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAgg                                                   >  6:88400/1‑127 (MQ=35)
                                                gTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                            <  2:205546/148‑1 (MQ=21)
                                                 tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg                           <  2:316570/148‑1 (MQ=21)
                                                  cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc                          <  2:58810/148‑1 (MQ=34)
                                                   aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTAAGGAAATGATGGTGGAAAAGGGTAAGGTAAAATTTTTTAGTTATga                          >  1:126690/1‑148 (MQ=11)
                                                     aaCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                       >  8:47913/1‑148 (MQ=21)
                                                     aaCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGa                      <  5:375926/149‑1 (MQ=34)
                                                       cTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                    <  3:210972/149‑1 (MQ=21)
                                                          aGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                  <  7:96913/148‑1 (MQ=34)
                                                            aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTaaaa                                        >  7:214490/1‑124 (MQ=25)
                                                            aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTaaaa                                        <  8:214490/124‑1 (MQ=25)
                                                               tcaaCATTCGAGCTGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg            >  6:316329/3‑149 (MQ=2)
                                                               tcaaCATTCGAGCTGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg            >  1:160404/3‑149 (MQ=2)
                                                               tGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc             <  7:20398/148‑1 (MQ=21)
                                                               tGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc             >  8:20398/1‑148 (MQ=21)
                                                                      tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAt     >  5:352964/1‑149 (MQ=18)
                                                                          gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt  >  8:36640/1‑149 (MQ=14)
                                                                          gCTCGAATGTTGATTTCTGGATTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt  >  4:349997/1‑149 (MQ=14)
                                                                                                                                              |                                                                                   
CGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGT‑‑TTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTT  >  NC_002947/699582‑699803

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: