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,2920CT95.5% 69.0 / NA 22intergenic (+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 (9/12);  minor base G (0/1);  total (9/13)
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.96e-01

GGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTG  >  NC_000913/2926046‑2926535
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
ggAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATa                                                                                                                                                                                                                                                 <  1:226463/251‑1 (MQ=255)
                                    gccagAAAAATTGAGCGTTGACGCCCGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGTGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGcgtg                                                                                                                                                                                                             <  2:75842/251‑4 (MQ=255)
                                                  gCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCGCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtg                                                                                                                                                                                               <  2:45028/251‑1 (MQ=255)
                                                                  taCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGa                                                                                                                                                                                 >  2:74401/1‑249 (MQ=255)
                                                                     ccGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTc                                                                                                                                                                                                                                      <  1:133467/193‑1 (MQ=255)
                                                                     ccGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTc                                                                                                                                                                                                                                      >  2:133467/1‑193 (MQ=255)
                                                                          gtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATg                                                                                                                                                                       <  2:236405/251‑1 (MQ=255)
                                                                                           aaCCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGTGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAGCGCACCGCca                                                                                                                                                       <  2:108695/250‑1 (MQ=255)
                                                                                                                              gACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCAc                                                                                                                   >  1:155688/1‑251 (MQ=255)
                                                                                                                                   caAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGcc                                                                                                              >  1:124427/1‑251 (MQ=255)
                                                                                                                                                   aTAAATTTTTTCTCAATTTGGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCAGCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTa                                                                                              <  2:115198/251‑1 (MQ=255)
                                                                                                                                                       aTTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGttt                                                                                          >  2:81745/1‑251 (MQ=255)
                                                                                                                                                                cAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATaa                                                                                 <  1:95155/251‑1 (MQ=255)
                                                                                                                                                                       gCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCaa                                                                           >  2:144470/1‑250 (MQ=255)
                                                                                                                                                                                      gCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGt                                                           <  1:101903/251‑1 (MQ=255)
                                                                                                                                                                                        aGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGttt                                                         <  1:179253/251‑1 (MQ=255)
                                                                                                                                                                                                     aGGTCATCAGGCAGGGCTTTTGTAATAAAAGGGAATGACGATATTCGACTCTTAAGGATTTTTTCTTTATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAACTGGAAGTGGATATTCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCGCGTTTTGAAAATTTCGat                                            <  2:160674/251‑1 (MQ=255)
                                                                                                                                                                                                           tcaGGCAGGGCTACTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCCACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAAc                                      >  2:159159/1‑251 (MQ=255)
                                                                                                                                                                                                                     cTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGcc                            >  1:129267/1‑251 (MQ=255)
                                                                                                                                                                                                                       aTTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTGGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGt                          <  1:9042/251‑1 (MQ=255)
                                                                                                                                                                                                                            aataaaaGGAAATGACGATATTCGTCTCATAAGGATTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACTTCTTGCGCCGTGAAc                      <  2:267762/246‑1 (MQ=255)
                                                                                                                                                                                                                                                ttCGTCTCATAAGGAGTTTTTCTTGATTACACTTCTTAGCCCGCTTGGCTCCGTGGATATGTTGTCGCAGCTGTAAGTGGCTATGCTTCAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGACCGCGGCGTAAAGATGGAACTg  >  2:160495/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GGAACGCATCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTG  >  NC_000913/2926046‑2926535

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