Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,314,321 (G)5→6 intergenic (+124/+23) yjdP → / ← phnP protein YjdP/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,314,3211.G98.1% 172.5 / NA 53intergenic (+124/+23)yjdP/phnPprotein YjdP/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (28/24);  minor base A (0/1);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCGGTGTCAGACAG  >  NC_000913/4314029‑4314619
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                           
cgGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCNCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCACTAACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGNACAATCGGTCCACTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGGGCCAgcg                                                                                                                                                                                                                                                                                                       >  2:700875/1‑300 (MQ=255)
                   aGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAAGGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCtcatca                                                                                                                                                                                                                                                                                    <  2:498087/300‑1 (MQ=255)
                       acgaGCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGcc                                                                                                                                                                                                                                                                                <  2:482638/300‑1 (MQ=255)
                                     cAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccat                                                                                                                                                                                                                                                                  >  2:734226/1‑300 (MQ=255)
                                       atatGACGCCCGCCTCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCTGAGAATGGAAGATGAGTATGGTCTATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAAAGGGTTAGGGTGAGGGGGGCCAGCACACACCTCCCCTCTCATCACGCCACCCCAATCTccaccc                                                                                                                                                                                                                                                                <  2:490205/300‑4 (MQ=255)
                                                                  aGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATNGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATNGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCCCAAACCCTGACg                                                                                                                                                                                                                                     <  2:616624/300‑1 (MQ=255)
                                                                               cagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGGGACAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCCTCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTc                                                                                                                                                                                                                        >  2:12585/1‑300 (MQ=255)
                                                                               cagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCACTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCAACCCAATCTCCATACCATCAAACCCCAACTCAAACCCTGACGGTAGTGCGTTTTc                                                                                                                                                                                                                        >  2:580987/1‑300 (MQ=255)
                                                                                 gCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                                                                                                                                                                                                      >  1:85033/1‑300 (MQ=255)
                                                                                 gCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCa                                                                                                                                                                                                                      >  2:724213/1‑300 (MQ=255)
                                                                                  cTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGATTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACANTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCACTATAATCAAGCCACCCAAATCTCAATCCCATCAAACCCCAACTCAAACCCTGACAGTAGTGCGTTCTCCAt                                                                                                                                                                                                                     >  2:180909/1‑300 (MQ=255)
                                                                                              cagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCNATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATc                                                                                                                                                                                                         <  1:519348/300‑1 (MQ=255)
                                                                                                gcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCaa                                                                                                                                                                                                       >  1:213742/1‑300 (MQ=255)
                                                                                                  gccagcNTGAGCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGNGGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGTAGAAAACGAGCGCAATAGTCAAGATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCTATCAAACCCCACCTCAAACCCTGACGGTAGTGTGTTTTCCATCAGCCATGCATCagac                                                                                                                                                                                                     <  2:633774/300‑4 (MQ=255)
                                                                                                     agcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGTAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACtgg                                                                                                                                                                                                  <  2:660921/300‑1 (MQ=255)
                                                                                                                         cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAAt                                                                                                                                                                             <  1:102918/300‑1 (MQ=255)
                                                                                                                         cAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATNCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCNGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAAt                                                                                                                                                                             <  2:608328/300‑1 (MQ=255)
                                                                                                                         cAACGGGAGGAGGAAGAGCGGAGACTGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATGCAATATCGAGCACAATCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTTGCTGATATGGTTCAGAAt                                                                                                                                                                             <  2:396964/300‑1 (MQ=255)
                                                                                                                         aaaCTTGAGGATGAGGGGCGTGGAATGGAAGGTGAGTATGTTTGCTGGGGTTGGGTATGGAGCAGGCATCTCCCCTCGCCCCACTGTGTAGAGGTATAAAACGCGCGCAATATTCAATATCAAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGGGCCAGCGTGCACATTCCCTCTCATCACGCCACCCTAATCTCCATCCCATCAAACTCCACCTCAAACCCCGACTTTAGTGCCTTTTCCATCAGCCACGCATCAAACTGGTGGCCGatacgattcacacc                                                                                                                                                                             <  2:47226/299‑12 (MQ=255)
                                                                                                                                    tGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCg                                                                                                                                                                  <  1:399681/300‑1 (MQ=255)
                                                                                                                                               gaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCAAATGGGGAGGGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAAGGTAGTGTGTTTTCCATCCGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGCGGCGAGCAGATAAcc                                                                                                                                                       >  2:32970/1‑300 (MQ=255)
                                                                                                                                                     aaGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTc                                                                                                                                                 <  1:15292/300‑1 (MQ=255)
                                                                                                                                                          gAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAgcgc                                                                                                                                            >  1:211738/1‑300 (MQ=255)
                                                                                                                                                          gAGTACGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAGTGGGGAGAGGGAGAAAACGAGCGCAGTATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCACCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAgcgc                                                                                                                                            <  2:46866/300‑1 (MQ=255)
                                                                                                                                                                  gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACgg                                                                                                                                    >  1:340023/1‑300 (MQ=255)
                                                                                                                                                                    gATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCAACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTa                                                                                                                                  >  1:494454/1‑300 (MQ=255)
                                                                                                                                                                        aGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATTAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTa                                                                                                                              <  1:276733/300‑1 (MQ=255)
                                                                                                                                                                        aGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTa                                                                                                                              >  1:738899/1‑300 (MQ=255)
                                                                                                                                                                          gTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGTAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGNGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTaaa                                                                                                                            <  2:662790/300‑1 (MQ=255)
                                                                                                                                                                          gTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGTGTAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTaaa                                                                                                                            <  1:202121/300‑1 (MQ=255)
                                                                                                                                                                          gTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTaaa                                                                                                                            <  1:263188/300‑1 (MQ=255)
                                                                                                                                                                            tGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATc                                                                                                                          <  1:49493/300‑1 (MQ=255)
                                                                                                                                                                            tGGGTANGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGNGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGTGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATc                                                                                                                          <  2:629134/300‑1 (MQ=255)
                                                                                                                                                                              ggTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATcac                                                                                                                        >  1:476537/1‑300 (MQ=255)
                                                                                                                                                                                         ggCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTAcgc                                                                                                             >  1:482337/1‑300 (MQ=255)
                                                                                                                                                                                           cATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCgg                                                                                                           <  1:383384/300‑1 (MQ=255)
                                                                                                                                                                                                   cTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCcg                                                                                                   <  1:20978/300‑1 (MQ=255)
                                                                                                                                                                                                       ccccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAGATCACAGTGATTACGAGGTGCAGCcgcgcg                                                                                               >  2:371698/1‑300 (MQ=255)
                                                                                                                                                                                                       ccccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACACGTGGAGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCcgcgcg                                                                                               >  2:474726/1‑300 (MQ=255)
                                                                                                                                                                                                        cccAATGGGTAGAGGGAGAAAACGAGCGCGATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTACATCcgcgcgc                                                                                              <  2:105952/300‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaCGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATcc                                                                            <  1:735764/300‑1 (MQ=255)
                                                                                                                                                                                                                                       tatTCAATATCGAGCACAATCGGTCCCCTCGCCCCCCTGGGGAGAGGGTTAGGGTGAGGGGGACCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGTTAGTGCGTTTTCCATCAGCCACGCATCAAACTGTGGTCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCt                                                               <  2:245274/300‑1 (MQ=255)
                                                                                                                                                                                                                                        gttCAATATCGAGCACAATCGGTCCCCTCGCNCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCANACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGTCGGGTGACTGCCATCCATCACCATTACCTg                                                              >  2:548495/2‑300 (MQ=255)
                                                                                                                                                                                                                                        atTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTg                                                              >  1:677848/1‑300 (MQ=255)
                                                                                                                                                                                                                                        atTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTg                                                              >  2:478535/1‑300 (MQ=255)
                                                                                                                                                                                                                                                               ccccTCGCCCCTCTGGGGAAAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGCTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGTGCGGCGGGTGACTGCAATCCATCAACATTACCTGCGGCTGATTATTGCGTAAAAAtt                                       >  2:145637/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ccTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGTGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTACATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCNCGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGttt                              >  2:671613/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ccTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGttt                              >  1:117170/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                        ccTCTGGGGAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCAACATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGttt                              >  2:219804/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             ggggAGAGGGTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTcc                         >  1:627808/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                    gggTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGc                  >  1:452157/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                    gggTTAGGGTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGAGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGAAAGc                  >  2:484434/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                          gggTGAGGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCgg            >  1:21198/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                              aggggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCGGTGTCa       >  1:108992/3‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                 ggggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCGGTGTCAGAc    <  1:54291/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   ggCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGTTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCGGTGTCAGACAg  <  2:145978/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                           
CGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATCACAGTGATTACGCGGTGCATCCGCGCGCGGCGGGTGACTGCAATCCATCACCATTACCTGCGGCTGATTATTGCGTAAAAATTTCAGCGTTTTTTCCGGCAAGCCTGCGGTGTCAGACAG  >  NC_000913/4314029‑4314619

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