Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,292,312 C→T P324P (CCC→CCT yraK → putative fimbrial adhesin YraK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,292,3120CT98.3% 194.8 / NA 59P324P (CCC→CCTyraKputative fimbrial adhesin YraK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/28);  minor base G (0/1);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTT  >  NC_000913/3292075‑3292557
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
cAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCg                                                                                                                                                                                                                                          >  1:200056/1‑251 (MQ=255)
 aaTGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCgg                                                                                                                                                                                                                                         >  1:271130/1‑251 (MQ=255)
    gccatgccaGACGGTTTTACTCCAGTAGACTTTGATATCGCTTATGACGGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCAGCGACGGGGCAACTGGGGTAGTAGACCAGTACAACCGGGTCGCCAGGCTAAGAATTTCATACAATGTGCCCGATGTTGGGATTCGTGTTGAAAATCTCGGCGGCGGAGTTGTAATTATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCa                                                                                                                                                                                                                                      <  1:343087/251‑1 (MQ=255)
         gccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGGGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTAGTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCGTCCTGTTGATCCTTCCGGGCATGGc                                                                                                                                                                                                                                  <  2:286768/250‑1 (MQ=255)
            aGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGAGCTACTTCAAAGATTAAAAACTCGTTGCAACTGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGACGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCg                                                                                                                                                                                                                              >  1:34297/1‑251 (MQ=255)
               cGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTc                                                                                                                                                                                                                            >  1:495155/1‑250 (MQ=255)
                ggTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCaa                                                                                                                                                                                                                          <  2:495155/251‑1 (MQ=255)
                ggTTTTACTCCAGTAGACTTTGAGATCACTTATGACTGTGGGGATACTTCAAAGATTAAAAACTCGTTGCAAAGGCGCATCGACGGTTCAACTGGGGGTGTAGACCAGTTCAACCTGGTCGCCAGGCGAGGAAGTTCAGACAGTGTGCCCGATGGCGTTTTTCGTAGTGAAAATCGCTGCGGCTGAGTTGCAAATATTCCTTGTCATAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCaa                                                                                                                                                                                                                          <  1:294236/251‑1 (MQ=255)
                                    tGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCGTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGtt                                                                                                                                                                                                       <  2:277349/250‑1 (MQ=255)
                                       tatCACTTATGACTGTGGTGATACTTCAAAGATTAAACACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATc                                                                                                                                                                                                   >  2:276589/1‑251 (MQ=255)
                                          cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTg                                                                                                                                                                                                 <  1:144712/250‑1 (MQ=255)
                                           actgatgacgttggggATACTTCAAAGATTAAAAACCCGTTGAAAATGCGCACCGACGTTACAACTGGGGTGGTAGACCAGTACAACCTGGTCGCAGGGGGAAGAATTTCAGACAATGGGCCCGAGGTCGGTATTGGTATTGAAAATCTCGGCGGCGGAGTGGAAAAGATTCCTTGTCAGAACGGTATCCTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTc                                                                                                                                                                                              <  2:343832/236‑1 (MQ=255)
                                           aCTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCCGTTACTCTGGt                                                                                                                                                                                               >  2:134348/1‑251 (MQ=255)
                                              tATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTAGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCg                                                                                                                                                                                             <  1:193922/250‑1 (MQ=255)
                                               aTGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCggt                                                                                                                                                                                           >  1:574550/1‑251 (MQ=255)
                                               aTGACAGTGGTGATACTTCAAAGCTTAAAAACTCGTTGCAAATGCGCATCGACGGTCCAACTGGGGTAGTAGACCAGAACAACCTGGTCGCAAGGCGCAGAAGTTCAGACAATGTGCACGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGAAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTACTCTGGTCggt                                                                                                                                                                                           >  2:350347/1‑251 (MQ=255)
                                                  aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCggtggg                                                                                                                                                                                        >  1:333981/1‑250 (MQ=255)
                                                    tgtgGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGAAGGTACCACTGGGGTAGTAGACCAGTACAACCTGGTCGCACGGCGAAGCAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGgtggtg                                                                                                                                                                                       >  1:83045/1‑250 (MQ=255)
                                                       ggtgATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCAGCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGGATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCt                                                                                                                                                                                   <  2:393008/251‑1 (MQ=255)
                                                              cTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAAc                                                                                                                                                                             >  2:225604/1‑250 (MQ=255)
                                                                           aaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTg                                                                                                                                                                                                 <  1:143970/217‑1 (MQ=255)
                                                                           aaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTg                                                                                                                                                                                                 >  2:143970/1‑217 (MQ=255)
                                                                             aactcgggCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGc                                                                                                                                                             <  1:512036/244‑1 (MQ=255)
                                                                                cgtagCAAATGCGAATAGACGGCAGACCTGGGGTCGATGATCAGTACAACCTAGTCGCAAGGCGAAGAAGATCAGACAATGTGCCCGATGTCGGTATTCGTATTGACAATATCGGCGGCGGAGGTGAAACTATTCCTTTTCAGAAAGGTCTCCTGCATGTTGATTCTTCCGGGCATGGCACCGTCGATATGCAAGCCTGGCCCGTTCAGCTGGTCGGTGGTGCGCTGGAACCAGGAAAATTTCTGGGcacc                                                                                                                                                          >  2:283159/5‑250 (MQ=255)
                                                                                cGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGcac                                                                                                                                                           <  2:572517/250‑1 (MQ=255)
                                                                                      aaaTGCGCATCGACGGTACAACTGGGGTAGTCGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGAATTCGTATTGAAAATCTCAGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGGATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGCACCAGGAAAATTTCAGGGCACGGccacc                                                                                                                                                    >  1:297794/1‑251 (MQ=255)
                                                                                          gcgcATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGGCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCTGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCCCCATTa                                                                                                                                                <  2:164107/251‑1 (MQ=255)
                                                                                          gcgcATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGCAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCAtt                                                                                                                                                 >  2:144070/1‑250 (MQ=255)
                                                                                             cATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCg                                                                                                                                             <  1:421385/251‑1 (MQ=255)
                                                                                                cGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTc                                                                                                                                           >  1:430371/1‑250 (MQ=255)
                                                                                                    ggTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTGCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGGTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATGTCAGGGCACAGCCACCATTACCGTCATCGt                                                                                                                                      <  2:83045/251‑1 (MQ=255)
                                                                                                     gTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTg                                                                                                                                     >  1:341743/1‑251 (MQ=255)
                                                                                                               gggagaagaacAGTACAACCGGCTCGCCAGGCGAAGAAGTTCAGACAATGAGACCGATGCCAGAGTTCGTATTGAAAATCTCGGCGGCGGAGTTGAAAATAGTCATTTTCAGAACGGTATCATTCCTGATGATACTTCCGGGCATGGCACCGTCAATAAGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTAGAAACAGGTAAATTTCAGGACACAGCCACCATTCCCGTCATCGAGCGGGAAaagc                                                                                                                           >  2:186529/11‑250 (MQ=255)
                                                                                                                   gtagACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGAc                                                                                                                       <  2:532524/251‑1 (MQ=255)
                                                                                                                             aCAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGATGTGCAAATATTCCTTGTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGttt                                                                                                              <  1:571258/250‑1 (MQ=255)
                                                                                                                                        gcaaGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTGCCGGGCAGGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCGGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTaacgcta                                                                                                  <  2:324326/248‑3 (MQ=255)
                                                                                                                                        gCCAGGCGAAGAAGTTCAGACAAGGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTACCTGTTGATCCTTCCGGGCATGGCACCGTCAAGATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTGCTGGGCACAGCCGTCACTACCGTCATCGTGCGGTAAAAGAAGACGTTGTGTGTTGATGaacgcaa                                                                                                  <  2:34297/251‑1 (MQ=255)
                                                                                                                                           aGGCGAAAAAGTTCAGACAATGTGCCCGAAGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTACTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGACATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTaacgccacgc                                                                                               >  2:377161/1‑251 (MQ=255)
                                                                                                                                             gCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCtt                                                                                             >  1:271500/1‑251 (MQ=255)
                                                                                                                                                aagaagTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGATCGGTTTCCTTCCTGTTGATCCTTCCGGGCATGGCCCCTTCAATATGCGCGCCTGGCCTGTTAATCTGTTCGGTGGTGAGCTGGTAACAGGAAAATTTCAGGTCTCAGCCCCCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCATCGCTTCgg                                                                                          <  1:222181/251‑1 (MQ=255)
                                                                                                                                                   aagTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCAc                                                                                       >  1:369816/1‑251 (MQ=255)
                                                                                                                                                        cAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATAGTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGtt                                                                                   <  2:182707/250‑1 (MQ=255)
                                                                                                                                                         aGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCCATATGCGCGCCTGGCCCGTTAATCTGGTCGGTGGTGAGCTGGCAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGTTTAACGCCACGCTTCGGCACTGTTAc                                                                                 >  1:339180/1‑251 (MQ=255)
                                                                                                                                                            cAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTAcccc                                                                              >  1:399219/1‑251 (MQ=255)
                                                                                                                                                                 tgCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCTAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCctgctg                                                                         >  2:275722/1‑251 (MQ=255)
                                                                                                                                                                         gTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACt                                                                                      >  2:188135/1‑230 (MQ=255)
                                                                                                                                                                         gTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACt                                                                                      <  1:188135/230‑1 (MQ=255)
                                                                                                                                                                          tCGGTATTTGTAATGAAAATCTAGGCGGCGGAGATGAAAATATTGATTTTCAGAACGGAATACTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGACTGGCAAGATAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGACACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGATGCTCCAgcgc                                                                >  2:398533/1‑251 (MQ=255)
                                                                                                                                                                                        gAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAgcg                                                  >  2:471039/1‑251 (MQ=255)
                                                                                                                                                                                                cggcggcggAGTTGCATATATTCCTTTTCAGAACGGTATCCTGCCTGTTGATCCTTGCGTGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAAGCTGGTCGGTGGTGAGCTGGAAACAGTAAAAGTTCAGTGCACAGCCTCCATTACCGTCAGCGTTCGGTAATATAAGACGTTGCGTTTAGATTAACGTAACTCTACGGCACTGGTACCCCGGCTGCTCCAGCGCATACTTAGACAGCGCAttt                                             <  2:32531/248‑1 (MQ=255)
                                                                                                                                                                                                  gcggcggAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTttctt                                        <  1:461096/251‑1 (MQ=255)
                                                                                                                                                                                                  gcggcggAGTTGCAAATATTCCGTTTCAGAACGGTTTCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCTGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTttctt                                        <  2:333981/251‑1 (MQ=255)
                                                                                                                                                                                                   gggcggAGTTGCAAATATTCCTGTTCAGAACGGTATCCTTCCTGTTTATTCGTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCGCAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGAGTAACGCAACGCTTGGGCTCTGTTACCCCTGCTGCTCCATCGCATACGTATACAGCCCAGTTttcttc                                       <  2:339180/250‑1 (MQ=255)
                                                                                                                                                                                                       ggAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGc                                   >  2:138834/1‑251 (MQ=255)
                                                                                                                                                                                                           ttGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATCTGCGAGCCTGGCCAGTTAATCTGGCCGGTGGTGAGCTGGAAACAGGCAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACATATACAGCGCATTTTTCTCCCCGCCGTg                               >  1:261704/1‑251 (MQ=255)
                                                                                                                                                                                                            tGCAAATATTCCTTTTCAGAACGGTATCCTTCCGGTTGATCCTTCCGGGCATGGCACCGGCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTg                               <  2:271130/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCgg                    <  1:144070/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGcgccgcc         <  1:276589/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGACACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTCTACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCttttt  >  2:529623/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
CAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTT  >  NC_000913/3292075‑3292557

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: