Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,255,974 G→T 7.7% intergenic (+34/‑670) asd → / → PP_1990 aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,255,9740GT 7.7% 152.4 / 4.1 52intergenic (+34/‑670)asd/PP_1990aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base T (2/2);  total (27/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 = 2.12e-01

CGCAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAG  >  NC_002947/2255827‑2256120
                                                                                                                                                   |                                                                                                                                                   
cgcAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGt                                                                                                                                                    >  3:136119/1‑149 (MQ=255)
       tGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTTGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTg                                                                                                                                             >  1:186833/1‑149 (MQ=255)
         aaCATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCg                                                                                                                                           >  6:120476/1‑149 (MQ=255)
          aCATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCgg                                                                                                                                          >  5:171708/1‑149 (MQ=255)
               gggTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGTCACCAGTGCGGTAttt                                                                                                                                     <  3:202412/149‑1 (MQ=255)
               gggTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTAttt                                                                                                                                     >  5:171117/1‑149 (MQ=255)
                   cGCAGTACCTGGGCGCGGTCACCGGGGGGGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTAttttttt                                                                                                                                 <  6:296446/149‑1 (MQ=255)
                    gCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTg                                                                                                                                                   >  6:180033/1‑130 (MQ=255)
                    gCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTg                                                                                                                                                   <  5:180033/130‑1 (MQ=255)
                    gCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTa                                                                                                                                <  6:171117/149‑1 (MQ=255)
                      attACCTGGGCGCGTTAACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTATGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGATAGTGCGGTATTTTTTTATg                                                                                                                              >  2:204662/3‑149 (MQ=255)
                            tGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGa                                                                                                                         >  6:36132/1‑148 (MQ=255)
                               gcgcGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAg                                                                                                                      >  3:361631/1‑148 (MQ=255)
                                 gcgTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTAt                                                                                                                   <  3:156961/149‑1 (MQ=255)
                                 gcgTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTAt                                                                                                                   <  5:120476/149‑1 (MQ=255)
                                    ttCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCa                                                                                                                <  6:31502/149‑1 (MQ=255)
                                       aCCGTGGGTGACCAGCTGCTGTGGGGCACCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTAGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCAttt                                                                                                             <  1:197230/149‑1 (MQ=255)
                                              gTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTcgc                                                                                                      <  3:157766/149‑1 (MQ=255)
                                                gACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCgg                                                                                                    >  8:170349/1‑149 (MQ=255)
                                                gACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCgg                                                                                                    <  3:210945/149‑1 (MQ=255)
                                                 aCCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCgg                                                                                                    <  3:284414/148‑1 (MQ=255)
                                                   cAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGcc                                                                                                 >  1:358188/1‑149 (MQ=255)
                                                      ctgctgTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTaa                                                                                              <  1:120062/149‑1 (MQ=255)
                                                          tgtgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGc                                                                                          >  7:87266/1‑149 (MQ=255)
                                                            tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGcc                                                                                         <  4:361631/148‑1 (MQ=255)
                                                            tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGc                                                                                          <  2:358188/147‑1 (MQ=255)
                                                            tgGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCg                                                                                        <  5:36132/149‑1 (MQ=255)
                                                             ggggCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCg                                                                                        >  2:181036/1‑148 (MQ=255)
                                                                 cgccgcCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATaa                                                                                   <  4:136119/149‑1 (MQ=255)
                                                                                  gccgcATGTTGCGGATTCTGCTGTATCGTTGATCGTTCAGGCTTCACCTGAAAAATACCTCACTGGTGACAGTTCGGTATTTTTTTATGCCTTAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTaa                                                                  <  1:204662/149‑1 (MQ=255)
                                                                                   ccgcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGTCTAATCGCCGGATAATTCATGCCGCTGTTTAAt                                                                 <  1:239148/149‑1 (MQ=255)
                                                                                     gcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATgg                                                               <  4:53791/149‑1 (MQ=255)
                                                                                     gcATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATgg                                                               >  5:238283/1‑149 (MQ=255)
                                                                                      cATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGc                                                              <  6:171708/149‑1 (MQ=255)
                                                                                               gATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCt                                                     >  7:230100/1‑149 (MQ=255)
                                                                                                        ggAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGa                                              >  5:13472/1‑147 (MQ=255)
                                                                                                              ttGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccg                                     >  8:87826/1‑149 (MQ=255)
                                                                                                                 aTCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgtt                                   <  1:181036/148‑1 (MQ=255)
                                                                                                                  tCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGACGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgc                                 >  5:254684/1‑149 (MQ=255)
                                                                                                                    gTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgccg                               >  6:335497/1‑149 (MQ=255)
                                                                                                                    gTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgcc                                <  8:87266/148‑1 (MQ=255)
                                                                                                                     ttCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgccgt                              >  5:179756/1‑149 (MQ=255)
                                                                                                                     ttCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAgccgttgccgt                              >  5:83478/1‑149 (MQ=255)
                                                                                                                                ccTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGATGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCAGTTGCGGttgctt                   >  2:281676/1‑149 (MQ=255)
                                                                                                                                    aaaaaTACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctttt                 <  8:230100/147‑1 (MQ=255)
                                                                                                                                        aTACCGCACTGTCACCAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgcttt            <  2:115379/148‑1 (MQ=255)
                                                                                                                                        aTACCGCACTGTCACCAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgcttttgcttt            >  1:115379/1‑148 (MQ=255)
                                                                                                                                        aTACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCCGTTGCCGTTGCCGttgcttttgcttt            <  1:4974/149‑1 (MQ=255)
                                                                                                                                         tACCGCACTGTCACCAGTGCGGTATTTTTTTATTCCTGAAAGTATGCATATGGGTCGCGGGCCTAATCGCCGGAAAATTCATGCCGCTGTTTAATGTCGAAATGTCTGCTCTGACAGAGAAGTTTCCGTTGCAGttgcttttgcttttg          >  8:318543/1‑149 (MQ=255)
                                                                                                                                             gCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTc      >  3:209812/1‑149 (MQ=255)
                                                                                                                                              cACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTc      >  7:85459/1‑148 (MQ=255)
                                                                                                                                                 tGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAg  >  8:157354/1‑149 (MQ=255)
                                                                                                                                                   gTGACAGTGCGGTATTTTTTTATGCCTGTAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGttgctt                   <  3:2093/130‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCAAGCTGAACATGGGGTCGCAGTACCTGGGCGCGTTCACCGTGGGTGACCAGCTGCTGTGGGGCGCCGCCGAACCGCTGCGCCGCATGTTGCGGATTCTGCTGGAGCGTTGATCGTTCAGGCTTCACCTGAAAAATACCGCACTGGTGACAGTGCGGTATTTTTTTATGCCTGAAAGTATGCATTTGGGTCGCGGGCCTAATCGCCGGATAATTCATGCCGCTGTTTAATGGCGAAATGTCTGCTCTGACAGAGCCGTTGCCGTTGCCGTTGCTTTTGCTTTTGCTTCTAAG  >  NC_002947/2255827‑2256120

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: