Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,292 C→T intergenic (+96/‑16) queF → / → ppnN 7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,2920CT97.0% 106.7 / NA 33intergenic (+96/‑16)queF/ppnN7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/16);  minor base A (1/0);  total (17/16)
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.83e-01

AACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTG  >  NC_000913/2926048‑2926529
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
aaCGCATCTTTAATGACCTGTTACGCTGCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAAGTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAg                                                                                                                                                                                                                                         <  1:23726/251‑1 (MQ=255)
              gACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGAt                                                                                                                                                                                                                           <  1:101679/251‑1 (MQ=255)
                           tCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAg                                                                                                                                                                                                               <  1:229566/250‑1 (MQ=255)
                              gccagccaGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGccc                                                                                                                                                                                                            <  1:172648/250‑1 (MQ=255)
                                    cagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGc                                                                                                                                                                                                     >  1:173262/1‑251 (MQ=255)
                                           aTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATgg                                                                                                                                                                                              <  2:249472/251‑1 (MQ=255)
                                           aTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGGGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATgg                                                                                                                                                                                              >  2:148666/1‑251 (MQ=255)
                                                  gTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                       >  1:230550/1‑251 (MQ=255)
                                                  gTTTACGCACGTTATACCCGTCGAGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGGTTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCCATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGGCATCAGGCAGGGCTATTGGAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                       >  2:28700/1‑251 (MQ=255)
                                                        gCACGATATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACCAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGCAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCAGCTTGGCTCCATGGATAtgtttccgc                                                                                                                                                                                  >  2:102008/1‑245 (MQ=255)
                                                           cGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAGATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCt                                                                                                                                                                              <  1:220639/251‑1 (MQ=255)
                                                           cGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCt                                                                                                                                                                              <  1:271630/251‑1 (MQ=255)
                                                              tataCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTgg                                                                                                                                                                            <  1:79182/250‑1 (MQ=255)
                                                                 aCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                                                                                                                                                        <  1:7609/251‑1 (MQ=255)
                                                                                  ggACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGc                                                                                                                                                        >  1:154133/1‑250 (MQ=255)
                                                                                                                              ccACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTTACGCAGAAGGTTGTGAAAGGTCATCATGCAGGGCTATTGTAGTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACtt                                                                                                           >  2:42830/1‑251 (MQ=255)
                                                                                                                                            cGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTACCGCCGAAGGCAGTGAGAGGTCATCAGGCACGGCTATTGGAATCCAAGGGAATGGCGATATTAGTCTTATTAGGAGGTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCGTGGATATGTTGTCGCAGCTGGAAGTGGAGATGCTTAAACGCACCGCCAGCAACGACCTCTATCAACTGTTTGGCAACTGTTCACTTGACGTACTGAACTc                                                                                             >  1:138516/1‑251 (MQ=255)
                                                                                                                                                       tttttCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTg                                                                                   >  1:175715/1‑250 (MQ=255)
                                                                                                                                                             tcAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaaa                                                                                                                                                           >  2:52578/1‑172 (MQ=255)
                                                                                                                                                             tcAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTaaa                                                                                                                                                           <  1:52578/172‑1 (MQ=255)
                                                                                                                                                             tcAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGAAATCAGGCAGCGCTCTTGAAATCAAAGGGAATGACGATCTTCCTCTAATAAGGAGTTTTTCTTGATTACACATATTAGCCCGATTGGCTCCATGGATATGTTGTCGCAACTGGAAGTGGATATGCTTAAACGCACAGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGCTTGACCGAt                                                                             >  1:298838/1‑250 (MQ=255)
                                                                                                                                                                    gtcttGCTGGATTCACGAAGAGGGTTGTGAAAGGTCATCAGGCAGGTCTATTGTAATCAAAGGGAATGACTGTATTCGTCTCAAAAGAAGTTTTTCTTGATTACACATATTAGCCCGCTTGTCTCCATGGATATGTTGTTGCAGCTGGAAGTCGATTTGCTTAATCGCACCGTCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCTTACTGAACTCCGGTAGTTTGACCGATAACAGCaa                                                                     <  1:240334/247‑1 (MQ=255)
                                                                                                                                                                        ggCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGaa                                                                 <  2:231721/250‑1 (MQ=255)
                                                                                                                                                                            ggATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGc                                                             >  1:173592/1‑251 (MQ=255)
                                                                                                                                                                                  aCGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGGCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGtctc                                                       <  1:256117/251‑1 (MQ=255)
                                                                                                                                                                                    gCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGt                                                     >  1:76648/1‑251 (MQ=255)
                                                                                                                                                                                        aagttggtgaacgttcttcaGGCAGGGCTATTGTAATCAAAGGAAATGACTATATTCTTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTTTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACATCAAAGAATTGCTGGCTCTTTTTg                                                 <  2:312954/234‑1 (MQ=255)
                                                                                                                                                                                          ggTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGa                                                >  2:254916/1‑250 (MQ=255)
                                                                                                                                                                                                gAAAGGTCACCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTATTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTc                                         <  1:299489/251‑1 (MQ=255)
                                                                                                                                                                                                  aagttcctccttcaggGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCTACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCTTTCGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCg                                        >  1:148297/13‑250 (MQ=255)
                                                                                                                                                                                                                           aTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTcc                                                                                            >  2:125375/1‑173 (MQ=255)
                                                                                                                                                                                                                           aTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTcc                                                                                            <  1:125375/173‑1 (MQ=255)
                                                                                                                                                                                                                                       gACGATATTCGTCACATAATGCGTTTTTATTGATTACACATATTAGCCCGCTAGCCTCCATGGATACGTTATCGCAGCTGCAAGTGGATATGAATAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACACCGGTAGTTTGACAGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGCTATTAACGACTAGCGCCGTGAACGCGGCGTAAGGCTg  >  2:303050/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
AACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTG  >  NC_000913/2926048‑2926529

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