Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,322,605 C→T intergenic (‑183/+731) phnF ← / ← phnD putative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,322,6050CT98.4% 187.3 / NA 62intergenic (‑183/+731)phnF/phnDputative transcriptional regulator PhnF/phosphonate/phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (39/21);  minor base G (1/0);  total (40/22)
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.66e-01

CAAGTTTTGCGGCTATCTCTTGATAGCGTGTTGGGTAGCTGGTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGC  >  NC_000913/4322359‑4322846
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cAAGTTTGGCGGCTATCTCTTGATAGCGTTTTGGGTAGCTGGTCGGATGTGTAGACAAGGGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCg                                                                                                                                                                                                                                               <  2:121220/251‑1 (MQ=255)
  aGTTTGGCGGCTATCTCTTGATAGCGTGTTGGGTAGCTGGTCGGATGTGTAGACAAGTGCAAAGATATCAATGCCTGGCTTATCAGATAAAGGGCTTACGCAACCGTTGAGAGAGGAAACCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGTTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAg                                                                                                                                                                                                                                             <  1:255509/251‑1 (MQ=255)
     tttGCGGCTATCTCTTGGTAACGTGTTGTGTAGCTGGTCGGATGTGTAGTCAAGTGCATAGATATCAATGCCTCGCTTATCAGATAACGGGCTTACGCATCCGTTTAGAGAGGATATCCAGCAGGCTGACCGTGACGATGACAAGCACCATC‑AGGGAGCAGGTTTGTTGGAAATGGAAACCACGAATCGCTTACCACAGGGTGACGTCGATCCCGCCTGCGCCTACCATGCCGACGACGGTTGCCGAGctg                                                                                                                                                                                                                                          >  2:17725/1‑249 (MQ=255)
          ggCTATCTCTTGATAGCGTGGTGGGTAGCTGGTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGtt                                                                                                                                                                                                                                     <  1:279419/251‑1 (MQ=255)
                                gggTAGCTGGTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGTTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAg                                                                                                                                                                                                               <  1:177778/251‑1 (MQ=255)
                                    aGCTGGTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTCTCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCAGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAGGCGATAGAGGGAGTAgg                                                                                                                                                                                                           >  1:32707/1‑251 (MQ=255)
                                        ggTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGAGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGgagtt                                                                                                                                                                                                       >  2:298330/1‑249 (MQ=255)
                                        ggTCGGATGGTTAGCCAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGGGCTGCCGCAACCGTTGAGAGAGGAAAGCCAGCAGGTTGACGGTGACGATTATAAGCACCAGC‑AGGGCGCAGGTTTGGTGGAATGGGAATCCGGGATTCGCTTCCCACAGGGGGCCGCGGATCCCGCCTTCGCCGACCATGCTGACGACGGTTGCCGAGCGGACGTTGGAGTCGAAGCGATAGAGTGAGTAGGAGAt                                                                                                                                                                                                       <  2:52763/251‑1 (MQ=255)
                                           cggttgtgTAGACAAGTGCTTAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCTCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATcag                                                                                                                                                                                                    <  1:191338/247‑1 (MQ=255)
                                           cGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGCTTCGAAGCGCTAGAGGGAGTAGGAGATc                                                                                                                                                                                                      >  2:93921/1‑249 (MQ=255)
                                               tgtgTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGt                                                                                                                                                                                                >  1:198569/1‑251 (MQ=255)
                                                    aGACAAGGGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGCCGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACATGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGATGCGTTAGTGGGCGTAGGAGATCAGCAGTGGCAt                                                                                                                                                                                           >  2:140038/1‑251 (MQ=255)
                                                       cAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATAAAGAAGGCTGACCGTGACGATGATAAGCACCATA‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGACTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCAc                                                                                                                                                                                        >  1:33242/1‑251 (MQ=255)
                                                        aaGTGCATAGATATCCATGCCTCGCTTATCAGATAAAGGGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCCTC‑AGGGCGCAGGTTTGTTGGAACTGGACACCGCGAATCGCTTCCCACAGGGGGACGCCGATCCCGCCTACGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTAGCAGCGCTAGAGGGAGTAGGAGAACAGCAGAGGAATCAcc                                                                                                                                                                                       >  1:271171/1‑251 (MQ=255)
                                                        aaGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCAcc                                                                                                                                                                                       <  1:376240/251‑1 (MQ=255)
                                                          gTGCATAGATATCAATGCATCGCTTATCAGATCAAGTGCTTACGAAACCGTTGAGAGAGGAAATACAGCAGGACGACCGTGACGATGATAAGCAACATC‑AGGGCGCAGGATTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGACGACGACGGTTGCCGGGCGGACGGTGGATTCGAAGAGAGAGAGGGAGTAGGAGATCAGAAGAGGCATCACCtg                                                                                                                                                                                     >  1:93447/1‑251 (MQ=255)
                                                             cATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACGGGGACGAGGAAAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTgg                                                                                                                                                                                  <  2:362275/251‑1 (MQ=255)
                                                               tAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGgc                                                                                                                                                                                 >  1:169781/1‑250 (MQ=255)
                                                                                  taccAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAAACCAGCAGGCTGACCGTGACGATGATAAGCCCCATC‑AGGGCGCAGGTTTGGTGGAACTGGAACCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGAtc                                                                                                                                                               <  2:83104/246‑1 (MQ=255)
                                                                                     cAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGCTGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGCCGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTCGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTc                                                                                                                                                          >  2:130238/1‑251 (MQ=255)
                                                                                     cAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTc                                                                                                                                                          >  1:313708/1‑251 (MQ=255)
                                                                                      aGATAAAGTGCTTACGCAAACGTTGAGAGAGGAACTCCAGCAGGCTGACAGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGAACCCGCCTGCGCCGACCATGCCGACGACGGTGGCCGAGCGGACGTTGGATTCGAAGCGATGGCGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGGAGAGGATCGCTTCg                                                                                                                                                         >  2:266401/1‑251 (MQ=255)
                                                                                                cTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTgg                                                                                                                                               >  1:217650/1‑251 (MQ=255)
                                                                                                 ttACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATACGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTCGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGc                                                                                                                                              >  2:28257/1‑251 (MQ=255)
                                                                                                 ttACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGc                                                                                                                                              >  1:253498/1‑251 (MQ=255)
                                                                                                 caccgctaCCGTGGAGAGAGGAAATCCAGCTGGTTGAACGGGACGATGATAAGCACCAGC‑AGGGCGCAGTTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCGCAGTGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCTGTCGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGTAGCGCCTGTGGCAGCACGCCGGGGAGTATCGCTTAGAGTTTGTTGGc                                                                                                                                              <  2:233866/244‑1 (MQ=255)
                                                                                                    cGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGAGGATAAGCCCCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCAcc                                                                                                                                           <  1:128027/251‑1 (MQ=255)
                                                                                                      cAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCgg                                                                                                                                         >  2:252248/1‑251 (MQ=255)
                                                                                                          cGTTGAGCGAGGACATCCAGCAGGCTGACCCGGACGATGATATGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTACCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGTCGACGGTTGTCGAGCGGACGTTGGATTCGATGCGATAGAGGGAGTAGGAGTTCAGCATCGGCATCACCTGTGGCAGCACTCCGTCGAGGCTCTCTTCGTGCTTGTTGGCACCGGTgtc                                                                                                                                     >  2:249345/1‑249 (MQ=255)
                                                                                                          cGTTGAGCGAGGAAATACAGCAGGCTGACCGTGACGATGCTAAGCAC‑ATCCAGCGAGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTACCAACGGGTGACGCCGATCCCGCAAGCGCCGACCATGCCGACGACGGATGCCGCGCGGCAGTTGGATTTGAAGCGCAAGAGGGAGCAGGAGATCAGCCGTGGCATCACCTGTGGACGCACGCCGTAGAGGATCTCTTCGAGCTTGTGGGCaacgcaggc                                                                                                                                     >  2:402379/1‑243 (MQ=255)
                                                                                                                gagagTAAATCCAGCAGGCTGACCTTGACGATGAAAAGCCCCATC‑AGGGCGCAGGTTTGTTGGAACTGAAAACCGCCAATCGCTTCCCCCGGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGTTGGCGCGAAt                                                                                                                               <  2:64596/251‑1 (MQ=255)
                                                                                                                      taaTCCAGCAGGCTGACCGTCACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGATCTTGTTGGCACCGGTGGCGCGAATGCCTTc                                                                                                                         >  2:85025/2‑251 (MQ=255)
                                                                                                                               cggCTGACCGTGACGATGATAAGCCCCACC‑AGGGCGCAGGTTTGTGGGAACTGGCACCCGCGAATCGCTTCCCACAGGTTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTGCTAGCTTGTTGGCACCGGTGGCTCGAATGCCTTCCTccgggcc                                                                                                                <  2:24566/250‑1 (MQ=255)
                                                                                                                                   tGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACcgggccgggc                                                                                                            >  1:306680/1‑251 (MQ=255)
                                                                                                                                    taCCGTGACGATGATAAGCACCATC‑AGGGCTCAGGTTTGTTGGACCTGGAAACCGCGAATCGCTTCCCTCAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCTCCGGTGGCGCGAATGCCTTCCACcgggccgggc                                                                                                            <  1:237642/249‑1 (MQ=255)
                                                                                                                                        gTGACGCTGAGAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGCGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGGGGCGCGAATGCCTTCCACCGGGCCGGGCTCaac                                                                                                       >  2:41189/1‑250 (MQ=255)
                                                                                                                                        gTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACAGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGACGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGCGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAAt                                                                                                       >  1:378587/1‑251 (MQ=255)
                                                                                                                                         tGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCCCCGGGCCGGGCTCAATc                                                                                                      >  1:313296/1‑251 (MQ=255)
                                                                                                                                         tGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATc                                                                                                      >  1:220017/1‑251 (MQ=255)
                                                                                                                                                    aGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGc                                                                                           >  2:201173/1‑251 (MQ=255)
                                                                                                                                                       accttC‑AGGGCGCAGTTTTGTTGGAACTGGAAACCGCGAACCTCTTCCCACAGGGTGACGCCGAGCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTc                                                                                        <  1:68191/247‑1 (MQ=255)
                                                                                                                                                        cccatgggggCGCAGGTTGGTTGGACCGGGAACCCGCGAATCGCTTCCCAAAGGGGGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTc                                                                                        <  2:182809/244‑1 (MQ=255)
                                                                                                                                                         cATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCgg                                                                                      >  1:147042/1‑251 (MQ=255)
                                                                                                                                                                ggCGCATGTTTGTTGGAACTGGAACCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGTCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAagca                                                                                <  1:379755/251‑1 (MQ=255)
                                                                                                                                                                    cAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCtt                                                                            >  1:308509/1‑251 (MQ=255)
                                                                                                                                                                           gttgGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGc                                                                     >  1:7244/1‑251 (MQ=255)
                                                                                                                                                                                aaCTGGAACCCGCGACTCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGTCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGTAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTGCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGcc                                                                <  1:170657/251‑1 (MQ=255)
                                                                                                                                                                                 agtGGAAACCGCGAACCGCTTCCCACAGGGGGACGCCGATCCCGCCTGCGCCGACCAGGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGGAGGAGATCAGCAGTGGCATCACCGGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCg                                                               <  2:198569/249‑1 (MQ=255)
                                                                                                                                                                                     gAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGCTAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCCCCGCTTCGGAAAGCAGCTTGGAGAGCACGCCggggg                                                           >  2:29470/1‑248 (MQ=255)
                                                                                                                                                                                        aCCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCAACGGCGGCCCGAATGCCTACCACCGGGCCGGGCCCAATCGCTTCCACCGCTTCGAAAAGCAGCTTGGAGAGCACCCCGgtggtgg                                                        >  2:315661/1‑250 (MQ=255)
                                                                                                                                                                                                   cTTCCCACAGGGTAACGCCGATCCCGCCTGCGCCGACCCTGCCGACGACGGTTGGCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGAACCGGAGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTACACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATCAACAgg                                             >  2:258925/1‑250 (MQ=255)
                                                                                                                                                                                                     tCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATCGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGCGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAAAAgcgc                                           >  2:118419/1‑251 (MQ=255)
                                                                                                                                                                                                     tCCCACAGGGTGACGCCGATCCAGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAgcgc                                           >  2:364783/1‑251 (MQ=255)
                                                                                                                                                                                                           aGGGTGACACCGATCCCGCCTGCGCCGACCAAGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCAACTGTGGCAGCACGCCGTCGAGGATCTCTTCGAGATTGTTGGCACCGGTGGCGCGAGTGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGCAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGcgccagcac                                     >  2:203971/1‑251 (MQ=255)
                                                                                                                                                                                                                 aCGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGCACGATGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGAGGCATCACCTGTGGCAGCACGCAGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGCATGACTTCCACCGGGCCGGGCTCAATCGATTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGGGTGGATAAACAGcgccagcacgccagc                               >  1:34511/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccGACCCCGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGccagcacgccagcacg                            <  1:130238/251‑1 (MQ=255)
                                                                                                                                                                                                                        tCCCGCCTGCGCCGACCATGCCGCCGACGGTTGCCGAGCGGACGTTGAATTCGAAGCGGTAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGGGGTGTGGATAAACAGCGCCagcacgccagcacgcc                          >  1:227733/1‑249 (MQ=255)
                                                                                                                                                                                                                           cGCCTGCGCCGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCCTCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGCTAAACAGCGCCAGCACGCCAGCCCGCCAGCg                      >  2:343562/1‑250 (MQ=255)
                                                                                                                                                                                                                           cGCCTGCGCCGACCACGCCGACGACTGTTGCCGAGGGGAGGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGa                     <  1:140038/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tccGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGGGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCCGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCCGCACGCCAGCCCGCCAGCGCCCGggcc              >  2:180479/2‑251 (MQ=255)
                                                                                                                                                                                                                                    cGACCATGCCGACGACGGTTGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGGGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGggccg             >  2:307924/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gcggacgGTTGCCGAGCGATCGTTGATTTCGAGGAGATAGAGGGAGTAGGGTATGAGCAGAGGCATCCCCTGTGGCAGCTCGCCGTAGAGGATCGCTGCGAGCTTTTTGGCACCGGTGGCGCGAATGCCTTCCACCGTGACGGGCTCAATCGCTTCCGCCGCTCCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGc  <  2:1911/248‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CAAGTTTTGCGGCTATCTCTTGATAGCGTGTTGGGTAGCTGGTCGGATGTGTAGACAAGTGCATAGATATCAATGCCTCGCTTATCAGATAAAGTGCTTACGCAACCGTTGAGAGAGGAAATCCAGCAGGCTGACCGTGACGATGATAAGCACCATC‑AGGGCGCAGGTTTGTTGGAACTGGAAACCGCGAATCGCTTCCCACAGGGTGACGCCGATCCCGCCTGCGCCGACCATGCCGACGACGGTCGCCGAGCGGACGTTGGATTCGAAGCGATAGAGGGAGTAGGAGATCAGCAGTGGCATCACCTGTGGCAGCACGCCGTAGAGGATCTCTTCGAGCTTGTTGGCACCGGTGGCGCGAATGCCTTCCACCGGGCCGGGCTCAATCGCTTCCACCGCTTCGGAAAGCAGCTTGGAGAGCACGCCGGTGGTGTGGATAAACAGCGCCAGCACGCCAGCACGCCAGCGAACGGGCCGAGGCCGACGGC  >  NC_000913/4322359‑4322846

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