Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,907 C→T intergenic (+197/+46) fimH → / ← gntP type 1 fimbriae D‑mannose specific adhesin/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9070CT98.9% 304.8 / NA 91intergenic (+197/+46)fimH/gntPtype 1 fimbriae D‑mannose specific adhesin/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (53/37);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCC  >  NC_000913/4549616‑4550201
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                      
gTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATCAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGAACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTATCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGtcat                                                                                                                                                                                                                                                                                                >  1:255784/1‑300 (MQ=255)
gTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGtcat                                                                                                                                                                                                                                                                                                >  1:299357/1‑300 (MQ=255)
gTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTNGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGATGGTCCGTTCCCCTCACCCTAACCCTCTCCCCGGCGGGGCGAGGGGACTGTCCGGGCATATTTTTAGACTTTGtcat                                                                                                                                                                                                                                                                                                >  2:246509/1‑300 (MQ=255)
gTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTAACTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGtcat                                                                                                                                                                                                                                                                                                >  1:478280/1‑300 (MQ=255)
gTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGAATTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCCATATTACCTGAAGCTAAAAAACTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTACAAGGAGGGGCGAGGGGACTGTCCGGGCACCTTTTTAGATTTTGtcat                                                                                                                                                                                                                                                                                                >  1:442073/1‑300 (MQ=255)
          aaCGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGc                                                                                                                                                                                                                                                                                      <  1:107796/300‑1 (MQ=255)
          aaCGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGc                                                                                                                                                                                                                                                                                      <  2:748819/300‑1 (MQ=255)
           aCGGCAAATTATGCACGTACCGGAGGGCAGTTGACTGCATGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGcc                                                                                                                                                                                                                                                                                     <  2:548024/300‑1 (MQ=255)
                  aTTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCAtt                                                                                                                                                                                                                                                                              >  1:595854/1‑300 (MQ=255)
                  aTTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTAACTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAAACCTCTCCCAGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCAtt                                                                                                                                                                                                                                                                              >  2:270159/1‑300 (MQ=255)
                           gTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGct                                                                                                                                                                                                                                                                      <  1:553528/299‑1 (MQ=255)
                              ccGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCCATATTAACTGAAGCTAAAAAACTGCATGTTAGCCCTTTGTCGGCCAGATAAGAAGAGCCAGCGTCGCATCTGGCATAAACAAAGCGCACTTCGCTGGTCTGTTCCCCTCACCCTCACCCTCTCCCCGGAGGAGCGAGGGGACTGTCCGGGCACCTTTTTAGATTTTGTAATCAGTTTGAGCCTGCCCTTGGTAGGttct                                                                                                                                                                                                                                                                    >  2:72462/1‑298 (MQ=255)
                              ccGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctct                                                                                                                                                                                                                                                                    >  1:420072/1‑298 (MQ=255)
                                ggAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                                                                                                                                                                                                 >  1:392282/1‑299 (MQ=255)
                                           gACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTTAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGc                                                                                                                                                                                                                                                     <  2:787923/300‑1 (MQ=255)
                                            aCTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCt                                                                                                                                                                                                                                                    >  1:491996/1‑300 (MQ=255)
                                            aCTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCAATTGGCAGGCTCTGGTGTCCTTTTACGCt                                                                                                                                                                                                                                                    >  1:627647/1‑300 (MQ=255)
                                             cTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTa                                                                                                                                                                                                                                                   >  1:547622/1‑300 (MQ=255)
                                                                     tgggCGTGACTTTTGTTTATCAAGAGAGAAATCACAGGACAGTGCGAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAATACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGTGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataa                                                                                                                                                                                                                           <  2:359139/298‑1 (MQ=255)
                                                                            gACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGcc                                                                                                                                                                                                                     >  2:465657/1‑299 (MQ=255)
                                                                                  tGTTTATCAATAAAGAAATCACAGGACATGGCTAATGCTGGTACGCAATAGTACCTGAAGCGAAGAACCTGCACGGTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCACATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGGTCCTCTCACCCTAACCTTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGct                                                                                                                                                                                                                                                                      <  1:553378/244‑1 (MQ=255)
                                                                                      tATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCNCCCCGGAGGGGCGAGGGGACTGTCCGTGCACANTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAg                                                                                                                                                                                                          <  2:603111/300‑1 (MQ=255)
                                                                                             aaaGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGANGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAcc                                                                                                                                                                                                    <  2:615743/299‑1 (MQ=255)
                                                                                             aaaGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAAACACGGAGTTGACCa                                                                                                                                                                                                   >  1:479457/1‑300 (MQ=255)
                                                                                                    tCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTcaa                                                                                                                                                                                             >  2:38671/1‑297 (MQ=255)
                                                                                                        aGGACATTGCGAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTCTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagcag                                                                                                                                                                                         <  2:435065/299‑1 (MQ=255)
                                                                                                               tGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGTCTGCAATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAAACACGGAGTTGACCAGCTCCAGCAGACCCCa                                                                                                                                                                                   >  2:39071/1‑298 (MQ=255)
                                                                                                               tGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCACCTGGCATAAAAAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCAGGAGGGACGAGGGGACTGTCCGGACACATTTTTAGATTTTGTCCTCAGTCTGAGCCTGCAATTGGAAGGCCCTGGTGCCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCa                                                                                                                                                                                   >  2:38418/1‑298 (MQ=255)
                                                                                                                        tGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTACTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACg                                                                                                                                                                        <  1:63849/300‑1 (MQ=255)
                                                                                                                                 atatTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTa                                                                                                                                                               >  1:57501/1‑300 (MQ=255)
                                                                                                                                      aCCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGAc                                                                                                                                                          <  2:666581/300‑1 (MQ=255)
                                                                                                                                         tGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCTTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTNTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAgg                                                                                                                                                       <  2:604899/300‑1 (MQ=255)
                                                                                                                                                        gCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCttt                                                                                                                                         >  1:596303/1‑299 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGAGGCGTCAACGTCGCATCTGGCATACACAAAGCGCACTTTGATGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGAAGGGGCGAGGGGACTGTCCGGGTACATTTTTAGATTTTGTCGTCAGTCTGAGCCTGCAATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTACTCATCAGCACAACCATCAGCCCACCCAAGGAGTTGAcc                                                                                                                                                                                                    >  2:710629/1‑237 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGNGGGGCGAGGGGACTGTCCGGGCACATTTTTAGNTTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaac                                                                                                                                     >  2:604468/1‑300 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGCGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCCACGTGTCTTTTCCTGCCAGGTCAAAGTAACCTTTGaac                                                                                                                                     >  1:54048/1‑300 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaac                                                                                                                                     >  1:176818/1‑300 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCAAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaac                                                                                                                                     >  1:177630/1‑300 (MQ=255)
                                                                                                                                                           cGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATCGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTCATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAACAGCTCCAGCAGACCCCAGGTTTTCAAAGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGaac                                                                                                                                     >  2:710635/1‑300 (MQ=255)
                                                                                                                                                                     ttGTAGGCCAGATAAGACGCTTCAGCTTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTATGGTCGAGTGGACTGTCCTGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGNTGTCCTTTTACGCTACCATGCTAATAATCAGCNCAATAATCAGCCCAACCACGGAGGTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAg                                                                                                                           <  2:599405/300‑1 (MQ=255)
                                                                                                                                                                     ttGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGTGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTTGCAGGCTCCGNTGTCCTTTTACGCTACCATGCTAATAATCAGCNCAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAg                                                                                                                           <  2:599993/300‑1 (MQ=255)
                                                                                                                                                                     ttGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGaaa                                                                                                                            >  1:249111/1‑299 (MQ=255)
                                                                                                                                                                     ttGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTCTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGCCAGGCCAAAATAACCTCTGACCAAACAgca                                                                                                                             >  2:165718/1‑296 (MQ=255)
                                                                                                                                                                     ttGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCAGGAGGGGCGAGGNGACTGTCCGGGCACATTTTTAGATTTTGTCATNAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAAAGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGaaa                                                                                                                            >  2:589344/1‑299 (MQ=255)
                                                                                                                                                                         agggcAGATAAGACGCGTCAGCGTCGCATCTGGCATAGGTAAAGCGCACTTGGCTGGTCTGTTCCCCTCACCCTAACCCTCTCTCCGGAGGGGCGATGGGATTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAACAATCAGCCCAACCATGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCGAGGTAACCTTTGaac                                                                                                                                     <  2:329814/282‑1 (MQ=255)
                                                                                                                                                                              aGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGt                                                                                                                  >  1:518657/1‑300 (MQ=255)
                                                                                                                                                                              aGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCAGGAGGGGCGAGGGGACTGTCGGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCACAACCAAGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAAACTTTGAACAACCAGAAAGATGCAtagt                                                                                                                  >  2:453088/1‑297 (MQ=255)
                                                                                                                                                                                   aGACGCGTCAGCTTCGTATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAAtg                                                                                                             <  2:424787/300‑1 (MQ=255)
                                                                                                                                                                                   aGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGTGGACTGTCCTTGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCNACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGAGGCATCGTTAAtg                                                                                                             <  2:565586/300‑1 (MQ=255)
                                                                                                                                                                                       gcgTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCTTAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTggg                                                                                                         <  2:299357/300‑1 (MQ=255)
                                                                                                                                                                                             gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAggg                                                                                                   <  1:600989/300‑1 (MQ=255)
                                                                                                                                                                                             gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAggg                                                                                                   <  1:627963/300‑1 (MQ=255)
                                                                                                                                                                                             gCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAggg                                                                                                   <  2:456650/300‑1 (MQ=255)
                                                                                                                                                                                               gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtt                                                                                                                                                                                                                                                           <  2:507875/146‑1 (MQ=255)
                                                                                                                                                                                               gTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtt                                                                                                                                                                                                                                                           >  1:507875/1‑146 (MQ=255)
                                                                                                                                                                                                    aTCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTgg                                                                                             <  1:417035/299‑1 (MQ=255)
                                                                                                                                                                                                       tGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGTGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCANAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAAcc                                                                                         <  2:571605/300‑1 (MQ=255)
                                                                                                                                                                                                            tAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCa                                                                                     >  1:490190/1‑299 (MQ=255)
                                                                                                                                                                                                             aaaCAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCa                                                                                     >  1:259424/1‑298 (MQ=255)
                                                                                                                                                                                                             aaaCAAAGCGCACTTTGCTGGTCCGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCCGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAAACAGGAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCagcagc                                                                                   >  1:259435/1‑300 (MQ=255)
                                                                                                                                                                                                             aaaCAAAGCACACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCAAGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGGTAATGTGGGTGAGGGTGTTGGAACCCGCAGc                                                                                   >  1:259315/1‑300 (MQ=255)
                                                                                                                                                                                                                    gtgtgttttgttggtttgtgcctttcatccTAACCCTCTCCTCTTAGGGTTGGGGGGACTGTCCGGGCGCATTTTTAGATTTTGTCGTCGGTCTGAGCCTGCCATTTGCAGGCTCTGGTGTCCTTTTGCGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCTAGCAGATCCCAGGTTTTCAGCGTGTCTTTTACTGACAGGTCAGAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAGTGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGcc                                                                            <  2:392750/272‑1 (MQ=255)
                                                                                                                                                                                                                           ttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCACCGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAAACAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACca                                                                     >  2:25916/1‑300 (MQ=255)
                                                                                                                                                                                                                           ttGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACca                                                                     >  1:390118/1‑300 (MQ=255)
                                                                                                                                                                                                                             gCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGTAGGGGCGAGTGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACTGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTGCTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACcagc                                                                   <  2:424528/300‑1 (MQ=255)
                                                                                                                                                                                                                             gCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACcagc                                                                   <  1:442800/300‑1 (MQ=255)
                                                                                                                                                                                                                               tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACcagca                                                                  <  2:645664/299‑1 (MQ=255)
                                                                                                                                                                                                                               tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAACCGTCGCCAGTACcagcag                                                                 <  2:627647/300‑1 (MQ=255)
                                                                                                                                                                                                                                ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTANCCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCagcagc                                                                >  1:561498/1‑300 (MQ=255)
                                                                                                                                                                                                                                ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCagcagc                                                                >  1:313798/1‑300 (MQ=255)
                                                                                                                                                                                                                                ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCAGTCGACAGTACCagcagc                                                                >  2:780589/1‑300 (MQ=255)
                                                                                                                                                                                                                                ggTCTGTTCCCCTCACCCTAACCCTCTCCACGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACcagca                                                                  >  2:552055/1‑298 (MQ=255)
                                                                                                                                                                                                                                 gTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCagcagc                                                                >  1:646355/1‑299 (MQ=255)
                                                                                                                                                                                                                                     gTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTCGGTGAGGGTGTTGGAAACCGCAGCCGTCGCCAGTACCAGCAGCGcc                                                             >  1:787659/1‑298 (MQ=255)
                                                                                                                                                                                                                                        ccccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGAtt                                                        >  1:79041/1‑300 (MQ=255)
                                                                                                                                                                                                                                              aCCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGc                                                   <  1:714193/299‑1 (MQ=255)
                                                                                                                                                                                                                                               cacTAACCCTCTACCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACNCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCa                                                 >  2:253886/3‑300 (MQ=255)
                                                                                                                                                                                                                                                   aaCCCTCTCCCCGTAGGGGCTAGGGGACTGGCCGGGCACATTTGTAGATTGTGTCATCAGGCTGAGCCTGCCATTTGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGTTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCa                                             <  2:442073/300‑1 (MQ=255)
                                                                                                                                                                                                                                                         ctcCCCGGAGGGGCGAGGGGGCTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAACAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTACAGCAGACCCCGGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAACGATGCATCGTTCATGTGGGTGGGGGTGTTGGAGCCTGCAGCCGCCGCCAGTACCAGCAGCGCCGGATTAACTCCAAcccactgac                                       >  2:55532/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                    ggCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGa                            >  2:445492/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                    ggCGAGGGGACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTAAGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAAAAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGa                            >  2:412756/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                           ggACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGa                     >  1:684439/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                            gACTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGAt                    <  1:246509/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              cTGTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCAGTCGACAGTACCAGCAGCGCAGGATTCACGCCAAACAGCTGACCAGTTACTGGATCTAgggtt                   >  2:425310/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                                gTCCGGGCGCATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCACCGTGTCTTTTACTGACAGGTCAACGTAACCTTTGAACAACCAGAGAGATGCATCGTTAATGTGGGTGAGGGTGTTNGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACTAGCTGACCAGTTGCTGGATCTAGGATTgca                >  1:172176/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTNGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTgca                >  1:172363/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                gTCCGGGCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTgc                 >  2:468356/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                    gggCACATTTTTAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCAc            <  2:550955/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            ttttAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATc    <  2:420072/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            ttttAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATc    <  1:72462/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ttAGATTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGTTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATccc  <  2:478280/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                      
GTCTGGGATTAACGGCAAATTATGCACGTACCGGAGGGCAGGTGACTGCAGGGAATGTGCAATCGATTATTGGCGTGACTTTTGTTTATCAATAAAGAAATCACAGGACATTGCTAATGCTGGTACGCAATATTACCTGAAGCTAAAAACCTGCACGTTAGCCCTTTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCTGGCATAAACAAAGCGCACTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCC  >  NC_000913/4549616‑4550201

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