Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 66,847 C→T G497G (GGG→GGA)  araA ← L‑arabinose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091366,8470CT97.5% 306.9 / NA 123G497G (GGG→GGA) araAL‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  major base T (56/61);  minor base G (0/3);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTT  >  NC_000913/66603‑67097
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTaa                                                                                                                                                                                                                                                        <  2:442348/251‑1 (MQ=18)
gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTa                                                                                                                                                                                                                                                         >  1:271955/1‑250 (MQ=18)
               aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                                                                                          <  1:524073/250‑1 (MQ=18)
               aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                                                                                                                                                                                                          <  1:970517/250‑1 (MQ=18)
                        cAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGcgcg                                                                                                                                                                                                                               >  1:551457/1‑251 (MQ=18)
                           cGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCt                                                                                                                                                                                                                            >  1:987664/1‑251 (MQ=18)
                               gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCGTTGCTTATCCGGCCTACCCAACCGCAGCGTGTAGGCGGGATGAGCCGCGTCAGCGTCGCATCAGGCGTTAACCACATGCTTCGACTACTTAGCGCAGAAATCCGTAATACACTTCGGTCCC‑GCGCAGcgcgggttta                                                                                                                                                                                                                         >  2:5616/1‑244 (MQ=9)
                               gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTa                                                                                                                                                                                                                         >  1:475848/1‑250 (MQ=18)
                                   cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCTTGTAGGCCTGCTCAGACTCTTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACATATCGACGACATCCGTAATACCCTTCGTTCCA‑GCGCAGCGCGGCTTTACACGc                                                                                                                                                                                                                    >  2:6759/1‑251 (MQ=17)
                                    aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCt                                                                                                                                                                                                                   >  1:439771/1‑251 (MQ=255)
                                     ggCGTTGAATGTCGGCTGCGCTTTGCTTACCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCTCATCAGGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTg                                                                                                                                                                                                                  <  2:135323/251‑1 (MQ=18)
                                       cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAACCGCGGGc                                                                                                                                                                                                                >  1:954197/1‑251 (MQ=18)
                                           gAATGCCGGTAGCTCTTTGCTTATCCGGCCAACCAAAACGGAACGCGGAGGCCTGAAAAGACGCGCCGGCGTCGCATCAGGTGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAAGCCGTAATACACTGCGTTCCA‑GCGCAGCGCGTTTTTAAACGCTGGCAGGc                                                                                                                                                                                                            <  2:571900/251‑1 (MQ=18)
                                           gAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAgg                                                                                                                                                                                                             <  2:799294/250‑1 (MQ=18)
                                                  ggATGCGCTTTGCTTATCCGGCCTACAAGATCGGACCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCATTACACCCCTTATGCGACTCCTTAGCGACACAATCCGAAATACGCTTGGTTCCA‑GCGCCTGGCGTCTTTAAACGCTAGCAGGCgtggtgc                                                                                                                                                                                                     >  2:475929/1‑247 (MQ=9)
                                                   gATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACCCCCCTGATGCGACTACTTAGCGTCGAAATCCGTAATAAACTTCGTTCCA‑GCGCACCGCGTCTTTAACCGCTGGCAGGCggttatcg                                                                                                                                                                                                    >  2:223981/1‑244 (MQ=17)
                                                         cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCA‑GAGCAGCGCGTCTTTAAACGCTGGCAGGCgggtattcttatc                                                                                                                                                                                               >  1:537598/1‑238 (MQ=18)
                                                         cTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTatc                                                                                                                                                                                               >  1:372730/1‑250 (MQ=18)
                                                           ttGCTGATCCGGCCTACTAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTGCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatca                                                                                                                                                                                              <  2:143570/249‑1 (MQ=18)
                                                             gCTTAGCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca                                                                                                                                                                                          <  1:600919/251‑1 (MQ=21)
                                                              cTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTCCACACCTGATGCGCCTACTTAGCGACGACATCCGTAATACCCTTCGTTCCA‑GCGCAGCGCGACTTTAAACGCTGGCCGGCggggatcgttgtcactcac                                                                                                                                                                                         >  2:691235/1‑233 (MQ=17)
                                                               ttATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCAcc                                                                                                                                                                                        >  1:158268/1‑251 (MQ=18)
                                                                tcaccgtccagcAAAATCGGAACGCGTAGGCCTGCTTAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTGAGCGACGAAATCCGTAATTCACTTCGTTCCA‑GCGCCGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGTTATCAATCAc                                                                                                                                                                                         <  2:67473/238‑1 (MQ=18)
                                                                     ggCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGCCTAATTAGCGACGAAATCCGTAATAAACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCggggatcggttttaatcaacctggtt                                                                                                                                                                                  >  2:878655/1‑226 (MQ=18)
                                                                      cccTACAAAATCGTTACGCGTAGGCCTGATAATGCGCGGCGGCGTCGCAGCAGGCGTTGAATGCCGGATGCGATTTGCTTACACGGGCTACAAAATCGCATCGTGTAGGCCTGATAAGACGCGTCAGCGTAGCATCAGACGTTACACACCTTAGGCGTCTATTTAGCGACGTAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCTCCGTTAttt                                                                                                                                                                                 <  2:835880/250‑1 (MQ=255)
                                                                                     acgcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCAtct                                                                                                                                                                  <  1:738691/250‑1 (MQ=25)
                                                                                       gcgTAGGCCTGATAAGACGCGCCAGCTTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCg                                                                                                                                                                <  1:215389/251‑1 (MQ=25)
                                                                                       gcgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCg                                                                                                                                                                <  1:172253/251‑1 (MQ=25)
                                                                                           aGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGGGCCTCTCGGCg                                                                                                                                                             >  1:611681/1‑250 (MQ=25)
                                                                                                 gtttaGACGAGCTAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGGTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgc                                                                                                                                                      <  2:704827/247‑1 (MQ=25)
                                                                                                 gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTACTACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCATTCACCGTGATTTCAATGTCGTGCATCTCGGTGAACTGgc                                                                                                                                                      >  1:515350/1‑251 (MQ=25)
                                                                                                  aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                                                                                     <  1:47930/251‑1 (MQ=25)
                                                                                                  aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                                                                                     <  1:621931/251‑1 (MQ=25)
                                                                                                           gcCAGCGTCGCATCAGGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCCGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTACTTTTCATGGCGTGCATATCGGCGGAATGGCGCATCTCGtt                                                                                                                                            >  1:756626/1‑251 (MQ=18)
                                                                                                           gcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGtt                                                                                                                                            <  1:693596/251‑1 (MQ=37)
                                                                                                              ctcttcgcctcctgCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCCATCCCCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCCTTggg                                                                                                                                         >  1:168354/14‑249 (MQ=25)
                                                                                                              aGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGGTCATCTCGGCGAACTGGCGCATATccttt                                                                                                                                           >  1:300319/1‑245 (MQ=25)
                                                                                                                cGTCGCATCTGGCGTTGAAGTCCGGATGCGCTTTGCGTTTCCGGTCTACCAAAGCGCAGCGTGTAGGGCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACATCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCTCAGAGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGTGATTTCAATGTCGGGCATCTCGGCGAACTGGCGCATATCGTTGAGGt                                                                                                                                       <  2:208428/251‑1 (MQ=255)
                                                                                                                 gTCGCATCAGGCGTTGAATGCCGGATGCGCTTTTCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCCAATCGTTGAGGtt                                                                                                                                      >  1:217463/1‑251 (MQ=25)
                                                                                                                     cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATAAATCACAGTGGTTTCAATGTAGTGCATCTCGGCGAAATGGGGCATATCGGTGAGGttaagg                                                                                                                                  >  1:844684/1‑247 (MQ=21)
                                                                                                                     cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCACCGGGATTTTTATATCGTGCATCTCGGCCGACTGGCGCATATTGTTGTGGttaggc                                                                                                                                  >  2:175604/1‑247 (MQ=21)
                                                                                                                        cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCa                                                                                                                               <  1:929328/251‑1 (MQ=37)
                                                                                                                        cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGGGATTTCAATGTCGTGGATCTCCGCGAACTGGCGCATATCGTTGAGGTTCAgcgc                                                                                                                                >  1:429794/1‑247 (MQ=25)
                                                                                                                             gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCCATCACCGTGATTTCAATGTCGTACATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcgccttgct                                                                                                                          >  1:599797/1‑242 (MQ=37)
                                                                                                                               ttgaggCCGGATGCGCTTTGCGTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCATGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTg                                                                                                                         <  2:440771/245‑1 (MQ=37)
                                                                                                                                 aaTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGTTAGTTATAAAACACCGGGGTTTCAATGTTGTGTATCTCGGCGCACTGGCGCCTATCCTTTAGGTGACGTGCATGGccgcc                                                                                                                       >  2:601392/1‑246 (MQ=16)
                                                                                                                                  aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTGa                                                                                                                        >  1:461821/1‑248 (MQ=37)
                                                                                                                                    gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACg                                                                                                                   <  2:313327/251‑1 (MQ=38)
                                                                                                                                      cggatggtcTTTGATTATCCGGCCTACAAAAACGCAGCGGGTAGGCCCGATAAGACGCGGCAGCGTCGCATCAGTCGTTACACATCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAAGCACCGTGATTTTAATGTCGTGCATCTCGGCGAACTGTCGCATATCGTTGAGGTTCAGCGCATGGCTGAATACGGt                                                                                                                 <  2:687803/243‑1 (MQ=255)
                                                                                                                                       ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTa                                                                                                                <  1:173395/251‑1 (MQ=38)
                                                                                                                                         attagcTTTGGTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATTCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATg                                                                                                              <  2:327825/247‑1 (MQ=37)
                                                                                                                                           gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATgg                                                                                                             <  1:137351/250‑1 (MQ=38)
                                                                                                                                            cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTg                                                                                                           <  1:218024/251‑1 (MQ=38)
                                                                                                                                             gcTTGGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTgc                                                                                                          <  2:372730/251‑1 (MQ=38)
                                                                                                                                                     tATCCGGCCTACAAAATCTCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGcca                                                                                                    <  1:77961/249‑1 (MQ=255)
                                                                                                                                                      aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGccacca                                                                                                 <  2:810264/251‑1 (MQ=255)
                                                                                                                                                          ggCCTAAAAAATCGTAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTGACACACCTGATGCGACTACTCAGCGACGAATTCCGTAATACACTTCGTTCCA‑GCGCATCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGa                                                                                             <  2:285197/251‑1 (MQ=255)
                                                                                                                                                           gCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTTTCAATCACCGTGATTTCAAATTCGTGCATCTCGGCGAACTGGGGCGTAACGTTGAGGTTCCACGCACGGCTTGAGACGGTATGGTGCGCGccccccgaggg                                                                                            >  2:595543/1‑242 (MQ=255)
                                                                                                                                                            ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCTCATATCGTTGAGGGTCCGAGCATGGCTGCAGACGGTATGGTGCGCGCCCCCCGCGAgg                                                                                           >  2:347618/1‑251 (MQ=255)
                                                                                                                                                               aCAAAAGCGGAGCGTGTAGGCCTGATAAGCCGCGTCAGCAGAGCAGCAGGCGTTACACAACGGAGGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTGCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGAATATCGTTGAGGTTCAGCGCAGGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATc                                                                                        <  2:319158/251‑1 (MQ=255)
                                                                                                                                                                 aaaaACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCa                                                                                      <  2:135556/251‑1 (MQ=255)
                                                                                                                                                                  aaaTCGCCGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCAACAGGCGTTACGCACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑TCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGGTGTCAATCACCGTGATTGCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGATGATCCAc                                                                                     <  2:205031/251‑1 (MQ=255)
                                                                                                                                                                  aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGAGATTTATATGTCGGGCATCTCTGCGAACTGGCGCCTATCGTTGGGGGTCAGCGCATGGCTTGAGACGGGATGGCGCGCGccccccgcgcgagtacac                                                                                     >  2:758493/1‑235 (MQ=255)
                                                                                                                                                                      cGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACCCTTCGTTCCA‑GCGCATCGCGTCTATAAACGCTGGCAGGCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGATGTTCAGCGCATGGCTGAAGACGGTTTGGTGCGCGCCACCAGCGAGGATCCACGCtt                                                                                 <  1:584342/251‑1 (MQ=255)
                                                                                                                                                                       gCAGCGTGTAGGCCTGATAATACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCCCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc                                                                                <  1:911505/251‑1 (MQ=255)
                                                                                                                                                                         aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAGCGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCTCGCTTCgg                                                                              <  2:708176/251‑1 (MQ=255)
                                                                                                                                                                          gCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAAACACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCgg                                                                              >  1:449839/1‑250 (MQ=255)
                                                                                                                                                                                aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGt                                                                      <  2:632408/251‑1 (MQ=255)
                                                                                                                                                                                  gcctgATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTg                                                                    >  1:504026/5‑251 (MQ=255)
                                                                                                                                                                                   cctgATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATAGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTgg                                                                   <  1:241907/248‑1 (MQ=255)
                                                                                                                                                                                       aTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGCACTGGCGCATATCGTTGAGGTTAAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCCGCGAGGATCCACGCTTCGGAA‑GCAGTTGGc                                                                  >  1:629216/1‑248 (MQ=255)
                                                                                                                                                                                       aTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGc                                                                  >  1:702692/1‑248 (MQ=255)
                                                                                                                                                                                               cgTCAGCGTCGCATCAGGCGGTACTCATCTGATTCGTCTTCTTAGCGGCGAAAGCCGTAATTCACTTCGTTCCA‑GCTCATCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTTATTTCAATGTCGTGCATCTCGGCGAACTGGCGCGTATCGTTGAGGTTCAGCTCATGGCTGAAGACGGTTTGGTGCGCGCCACCCGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTg                                                       <  2:353220/251‑1 (MQ=255)
                                                                                                                                                                                                 tCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCCATTTCGTGCCTCTCGGCGGACTGGCGCCTATCGTTGAGGTTCAGCGCATGGCTGAAGGCGGTATGGTGCGGGCCACCAGCGGGGATCCCCCCTTCGGAA‑GCAGTTGGCCGATCCGGTTgtg                                                     >  2:520113/1‑249 (MQ=255)
                                                                                                                                                                                                   aGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCt                                                   <  1:269609/251‑1 (MQ=255)
                                                                                                                                                                                                      gTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGTGCTTTCAATGTCGTGCATCTCGGGGAACTGGCGCATATCGTTGAGGTTCCGCGCATGGCTGAAGACGGTATGGTTCGCGACACACGCGAGGATCCACGCTTCGGAA‑GAAGTTGGCAGATCCGGTTGCGCtttt                                                >  1:645861/1‑250 (MQ=255)
                                                                                                                                                                                                         gCATTAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCcac                                             <  1:717794/251‑1 (MQ=255)
                                                                                                                                                                                                         gCATCAGGCGTTACACACCTGAGGCGATTAATTAGCGACGAAATCCGTAATACCCTTCGTTCCC‑GCGCAGCGAGGCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCTCCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGTCAGATCCGGTTGCGCTTTCcac                                             <  2:108475/251‑1 (MQ=255)
                                                                                                                                                                                                           aTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTTTCAATCACCGCGGTTTCAATGTCGTGCCTCTCGGCGAAATGGCGCCTATCGCTTAGGTTTAGCGCATGGCTGAAGACGGTAAGGGGCGCGCCCCCAGCGAGGATCCACCCTTCGGAA‑GAAGTTGGCAGATCCCGATGCGCTTTcccccc                                           >  2:200542/1‑247 (MQ=255)
                                                                                                                                                                                                                gCGTTACACCCCTGATGCGTCCACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                      <  2:519981/251‑1 (MQ=255)
                                                                                                                                                                                                                gCGTTACACACCTGATGCGACTTCTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGCCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTTAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTTGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                      <  2:104429/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAtt                                     <  2:167406/251‑1 (MQ=255)
                                                                                                                                                                                                                  gTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGGATTGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTc                                    >  1:868891/1‑251 (MQ=255)
                                                                                                                                                                                                                   ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTTTCGTTATCAATCACCGTGATTTCAATGTCTTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCCACGCCTGGCTGAAGACGGTATGGTGCGCGCCCCCCGCGGGGATCCACGCTTCGGAA‑GCCGTTTGCAGATTCGGTTGCGCCTTCCACAGCGCATTct                                   >  1:349276/1‑250 (MQ=255)
                                                                                                                                                                                                                     acacacctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGcc                                 <  1:921060/243‑1 (MQ=255)
                                                                                                                                                                                                                             gATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCACCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGTATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                         <  2:406851/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggcgacttactAGCGACGAAATCCCTAATACACTTTGGTCCAAGC‑CAGGGCGGCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTGCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGtt                       <  2:954197/241‑1 (MQ=255)
                                                                                                                                                                                                                                 cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGTGGTTCAGCGCATGGCTGGAGGCGGGATGGTGCGCGCCCCCAGCGAGGATCCACGCCTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCCTTTGCCACCGGCAGtttt                     >  2:255382/4‑250 (MQ=255)
                                                                                                                                                                                                                                 cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGTAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGGAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAggttt                     >  1:704961/4‑247 (MQ=255)
                                                                                                                                                                                                                                 cgacTACTGAGCGACTAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                     <  2:611681/248‑1 (MQ=255)
                                                                                                                                                                                                                                  gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTCTCGTTATCCATCACCGGGATTTCAATGTCGAGCATCCCGGCGAACTGGCGCATCTCGGTGAGGGTAAGCGCATGGCTGAAGACGGTATGGTGCGCGCCCCCAACGGGGGTCCACCCTTCGGAACGGAGT‑GGCAGGTCCGGGTGGGCTTTCCACCGCGCCATCGCCACCAGCaacttcg                    >  2:486716/3‑245 (MQ=255)
                                                                                                                                                                                                                                   acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCgg                   <  1:408356/250‑1 (MQ=255)
                                                                                                                                                                                                                                   acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGGGGTTTCCATGTCGTGCCTCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGCAGACGGTATGGTGTCCGCCACCGGCGGAGTGCCCCGCCTCGGGA‑GGAGTTGGCAGAGCCCGTTGCGCTTTCCACCGCCCATTCTCCCCCGGccgattccc                   >  2:274558/2‑243 (MQ=255)
                                                                                                                                                                                                                                      aCTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCAACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCa                 <  1:668344/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGGTTTCAATGTCGTGCCTCTCGGCGAACTGGCGAATATCGTTGAGGTTCAGAGCGTGGCTGGAGACGGTAGGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑ACAGTTGGCAGATCCGGTTGCGCATTCCACAGCGCCATCGCCACAGGCCGTTTCGGcgcg               >  1:295358/1‑248 (MQ=255)
                                                                                                                                                                                                                                        ttAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGATGGCAGGCGGGGATCGCTAACAATCACCGCGACTTCAATGACCAGCATCTCGGAGAACCGGCGAATCAAGTTGAGGTTCAGCACATGGATGAAGACGGTAAGGTGCGGGCCACCAGCGAGGATCCACGATTCGGAA‑GGCGTTGGCAGAACCGGTTCCGGTTTCCACAGCGCATTCGCCCACGGGAGTTTCGGCaagg              >  1:467742/1‑248 (MQ=255)
                                                                                                                                                                                                                                          aGCGACGAAATCCTGAATACACGTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGa             <  2:296435/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gaAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt      <  1:420463/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gaAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt      <  1:36424/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    cccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGttt  <  2:495778/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
GCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTT  >  NC_000913/66603‑67097

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