Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,458,124 T→G 8.5% intergenic (+17/+23) PP_1276 → / ← algA hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,458,1240TG 8.5% 151.7 / 8.1 59intergenic (+17/+23)PP_1276/algAhypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base G (2/3);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

TCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGAT  >  NC_002947/1457980‑1458269
                                                                                                                                                |                                                                                                                                                 
tCGAAATGACTTCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGcc                                                                                                                                              >  1:1583634/1‑150 (MQ=255)
  gAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCa                                                                                                                                            >  2:1344537/1‑150 (MQ=255)
         cTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                     <  1:1267142/150‑1 (MQ=255)
            ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTc                                                                                                                                  >  1:962628/1‑150 (MQ=255)
              acaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGCCAg                                                                                                                                <  2:1258797/150‑1 (MQ=255)
               caAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCTGACTCGCCCCTGGGCTGACCAGGGGCGAgcctgc                                                                                                                               <  1:948224/150‑6 (MQ=38)
               caAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCACGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGc                                                                                                                               <  1:1428972/150‑1 (MQ=255)
                aaGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCa                                                                                                                              >  2:773410/1‑150 (MQ=255)
                     tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTAc                                                                                                                         >  2:1445449/1‑150 (MQ=255)
                     tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTAc                                                                                                                         >  1:1375082/1‑150 (MQ=255)
                              aaTCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCACATTCGACATCTACGTCATAGCCCCGCGTGTTGAATATGGAGGCTGGCCCCTGGTCAGCCCATGGGCGACTAAGCATCTACAGCTCGATg                                                                                                                >  2:385574/1‑150 (MQ=255)
                                   ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGt                                                                                                              <  1:1175025/147‑1 (MQ=255)
                                   ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGt                                                                                                              >  2:1175025/1‑147 (MQ=255)
                                         tttGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGATCGCCCCTGGTCAGCCCAGGGGCGAGTCAGTATTTACTGCGCGATGGTCTTCACCGa                                                                                                     <  1:1066877/150‑1 (MQ=255)
                                             aaGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAggg                                                                                                                                         <  1:217966/110‑1 (MQ=255)
                                             aaGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAggg                                                                                                                                         >  2:217966/1‑110 (MQ=255)
                                                gggAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAAcgccg                                                                                              <  2:454147/150‑1 (MQ=255)
                                                 ggaggaTGATTTGCAGCGAACTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAAcgccgc                                                                                             <  2:1106431/150‑1 (MQ=255)
                                                            tGCAGCGACCTTTCAGGGTGGTGGGTCGCCCATTCAACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCCGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGGCTTCACCGAAACGCCGCGCTCGATgggg                                                                                  <  1:385574/150‑1 (MQ=255)
                                                                gCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTgg                                                                              >  2:53296/1‑150 (MQ=255)
                                                                 cGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGa                                                                             <  1:892269/150‑1 (MQ=255)
                                                                     cTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATg                                                                                                                >  2:869268/1‑111 (MQ=255)
                                                                     cTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATg                                                                                                                <  1:869268/111‑1 (MQ=255)
                                                                     cTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTg                                                                         >  1:1017215/1‑150 (MQ=255)
                                                                      tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                     >  2:1558756/1‑89 (MQ=255)
                                                                      tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                     <  1:1558756/89‑1 (MQ=255)
                                                                      tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTgc                                                                        >  2:750281/1‑150 (MQ=255)
                                                                      tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTgc                                                                        >  2:821148/1‑150 (MQ=255)
                                                                      tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTgc                                                                        >  1:1174715/1‑150 (MQ=255)
                                                                           ggttggtggGTGGCCCATTCGACATCTACGTCCCAGCCAACAGCTTTTAATAAGCATGCTCGCCCCTGGTCAATATAAGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCa                                                                   <  1:773410/147‑1 (MQ=255)
                                                                                  gggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACAt                                                            >  2:872116/1‑150 (MQ=255)
                                                                                   ggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGAt                                                                                      >  1:730712/1‑123 (MQ=255)
                                                                                   ggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGAt                                                                                      <  2:730712/123‑1 (MQ=255)
                                                                                        cccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGTGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACAACCTCGa                                                      <  1:1344537/150‑1 (MQ=255)
                                                                                         ccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                     >  1:340124/1‑150 (MQ=255)
                                                                                            ttCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGAc                                                                            >  1:714683/1‑124 (MQ=255)
                                                                                            ttCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGAc                                                                            <  2:714683/124‑1 (MQ=255)
                                                                                                  aTCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGc                                                                                                                                               >  2:851150/1‑51 (MQ=255)
                                                                                                  aTCTACCGCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGc                                                                                                                                               <  1:851150/51‑1 (MQ=255)
                                                                                                     tACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTc                                         >  2:72123/1‑150 (MQ=255)
                                                                                                             aGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGc                                 <  1:408142/150‑1 (MQ=255)
                                                                                                              gCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                     <  2:1094436/129‑1 (MQ=255)
                                                                                                              gCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                     >  1:1094436/1‑129 (MQ=255)
                                                                                                               ccccGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGccc                               <  1:1445449/150‑1 (MQ=255)
                                                                                                                   gAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGt                           <  1:53296/150‑1 (MQ=255)
                                                                                                                          gAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTccc                    >  1:1499292/1‑150 (MQ=255)
                                                                                                                                gCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGTGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTg              >  2:375125/1‑150 (MQ=255)
                                                                                                                                 cAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGc             <  2:962628/150‑1 (MQ=255)
                                                                                                                                   gacTCGCCCCTGGGCTGACCAGGGGCGAGCCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGccc                                                                    >  2:230473/3‑93 (MQ=37)
                                                                                                                                   gacTCGCCCCTGGGCTGACCAGGGGCGAGCCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGccc                                                                    <  1:230473/91‑1 (MQ=37)
                                                                                                                                      tCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTc        >  2:1507476/1‑150 (MQ=255)
                                                                                                                                      tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCg                                              <  2:774309/112‑1 (MQ=255)
                                                                                                                                      tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCAATACACATCCTCGAAGCGCTCg                                              >  1:774309/1‑112 (MQ=255)
                                                                                                                                        gCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAACCGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCGCGATAGCATCCTCGCCCAGGTATCTCCCGGACTGCACCTCga      >  2:1320151/1‑150 (MQ=255)
                                                                                                                                         ccccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat     <  1:72123/150‑1 (MQ=255)
                                                                                                                                          cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGcc                                <  2:46983/122‑1 (MQ=255)
                                                                                                                                          cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGcc                                >  1:46983/1‑122 (MQ=255)
                                                                                                                                          cccTGGTCAGCCCAGGGGCGAGTCAGCATGTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgatg    <  1:821148/150‑1 (MQ=255)
                                                                                                                                            cTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgatgat  >  2:611892/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGAT  >  NC_002947/1457980‑1458269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: