Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,003,462 C→A A49S (GCG→TCG)  fliC ← flagellar filament structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,003,4620CA96.5% 235.9 / NA 86A49S (GCG→TCG) fliCflagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (35/48);  total (35/51)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACC  >  NC_000913/2003220‑2003705
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
cTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCg                                                                                                                                                                                                                                             >  1:198879/1‑251 (MQ=255)
  gACCAGATAGGCGGTCAATTTCATCCAGACGGGAGTGAATTTCGGCCTGGATAGCAGACAGATCAGACTCAGAGTGAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTGCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCt                                                                                                                                                                                                                                           <  2:653922/251‑1 (MQ=255)
  gACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTGAGAGGTGAAACGGTTAGCAATCGACTGACCCGCt                                                                                                                                                                                                                                           <  1:654170/251‑1 (MQ=255)
       gATACGCGGTCAATTTCATCCAGACGGGATTTACTTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTATGCGCCACGGAGATCCCGTCGTTGGCTTTACGGGACGCCTGAGTCAGGCCTTTAATGTTAACGGTGAAACGGTTAGCCATCGACTGACCCGCTTCGTc                                                                                                                                                                                                                                      >  2:348408/1‑251 (MQ=255)
             cgGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCtt                                                                                                                                                                                                                                <  2:22058/251‑1 (MQ=255)
              ggTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTc                                                                                                                                                                                                                               >  1:301526/1‑251 (MQ=255)
              ggTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCAAACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGGGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAAAGTTAGAGGTGAAGCGGTTAGCAATCGAATGACCCGCTGCGGAATCCTTc                                                                                                                                                                                                                               >  2:281138/1‑251 (MQ=255)
                     ttCATCCAGACGGGCTTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGGGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCgcgcgtt                                                                                                                                                                                                                        >  2:345601/1‑248 (MQ=255)
                                  gATTTAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGAAAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAAGGTTAGAGGTGGAACGGTTATCAATCGACTGACCCGCTGAGTCATCCTTCGCGCTGTTAATACGCAAGcc                                                                                                                                                                                                           >  2:227167/1‑251 (MQ=255)
                                   aTTTAATTTCGTCCTGGATAGCAGACAGATCAGTCTCAGTGTTAGTACCGGTAGTGGCCTGTACCGTTAGTTCACGCACCCGCTGTAAGTTGTTGTTGATTTCGGATAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa                                                                                                                                                                                                          <  2:665852/251‑1 (MQ=255)
                                      tAATTTCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGATGTTGTTGATTTCGGACAGCACGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGGGTTACGGGCCGCCTGAGTCAGGCCTTTTATGTTAGAGGAGATACGGTTAGCAATCGACTGACCCGCTGCGGCATCCTTAGCGATGTTAATACGCAACCCaga                                                                                                                                                                                                        >  1:105918/1‑250 (MQ=255)
                                         tttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCCACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCagaaaac                                                                                                                                                                                                    >  1:659104/1‑248 (MQ=255)
                                         tttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa                                                                                                                                                                                                          >  1:443997/1‑245 (MQ=255)
                                         tgtcggccTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCa                                                                                                                                                                                                          <  2:443997/239‑1 (MQ=255)
                                          ttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAagac                                                                                                                                                                                                    >  1:569991/1‑250 (MQ=255)
                                          ttCGTCCTGGATAGCAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGACAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGCAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAgacgac                                                                                                                                                                                                    >  1:791681/1‑250 (MQ=255)
                                                  ggattgcagaCAGATCAGACGCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAAGCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTc                                                                                                                                                                                            <  2:708408/243‑1 (MQ=255)
                                                              gATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACAAGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGTCCGCCTGAGGCAGGCCTTTAATGTTATAGGTGAAACGGCTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTAGATAGAACTAGAc                                                                                                                                                                               >  2:81405/1‑251 (MQ=255)
                                                               atcagaatcagaGTTAGTACCGGTAGTGGCCTGTACCGTAAGTTCACGCACACGATGTAAGTTGTTGTTGATTTCGGACCGCGCGCCTTCGGAGAGATGAGCAACGGAGATACCGTCGTTGGCGGTACGGGCCGCCTGAGTCAGGCCTTTACTGTTAGAGGTGACACGGTTAGCAATCGACTGACCCGCTGCGTCCTCCTTCGCGCTGTTAATACGCAACAAAGAAGACAGACGCTCGATAGacaccgacc                                                                                                                                                                              >  2:706891/1‑243 (MQ=255)
                                                               aTCAGACTCAGAGTTAGTCCCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGAGAGCCTTCGGTGGTCTACGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAACCGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTTATACGCAAGCCAGAAGACAGACGCTCgttacaaccagacc                                                                                                                                                                              >  2:106290/1‑238 (MQ=255)
                                                                 cAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAgc                                                                                                                                                                            <  2:424227/251‑1 (MQ=255)
                                                                  aGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAgcg                                                                                                                                                                           >  1:848979/1‑251 (MQ=255)
                                                                    aCTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCa                                                                                                                                                                         <  1:74386/251‑1 (MQ=255)
                                                                          agTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGCTAGAACTCGACAGCGCCGACTgg                                                                                                                                                                   >  1:353981/1‑251 (MQ=255)
                                                                             tAGTACCGGTAGTGTCCTGTCCCGTCAGTTCCCGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGGGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGtt                                                                                                                                                                 <  2:6500/250‑1 (MQ=255)
                                                                              tggcccGGGAGTGGTCTGGACCGACAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCATGCCGTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCTGACTGGTTc                                                                                                                                                                <  2:367974/246‑1 (MQ=255)
                                                                               gTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCtt                                                                                                                                                              <  1:65969/251‑1 (MQ=255)
                                                                                       aGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGGCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGatat                                                                                                                                                      <  1:732084/251‑1 (MQ=255)
                                                                                             cTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGGTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATCGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTAtttt                                                                                                                                                >  1:206302/1‑251 (MQ=255)
                                                                                                 aCCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGt                                                                                                                                            >  1:301149/1‑251 (MQ=255)
                                                                                                 aCCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGACTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGt                                                                                                                                            >  1:699511/1‑251 (MQ=255)
                                                                                                    gTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGAt                                                                                                                                         <  1:846454/251‑1 (MQ=255)
                                                                                                     tCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATc                                                                                                                                        >  1:777136/1‑251 (MQ=255)
                                                                                                       aattccctcacaCGCTGTAAGTTGTTTATGATGTCGTACAGCGTGCCTGCGGGGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACAGACCCGCTGCGTCATCCTTAGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAg                                                                                                                                      <  2:191472/243‑1 (MQ=255)
                                                                                                       aGTTCACGCCCACGCTGTAAGTTGTTGATGATTTCGGCCAGCGCGCCTCCGGTGGTCTGCGCAACGGAGCTACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGCACTCGACAGCGCAGCCTGATACTTGTTGATTTTCTTTTGaatgctc                                                                                                                                        >  2:759240/1‑243 (MQ=255)
                                                                                                        gTTCACGCCCCCGCTGTAATGTGTTTGTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGATATGCCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGc                                                                                                                                     <  2:298350/251‑1 (MQ=255)
                                                                                                        gTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGc                                                                                                                                     <  1:345601/251‑1 (MQ=255)
                                                                                                         ttcactcccacGCTGTAGGTTGTTGTAGATTTCGGACAGCGCGCCTGCGGTGGTCTGCTAAACGGAGATACCGTCGTGGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGTCTGGTTCTTGTTGATATTATTTTGAGTGATCAGc                                                                                                                                     <  2:603922/242‑1 (MQ=255)
                                                                                                             cGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCgag                                                                                                                                  >  1:759141/1‑249 (MQ=255)
                                                                                                                 cacGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCGTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTg                                                                                                                            <  1:152116/251‑1 (MQ=255)
                                                                                                                          aGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTAtt                                                                                                                    >  1:660343/1‑250 (MQ=255)
                                                                                                                             gttttttgATTTCGGACAGCGCGCCTTCGGTGGTTTGCGCAACGGGGATACCGTCGTTGTCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGTCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAAt                                                                                                                 <  2:413601/245‑1 (MQ=255)
                                                                                                                                 gttgATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACtt                                                                                                            >  1:300559/1‑251 (MQ=255)
                                                                                                                                     aTTTCGGACAGCGCGGCTTCAGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGTGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAAAGGTTAGCAATAGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCACGCCAGAAGACGGAAGCTCGATAGAACTCGACGGCTCACACTGGCTCTTGTCTATTTTATTTTGAGGGATCAGCGAGGGGCTGTTGGttttcatagattggcc                                                                                                        >  2:321987/1‑236 (MQ=255)
                                                                                                                                       ttCGGAAAGCGCTCCCTCTGGGGTCTGCGCAACGGAGATACCGTCGTAGGCGTTACGGGCCGCCTGAGTCATGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCTACTGACCCGCTGCCGCATCCTTCGCGCTGTTAATCCGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATGTTGAGTGCTCAGCGAGAGGCTGTTGGTATTAATGACTTGTGcc                                                                                                       <  2:329404/250‑1 (MQ=255)
                                                                                                                                         cGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGCTATTATTTTGAGTGATCAGCGAGCGGCTGTTGGTATTAATGACTTGTGCCa                                                                                                      >  2:767179/1‑249 (MQ=255)
                                                                                                                                         cGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCa                                                                                                      <  1:767179/249‑1 (MQ=255)
                                                                                                                                         cGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATg                                                                                                    <  1:281138/251‑1 (MQ=255)
                                                                                                                                         cGGACAGCGCGCATTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATg                                                                                                    >  1:701532/1‑251 (MQ=255)
                                                                                                                                           gACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCAGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGAAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTGTTAATGAATTGTGCCATGAt                                                                                                  >  2:325287/1‑251 (MQ=255)
                                                                                                                                            aCAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGAtt                                                                                                 <  1:241540/251‑1 (MQ=255)
                                                                                                                                             cAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTc                                                                                                >  1:853556/1‑251 (MQ=255)
                                                                                                                                                    ccTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGAc                                                                                                                                                                      <  1:511548/174‑1 (MQ=255)
                                                                                                                                                    ccTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGGTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCACGAAGACAGACGCTCGATAGAACTCGACAGCGCAGAc                                                                                                                                                                      >  2:511548/1‑174 (MQ=255)
                                                                                                                                                      gtCGGTGGGTTGCGCAACGGAGATACCGTCGTTGGCGTTACGGACCGCCTGAGTCAGGCCTTTAATGTTAGAGGATAAACGGTTAGCAATCGACTGACACGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTTATATTATTTTGTGTGATCAGCGAGAGGCTGTTGGTAGTAATGACTTGTGCAATGATTCGTTATCCta                                                                                       <  2:457323/250‑1 (MQ=255)
                                                                                                                                                       gcGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTTCGGGCCGCCTTAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATAGTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCt                                                                                        <  2:569991/248‑1 (MQ=255)
                                                                                                                                                         ggtggtCAGCGCAACGGATATACCGTCGTTGGCGTTACGGGCCTCCTGAGTCAGGCTTTTAATGTTAGAGGTGAAACGGTTATCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTatat                                                                                    <  2:264328/251‑1 (MQ=255)
                                                                                                                                                           tggtCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCg                                                                                               >  1:266724/1‑238 (MQ=255)
                                                                                                                                                           tggtCTGCGCAACGGAGATACCGTCGTGGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCg                                                                                               <  2:266724/238‑1 (MQ=255)
                                                                                                                                                                  cgcAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTg                                                                                  <  1:517159/244‑1 (MQ=255)
                                                                                                                                                                  cgcAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGGGCCATGATTCGTTATCCTATATTg                                                                                  >  2:517159/1‑244 (MQ=255)
                                                                                                                                                                      aCGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTgt                                                                       >  1:671421/1‑250 (MQ=255)
                                                                                                                                                                           gaTACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACg                                                                  <  1:175187/251‑1 (MQ=255)
                                                                                                                                                                            aTACCGTCGTTGGCGTTACGGGCCGCCTGGGTCAGGCCTTTAAGGTGAGATGTGAAACGGTTAGCAAGCTACTGACCCGATGCGTCATCCTGCGCGCTGTTAATTCGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTTGTTCTTGTTGATATTATTTTGAGTGATCAGCTAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATTCTATATTGCAAGGCGTTGATTACg                                                                  <  2:777021/250‑1 (MQ=255)
                                                                                                                                                                                cgtcgtTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTg                                                             <  2:54791/251‑1 (MQ=255)
                                                                                                                                                                                      gtgcgttttgggCCGCCTGCGGCAGGCCATTAATGTTAGAGGTGAAACGGTTAGCAAGCGACTGACCCGCTGCTTCATCCTTAGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGttt                                                        <  2:688172/241‑1 (MQ=255)
                                                                                                                                                                                       ggCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTGCTTGTTGATAGTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTTCAAGTCGTTGATTACGTATTGGGTTTcc                                                      <  1:738402/251‑1 (MQ=255)
                                                                                                                                                                                                   cgccTGAGGCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTc                                          <  2:777136/251‑1 (MQ=255)
                                                                                                                                                                                                        taGTCAGGCTTTTAATGTGAGAGGTGAAACGGTTAGCAATCTACTGACCAGCTGCGTCATCCTGCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTAGACAGCGCAGACTGGTTCTTGTTGATAGTATTTGGAGTGATCAGCGAGAGGCTGTTGGTAGTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGGCGGTGATTACGTATTGGGTTGCCACCCGTCGGCTCAAt                                       <  2:353981/248‑1 (MQ=255)
                                                                                                                                                                                                        gAGTCAGGCCTTTAATGTTAGAGGTTAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTTGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAAt                                       <  1:871818/249‑1 (MQ=255)
                                                                                                                                                                                                        gAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATc                                      <  1:46515/250‑1 (MQ=255)
                                                                                                                                                                                                        gAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCg                                     <  2:759141/251‑1 (MQ=255)
                                                                                                                                                                                                           ccAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCAGCCTTCGAGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCg                                  <  2:650835/250‑1 (MQ=255)
                                                                                                                                                                                                               gCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCaa                              <  1:428753/251‑1 (MQ=255)
                                                                                                                                                                                                                    tAATGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGTCTCAATCGCCGTCAACCCt                          <  1:516341/250‑1 (MQ=255)
                                                                                                                                                                                                                     aaTGTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCGTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTg                         <  1:72035/250‑1 (MQ=255)
                                                                                                                                                                                                                        gTTAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCATACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGAGTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTAt                     <  2:393865/251‑1 (MQ=255)
                                                                                                                                                                                                                          tAGAGGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATc                    <  2:301149/250‑1 (MQ=255)
                                                                                                                                                                                                                          tAGAGGGGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATc                    <  2:198879/250‑1 (MQ=255)
                                                                                                                                                                                                                          gagagGTGAAACGGTTAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCTCCCGTCGGCTCAATCGCCGTCAACCCTGTTATc                    <  1:706891/249‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcttcgCAATCGACGGACCCGCTGAGTACTCCTTCGCGCCTTTAATACGAAAGCCAGAAGACCGACGCTCTATAGAACTCGACAGCGCTGTCTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATTCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCAGTCGGCTCAATCGCCGTCAACCCGGTTACCGTCTGTCGTaaa        <  2:105918/245‑1 (MQ=255)
                                                                                                                                                                                                                                      gtttAGCAATCGACTGCACCTCTGCGTCATCCTTCGCGCTGTTAATTCGCAAGCCAGAAGACATACGCTTGATAGAACTCGACAGCGCAGACTGGTTCTTGTGGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTTCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCTTaaaa       <  2:659104/249‑1 (MQ=255)
                                                                                                                                                                                                                                      ggTTAGAAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTaaaa       <  1:325287/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAaaca     <  1:395720/250‑1 (MQ=255)
                                                                                                                                                                                                                                         tAGCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAaaca     >  1:506028/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCAATCGACTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTAGATAGAACTCGACAGCGAAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTTCAATGATTCGTTATCCTATATTGCAAGTCGTTGATTCCGAATGTGGTTTCCACCCGTAGGCTCAATCGCCGTCTACCCTGTTATAGGCAGTCGTaaaaatacc  >  2:765276/1‑246 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CTGACCAGATACGCGGTCAATTTCATCCAGACGGGATTTAATTTCGTCCTGGATAGAAGACAGATCAGACTCAGAGTTAGTACCGGTAGTGGCCTGTACCGTCAGTTCACGCACACGCTGTAAGTTGTTGTTGATTTCGGACAGCGCGCCTTCGGTGGTCTGCGCAACGGAGATACCGTCGTTGGCGTTACGGGCCGCCTGAGTCAGGCCTTTAATGTTAGAGGTGAAACGGTTAGCAATCGCCTGACCCGCTGCGTCATCCTTCGCGCTGTTAATACGCAAGCCAGAAGACAGACGCTCGATAGAACTCGACAGCGCAGACTGGTTCTTGTTGATATTATTTTGAGTGATCAGCGAGAGGCTGTTGGTATTAATGACTTGTGCCATGATTCGTTATCCTATATTGCAAGTCGTTGATTACGTATTGGGTTTCCACCCGTCGGCTCAATCGCCGTCAACCCTGTTATCGTCTGTCGTAAAACAACC  >  NC_000913/2003220‑2003705

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