Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,518 T→C E11E (GAA→GAG)  araD ← L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,5180TC98.6% 438.5 / ‑4.7 139E11E (GAA→GAG) araDL‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (70/66);  total (71/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 139 total aligned reads displayed.

CGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAC  >  NC_000913/66277‑66763
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
cgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATAcc                                                                                                                                                                                                                                              >  1:341194/1‑251 (MQ=255)
cgccAATGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATAcc                                                                                                                                                                                                                                              >  2:8746/1‑251 (MQ=255)
  ccAATGGACGGGAATGACTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTg                                                                                                                                                                                                                                            <  1:615056/251‑1 (MQ=255)
   cgaTGGAAGGGAAGTCCTGATAGAGCAGTAGGTGAGTTGGCGTATCGGAGGAGGGCCTTTTCGCACCGTCAACCACTTCAACGGTTTCGATGCTAACCACGACCATAGCGTCAGCGGGCATGATGCTGTAATCGACGCCGGAAGGGTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTGTTGGCAGCGCCAGGTTGGCCTCTAATACCTga                                                                                                                                                                                                                                           <  2:41962/249‑2 (MQ=255)
    aaTGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGATTTTTCGCACCTTCAACCACTTCACCGGGTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACAAGGTTATGTTTTGGCAGCGCCAGGGTGGCCTCTCATAACTGGc                                                                                                                                                                                                                                          >  2:791142/1‑251 (MQ=255)
     aTGGACGGGAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCg                                                                                                                                                                                                                                         >  1:229833/1‑251 (MQ=255)
       ggacgggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGtt                                                                                                                                                                                                                                       <  1:426067/247‑1 (MQ=255)
         acgtgAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTg                                                                                                                                                                                                                                     <  1:216686/247‑1 (MQ=255)
          ggggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTga                                                                                                                                                                                                                                    <  2:238984/251‑1 (MQ=255)
            ggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATTCCTGGCGTTTgaga                                                                                                                                                                                                                                  <  1:229970/251‑1 (MQ=255)
            ggAATGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTgaga                                                                                                                                                                                                                                  <  1:465549/251‑1 (MQ=255)
               aTGCCTGATAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGGGAtct                                                                                                                                                                                                                               >  1:327825/1‑251 (MQ=255)
                 gCCTGATAGAGCCGTCGGTGAGTTGGCGTATCGGAGGCGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGATTGATCACAAAGACGCAGCACTCGCGATCAACGACGCTGACGTTGCCCCATGTTAGCGTGACCAGGTTATGTATTGGCAGCGCCAGGTTTGCCTCTACTACCTGGCGTTTGAGAtcttc                                                                                                                                                                                                                             >  2:741446/1‑251 (MQ=255)
                  ccTGATAGAGCAGTCGGTGAGTTGGCGTATCGGTGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGGTGGCCTCTAATACCTGGCGTTTGAGAtcttct                                                                                                                                                                                                                            >  1:273283/1‑251 (MQ=255)
                  catGATAGAGCAGTCGGGGAGGTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCACCCACTTCACTGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGATGTAATCGACGCCGGAAGGTTTGATCCCAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGAtcttct                                                                                                                                                                                                                            <  2:182277/249‑1 (MQ=255)
                   cTGATAGAGCAGTCGGAGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCAATGTGAGCGTGACCCGGTTATGTTTTGGCAGGGCCAGGTTGGACTCTAATACCTGGCGTTTGAGAtcttct                                                                                                                                                                                                                            >  2:459363/1‑250 (MQ=255)
                      aTAGAGCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATAAACGGCGCTGACGTTGCCCCATGTGAGCGTGCCCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTCATACCTGGCGTTTGAGATCTTCTAAc                                                                                                                                                                                                                         >  2:746966/1‑250 (MQ=255)
                        agagCAGTCGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATg                                                                                                                                                                                                                      <  1:715992/251‑1 (MQ=255)
                            cagtcGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCCGGTTGGCCTCTAATACCTGGCGGTTGAGATCTTCTCACATGTTg                                                                                                                                                                                                                   >  2:642175/5‑250 (MQ=255)
                             agtcGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGAc                                                                                                                                                                                                                 >  1:420781/4‑251 (MQ=255)
                             agtcGGTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGAc                                                                                                                                                                                                                 >  1:440771/4‑251 (MQ=255)
                               acGGTGCGTTGGCGTATCGGAGGAGGGCTTTGTCGCAACTTCAACCAATTCACCGGTTTAGATGCTAAACACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGCCGCCGCGCTCGCGATAAACGGCGATGACGTTGCCCCATGTGAGCGTGCCCAGGTTAGGTTTTGGCCGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGTTCTTCTAACATGTTGACTc                                                                                                                                                                                                               >  2:883853/2‑251 (MQ=255)
                                gggTGAGGTGGCGGATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTcc                                                                                                                                                                                                              <  2:383630/250‑1 (MQ=255)
                                  gTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGAACAGGTTATGATTTGGCAGCGCCAGGTTGCCCTCTATTCACGGGCGTTTGAGATCTTTTAACATGTTAACTCCtt                                                                                                                                                                                                            >  2:269609/1‑251 (MQ=255)
                                  gTGAGTTGGCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTGCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCACTGTGAGCGTGACCGGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGCGATCTTCTAACATGTTAAATCCtt                                                                                                                                                                                                            >  2:177175/1‑251 (MQ=255)
                                  gTGAGTGGGCGTATCTGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCtt                                                                                                                                                                                                            <  2:919026/251‑1 (MQ=255)
                                         ggCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGGGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCg                                                                                                                                                                                                     >  2:970517/1‑251 (MQ=255)
                                          gCGTATCGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCTTGCCgg                                                                                                                                                                                                    >  1:500272/1‑251 (MQ=255)
                                                cGGAGGAGGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCATGCCGGATgcgc                                                                                                                                                                                              <  1:446990/251‑1 (MQ=255)
                                                  gaggagGGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGCCCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCtt                                                                                                                                                                                            >  1:984478/1‑251 (MQ=255)
                                                   aggaggGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGTTCATGATTCTGTAATCGACGCCGGAAGGTTTGATCTCAAAGACGCCGCGCTCTCGATCAACGGCGCTGACGTTGCCCCAGGTGAGCGTGACCAGGTTATGGTGTGGCAGCTCCAGGTTGGCCTTTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCttt                                                                                                                                                                                           <  1:906768/251‑1 (MQ=255)
                                                    ggaggGCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAAACACGACCATATCGTCAGCGGTACTGATGCTGTAAGCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACAAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACATGGCGTTTGAGATCTTCTTACATGTTTACTCCTTCGTACCGGATGCGCTTTg                                                                                                                                                                                          >  2:621931/1‑251 (MQ=255)
                                                       gggCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCtt                                                                                                                                                                                       <  1:295607/251‑1 (MQ=255)
                                                         gCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTAt                                                                                                                                                                                     >  1:524771/1‑251 (MQ=255)
                                                         gCTTTTTCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTAt                                                                                                                                                                                     >  1:632408/1‑251 (MQ=255)
                                                            ttttCGCACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATcc                                                                                                                                                                                   <  1:29751/250‑1 (MQ=255)
                                                                       tCAACCGCTTCACCGTTTTGGAGGCTAACCCCGACCATATCGTCAGCGGGCATGATGCTGGAATCGACGCCCGACGGTTTGATCACAAAGTCGCCGCGCTCGCGTGCATCGGCGCTGCCGTTGCACCATGTGAGCGTTACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAAGACCTGGCGTTGGAGATCTTCTAACATGTTGACTCCTACGTGCCTGATGCGCTTTGCTTATCCGGCCTACaaaa                                                                                                                                                                        <  2:380241/250‑1 (MQ=255)
                                                                        cAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTTACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCTGATGCGCTTTGCTTATCCGGCCTACAAAAt                                                                                                                                                                       <  1:888156/250‑1 (MQ=255)
                                                                             aCTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGc                                                                                                                                                                 <  2:558049/251‑1 (MQ=255)
                                                                             aCTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGAAAGGATTGATCACACAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGGTATGTGTTGGCAGCGCCAGGTTGGCCTCTAAACCCAGGCGGTTGAGATCATCTAAAATGTTGACTCCTTCGTGCCGGATACTCTTTGCTGATACGGCCTCAAAAATCGccgc                                                                                                                                                                 >  2:620176/1‑248 (MQ=255)
                                                                                    cGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAgg                                                                                                                                                          <  2:813487/251‑1 (MQ=255)
                                                                                     ggTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGc                                                                                                                                                         >  1:770339/1‑251 (MQ=255)
                                                                                         tCGATGCTAACCACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTg                                                                                                                                                      <  1:575086/250‑1 (MQ=255)
                                                                                          cGATGCTAACCCCGTCCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATGCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGa                                                                                                                                                     <  2:315787/250‑1 (MQ=255)
                                                                                            aTGCTAACAACAACCATATCGTCAGCGGCCATGATGCTGTAATCGACGCCGGAAGGTATGATCACAAAGACGCCGCGCTCGCGATCACCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGCTATGTTTTGGCAGCGCCAGGTTGGCCTCTCATACCTGGCGTTTGAGATCTTCTAACCTGTTGACTCCTTCGAGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGTCCTgtta                                                                                                                                                   >  2:480321/1‑247 (MQ=255)
                                                                                                   ccACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACgcgc                                                                                                                                           >  1:971348/1‑251 (MQ=255)
                                                                                                   ccACGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACgcgc                                                                                                                                           >  1:108475/1‑251 (MQ=255)
                                                                                                     aCGACCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCa                                                                                                                                         <  1:658599/251‑1 (MQ=255)
                                                                                                        aCCATATCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCg                                                                                                                                      >  1:406851/1‑251 (MQ=255)
                                                                                                             atCGTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCa                                                                                                                                 <  1:118082/251‑1 (MQ=255)
                                                                                                               cGTCAGCGGTCATGATGCTGTACTCGACGCCGGAAGGTTTGATCACAAAGCCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTTAGCGTGACCTGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGAATGCGCTTTGCTTATCCGGCCTACAATATCGCAGCGTGTAGGccgtcttagccgctccctcctcgcatc                                                                                                                               >  2:36119/1‑226 (MQ=255)
                                                                                                                gTCAGCGGTCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATTCCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa                                                                                                                              <  1:411787/251‑1 (MQ=255)
                                                                                                                gTCAGCGGTCATGATGATGTAATCGACGCCGGAAGGTTTGTTCACAAAGACGCCGCGCTCGCGCTCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGGTCTTCTACCCTGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCAGGCCTACAAAATTGCATCGTGTAGGCCTGATAAGCAGCGCCAGCggcacatcc                                                                                                                              >  2:438982/1‑243 (MQ=255)
                                                                                                                   aGCGGTCATGATGCTGTAAGCGACGCCGGAATGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGGTGGCCTCTAATACCTGGCGTTTGAGATCTGCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGc                                                                                                                           <  1:431300/251‑1 (MQ=38)
                                                                                                                        tCATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTg                                                                                                                       >  1:1233/1‑250 (MQ=255)
                                                                                                                         caggatgctttaaacGACGTCGGAAGGTTTGTGCACAAAGAGGCCGCGCTCGCGTTCTACGGCTCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCGGGCGTTTGAGATCTTCTAACAGGTTGACTCCTTCGTGCCGGATGCTCTTTGCTGATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCTTCAGGCGTTGa                                                                                                                      <  2:515798/236‑1 (MQ=25)
                                                                                                                         cATGATGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGa                                                                                                                      >  1:704827/1‑250 (MQ=255)
                                                                                                                           tgatgCTGTAATCGACGCCGGAAGGTTTAATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATTCGGCCTACAACATCGAAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATg                                                                                                                   >  2:930553/5‑251 (MQ=38)
                                                                                                                           ggatgCTGTAATCGACGTCGGAAGGTTTGATAACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATg                                                                                                                   <  2:753752/247‑1 (MQ=255)
                                                                                                                            gatgttGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTTACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGc                                                                                                                  <  2:553181/246‑1 (MQ=255)
                                                                                                                             atgCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCAGGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGcc                                                                                                                 <  2:229833/249‑1 (MQ=255)
                                                                                                                             atgCTGTAATCGACGCCGGAAGGATTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGcc                                                                                                                 >  1:143570/3‑251 (MQ=255)
                                                                                                                                 tGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAt                                                                                                             >  1:810264/1‑251 (MQ=255)
                                                                                                                                     aTCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATgcgc                                                                                                         >  1:519981/1‑251 (MQ=255)
                                                                                                                                        gACGCCGGAAGGGTTGATCACAAAGTCGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCttt                                                                                                      >  1:319158/1‑251 (MQ=255)
                                                                                                                                         aCGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTg                                                                                                     <  1:292541/251‑1 (MQ=255)
                                                                                                                                          cGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGc                                                                                                    <  1:867018/251‑1 (MQ=255)
                                                                                                                                                 aGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGAATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                             >  1:205031/1‑251 (MQ=255)
                                                                                                                                                  ggTTTGATCACAAAGGCGCCGCGCTCGCGATCAACGTCGCTGACGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCTTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCgg                                                                                            <  1:462229/251‑1 (MQ=255)
                                                                                                                                                   gTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGc                                                                                           >  1:167406/1‑251 (MQ=255)
                                                                                                                                                      tgagaccaAAGCCGCCGCTCTCGCGATCAACGGCGCTGACGGTGCCCCATGTGAGCGTGACCAGGTTATTGTTTGGCACCGCCAGGTTGGCCTTTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACGCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGGTGAATGCCGGATGCGCTTTGCTTATCCGGCCt                                                                                         <  2:500272/244‑1 (MQ=21)
                                                                                                                                                          cacaAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACaaa                                                                                    <  1:841327/251‑1 (MQ=255)
                                                                                                                                                            caAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAt                                                                                  >  1:753943/1‑251 (MQ=255)
                                                                                                                                                             aaaGACGCCGCGCTCGCGATCAACGGCGCTGAAGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCCGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGACGTATGCGCTTTGCTTATCCAGCCCTCAAAATc                                                                                 >  2:172253/1‑251 (MQ=37)
                                                                                                                                                               aGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTGTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCgg                                                                               <  2:71857/251‑2 (MQ=255)
                                                                                                                                                                     cgcgCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCggaacgcg                                                                         >  1:809136/1‑244 (MQ=255)
                                                                                                                                                                        gcTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAg                                                                      <  1:838572/251‑1 (MQ=255)
                                                                                                                                                                          tCGCGATAAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGc                                                                    >  2:137351/1‑251 (MQ=255)
                                                                                                                                                                            gcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCt                                                                  >  1:512221/1‑251 (MQ=255)
                                                                                                                                                                             cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTTACCAGGGTATGTGTTGGCAGCGCCAGGTGGGCCTCTAATACCTGGCGTTTGAGATATTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTg                                                                 <  2:235722/251‑1 (MQ=21)
                                                                                                                                                                             cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTg                                                                 >  1:104429/1‑251 (MQ=255)
                                                                                                                                                                              gATCAACGGCGCTGCCGTTGCCCCAGGTGAGGGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCGCTAATACCGGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGa                                                                <  2:434521/251‑1 (MQ=21)
                                                                                                                                                                                  aaCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAg                                                            <  1:534189/251‑1 (MQ=255)
                                                                                                                                                                                     ggCGCTGACTTTGACCCATGTGAGCGGGACCAGGTTGTGTTTTGGCAGGGCCAGGTTGGCCTCTAATACCTGGCTGTTGAGATCTTCTAACATGTTGACGCCTTCTTGCCGGATGCGCGTTGCTTATCAGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGAcg                                                         <  2:418529/251‑1 (MQ=21)
                                                                                                                                                                                        gcTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGCGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACgcgc                                                      >  2:575701/1‑251 (MQ=255)
                                                                                                                                                                                           gACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGATTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGGTGCATGCCGCATGCGCTTTGCATATCCGGCCTACAAAATCGGAACACGTAGGCCTGATAAGACGCGCCa                                                    >  2:652328/1‑250 (MQ=25)
                                                                                                                                                                                             cGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACACAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGTCGTGGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGACCGCGTAGGCCTGATACGAcgctcacgat                                                 >  2:135576/1‑244 (MQ=25)
                                                                                                                                                                                              gTTGCCCCATGTGAGCGTGACCAGGTTCTGTTTTGGCTGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGt                                                >  1:771602/1‑251 (MQ=255)
                                                                                                                                                                                                  cccaaTGTGAGCGTGCCCAGGTTATGTTTTGGCAGAGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGCCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAACATCGGAACGCGTCGGCCTGATAAGACGCGCCAGCTTCGc                                             >  2:737221/5‑250 (MQ=21)
                                                                                                                                                                                                    ccATGTGAGCGTGACCAGGTTAGGTTTTGGCAGCGCCCGGTTGGCCTCTAATTCTGGGGGGTTGATATCTTCTAACATGTTGTCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCGGCGGGTAGGCCTGATAAGACGCTCCATCGTCGCATCAGGGGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATc                                          <  2:469732/251‑1 (MQ=18)
                                                                                                                                                                                                     cATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa                                         <  1:642175/251‑1 (MQ=255)
                                                                                                                                                                                                           ggCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGtt                                   <  1:746966/250‑1 (MQ=255)
                                                                                                                                                                                                                      gttatgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGACCGCGTAGGCCTGATAAGACGCTCCAGCGACGCATCACGCGTTGAATGCCgg                          >  2:662150/5‑249 (MQ=21)
                                                                                                                                                                                                                       ttGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACGTGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAA‑TCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGATGAATGCCGGATgc                      >  2:496030/1‑251 (MQ=255)
                                                                                                                                                                                                                        tatgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATgc                      >  1:387439/3‑251 (MQ=255)
                                                                                                                                                                                                                         atgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATAAGGCGTTGAATGCCGGATgcg                     >  1:988079/2‑251 (MQ=255)
                                                                                                                                                                                                                                ggCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGACTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCAGGATGCGCTTTGCt              >  1:571900/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTAt           >  1:513033/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg        >  1:342012/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcCAGGTTGGCCTCTAAGTCCTGGCGTTTGAGATCTTCTAACATGTTGACCCCTTCGTGCCGGACGCGCTTTGCTTACCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGTCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCg        <  1:620176/251‑1 (MQ=21)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
CGCCAATGGAGGGGAATGCCTGATAGAGCAGCCGGTGAGTTGGCGTGTCGGAGGAGGGCTTTTTCGTACCTTCAACCACTTCACCGGTTTCGATGCTAACCACGACCATATCGTCAGCGGTCATGACGCTGTAATCGACGCCGGAAGGTTTGATCACAAAGACGCCGCGCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAC  >  NC_000913/66277‑66763

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: ‑