Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,546,601 C→T I502I (ATC→ATT fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,546,6010CT98.3% 201.8 / NA 58I502I (ATC→ATTfimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/24);  minor base G (0/1);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546359‑4546847
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGCCAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAacacc                                                                                                                                                                                                                                                >  2:312093/1‑250 (MQ=255)
           cGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGtt                                                                                                                                                                                                                                      <  2:142958/250‑1 (MQ=255)
            ggTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTAt                                                                                                                                                                                                                                    <  1:501662/251‑1 (MQ=255)
                          cATGGGGGCCCTGGGCGCTCTGCCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGCGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGATACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                                                                                                                                                                                       >  2:353694/1‑250 (MQ=255)
                             gcgggCACTGGGCGCTCTATCTGTGGATCTGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGCCGGACAATCGGTGCGTTTTCTCTATAACAACTCGCTCCATGAATCAGGCACGAGTATTCAGTTATTGGGTTACCGTTATTCGTCCAGCGGATCTTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAACCACAGGACGGAGTTATTCAGGTTAAGCCGAAAt                                                                                                                                                                                                                   >  2:86158/3‑251 (MQ=255)
                              ggggCCCTGGGCGCTCTGTCTGTGGATATGACGCCGGCTAATTCCACCCTTCCCGATGACAGTCAGCAGGACGGATAATCGGTGCGTTGTCTCTATAACAAATCTCTCAATGAATCAGGTACGAATAGTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAGCATACAGTCGAATGAATGGCTACAACATTGAAACACAGGTCGTAGGTATTCAGGTTAAGCCTAAAtt                                                                                                                                                                                                                  <  2:464326/251‑1 (MQ=255)
                                  cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                                                                                                                                                                              <  1:66957/251‑1 (MQ=255)
                                  cACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTAGTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                                                                                                                                                                              <  1:514450/251‑1 (MQ=255)
                                   aCTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTGCGCTGAGACAACATACAGTCGAATGAATGGCGACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                                                                                                                                                                             <  1:24776/251‑1 (MQ=255)
                                            tctgtctgTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                                                                                                                                                                                     >  1:74940/1‑250 (MQ=255)
                                             ctgtctgtGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                                                                                                                                                                                     <  1:432115/249‑1 (MQ=255)
                                                   gtgGATATGAAGCAGGCTAATTCCACACGTCCCGATGACAGTCAGCATGACGGACAATCGTTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAGTATTATGTTAGAGAGTTCCCGTTATTCGACCAGCGACTATTTTAATTTCGCTGATACAACATACAGTCGACAGAATAGCTACAACATTGAAACACATGACGGACATATCCAGATTAAGCCGAAATTCACCGCCTATTAACACCTCGc                                                                                                                                                                                             >  2:420140/1‑251 (MQ=255)
                                                         atGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCTTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                                                                                                                                                                                       <  1:202384/251‑1 (MQ=255)
                                                           gACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAAGCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACa                                                                                                                                                                                     <  1:501868/251‑1 (MQ=255)
                                                              gCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                                                                                                                                                                                  <  1:455067/251‑1 (MQ=255)
                                                               cAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                                                                                                                                                                                 <  2:210016/251‑1 (MQ=255)
                                                                 ggCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                                                                                                                                                                               >  1:56379/1‑251 (MQ=255)
                                                                  gCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                                                                                                                                                                              <  1:11034/251‑1 (MQ=255)
                                                                        tctacacTTCCCGATTACAGTCAGCAAGACGGACACTCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATGCAGTTAGTGTGTTACCGTTATTCGACCAGCGGATATGTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATGGAAACACAGGACGGATTTATTCATGTTAAGTCGAAATTCACCGTCTAGTACAACCTCGCTTATAACAATCGCGGGAAAt                                                                                                                                                                         <  2:65791/247‑1 (MQ=255)
                                                                         ccACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt                                                                                                                                                                        >  2:109117/1‑250 (MQ=255)
                                                                          cacaCTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc                                                                                                                                                                      <  1:321514/251‑1 (MQ=255)
                                                                          cacaCTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc                                                                                                                                                                      <  2:56379/251‑1 (MQ=255)
                                                                                tCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc                                                                                                                                                                 >  1:489362/1‑250 (MQ=255)
                                                                                  ccggtGACATTCAGAATAACGGATACTCAGTGCGTTTTCTCTATAACACATCGCCCAATGCATCATGCACGAATATTTCGTTAGTGGGATACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACCACATACAGTCGCATGAATGCCTACAACCTTGAAACACAGGACTGGGTTATTCAGATTAACCCGAACTTCATCGACTATTACAACCTCGCTTATAACAATCGCGGGCAATTACAACTCAcc                                                                                                                                                              >  2:195006/5‑251 (MQ=255)
                                                                                          aGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTAAATGAATCAGGCACGAATATTCAGTTAGCGGGTTACCGTTATTCGACCCGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTc                                                                                                                                                       >  1:292096/1‑250 (MQ=255)
                                                                                             cAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCa                                                                                                                                                   <  1:470645/251‑1 (MQ=255)
                                                                                                      ggACAATCGGTGCGTTTTCTCTATAACAAATTGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAAGGGCTACAACATTGAAACACGGGAAGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGAGCTGATAACAAACGCGGGAAATTACAAATCACCGTTAATCAGCAACTCGGgcg                                                                                                                                          >  2:16646/1‑251 (MQ=255)
                                                                                                         cAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCCGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGcac                                                                                                                                       >  2:336442/1‑251 (MQ=255)
                                                                                                             cGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGTCACGACTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTGTAATTTAGCTGATACAACATACAGTCGAATGAATGGCTCCAACATTGAAACACAGGACGGAGTTATTCAGGTTACGCCGTAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGCAAATTACACATCACCGTTACTCAGCAACTCGGGCGCACATc                                                                                                                                    >  2:274471/1‑250 (MQ=255)
                                                                                                               gTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCaa                                                                                                                                  >  2:167233/1‑250 (MQ=255)
                                                                                                                  cGTTTTCGCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACt                                                                                                                              <  1:287669/251‑1 (MQ=255)
                                                                                                                          ctATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCATCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGATTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACACGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTAtt                                                                                                                         <  1:312832/248‑1 (MQ=255)
                                                                                                                                  aaTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                                                                                                                                                                                  >  2:262408/1‑151 (MQ=255)
                                                                                                                                  aaTCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                                                                                                                                                                                  <  1:262408/151‑1 (MQ=255)
                                                                                                                                      gCTCAATGAATCAGGCACGCATATCCAGCTAGAAGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                                                                          >  2:263948/1‑251 (MQ=255)
                                                                                                                                            tGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGAGGTAGCCATCAAACtt                                                                                                    >  2:310910/1‑251 (MQ=255)
                                                                                                                                                 cAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACAATGTATTTGAGTGGTAGCCATCAAACTTATTgg                                                                                               >  2:299564/1‑251 (MQ=255)
                                                                                                                                                              tCAGTTAGTGGGTTACCGTTAGTCGACCCGCGGCTATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATg                                                                                  >  2:103652/1‑251 (MQ=255)
                                                                                                                                                                tgTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAACCCCAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGTAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGAATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAATAATga                                                                                 >  2:100321/2‑249 (MQ=255)
                                                                                                                                                                      tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGCGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                                                                            >  1:286265/1‑249 (MQ=255)
                                                                                                                                                                       gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACCACATAAAGTCGAATGAATGGCTACCACCTTGAAACACAGGACGGCGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACCCTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                                                                         >  1:219185/1‑251 (MQ=255)
                                                                                                                                                                         gTTACCGTTATTCGAACAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCTCCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGTCATCAAACTTATTTGGGAACGAGTAATTTCGATGAGCaa                                                                       <  2:26859/251‑1 (MQ=255)
                                                                                                                                                                                      gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAAGGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt                                                          <  2:286265/251‑1 (MQ=255)
                                                                                                                                                                                         cAGCGGATCTTTTAATTTAGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCAAGGCTGGATTa                                                        >  1:394727/1‑250 (MQ=255)
                                                                                                                                                                                          aGCGGATATTTTAATTTCGCTGATACAACATACTGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAAAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa                                                       >  1:399417/1‑250 (MQ=255)
                                                                                                                                                                                             ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATAc                                                   >  1:365152/1‑251 (MQ=255)
                                                                                                                                                                                               atatTTTAATTTAGCTAATACAACATCCCGTCGAATGAATGCCTACAACATTGCAACACAGGACGGAGTTCTTCAGGTTAAGCCGACATTCACCGAATCTTACAACCTCGCTTATAAAAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGAGCGCACATCAACACTGTATTTGGGTGGAAGTCAACAAACTTATTGGGGAACGAGTAATGTCGAAGAGCAATTCCAGGCTGGATTAAATaatg                                                 >  2:490992/1‑248 (MQ=255)
                                                                                                                                                                                                  ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGt                                              >  1:147314/1‑251 (MQ=255)
                                                                                                                                                                                                   tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGtt                                             >  1:228399/1‑251 (MQ=255)
                                                                                                                                                                                                         ttCGCTGATACAACATACAGTCGGATGAATGGCTACAACCTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGGAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGCAATGGCGATGCGAAGTTCCAGGCTGGGTCAAATAATGCGTTCGa                                          >  2:425816/1‑248 (MQ=255)
                                                                                                                                                                                                          tCGCTGATACAACATACCGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGa                                          >  2:244306/1‑247 (MQ=255)
                                                                                                                                                                                                               gATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCa                                  <  1:310910/250‑1 (MQ=255)
                                                                                                                                                                                                                               gaatgaatGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATCCTGCGTTCGAAGATATCAACTGGACGCTCAGCtat                 <  1:468687/251‑1 (MQ=255)
                                                                                                                                                                                                                               gaatgaatGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCtat                 >  2:416836/1‑251 (MQ=255)
                                                                                                                                                                                                                                aatgaatgGCTTCAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCTAAATTCACCCACTATGTCAACCTCGCTTATAACAAACGCGTTAAAGTACAACTCACCGTTACTCAGCATCTCGGGCTCATATCAACACTTGATTTGAGTGTTAGCCATCAAACTTATTGGGGAACGATTAATTTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCATCTGGACGCTCAGCtata                <  1:205416/251‑1 (MQ=255)
                                                                                                                                                                                                                                           acaacaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCCTTACTCAGCAATTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACATATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACg                          <  1:514382/230‑1 (MQ=255)
                                                                                                                                                                                                                                              acaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAc  >  1:202582/1‑251 (MQ=255)
                                                                                                                                                                                                                                              acaTTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAc  >  1:453738/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GCTTTTAATTTCGGTATCGGGAAAAACATGGGGGCACTGGGCGCTCTGTCTGTGGATATGACGCAGGCTAATTCCACACTTCCCGATGACAGTCAGCATGACGGACAATCGGTGCGTTTTCTCTATAACAAATCGCTCAATGAATCAGGCACGAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAAATACTGCGTTCGAAGATATCAACTGGACGCTCAGCTATAGCCTGACGAAAAAC  >  NC_000913/4546359‑4546847

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