Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,182,8090GA75.6% 135.8 / 47.8 79H7H (CAC→CAT) potDspermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (12/7);  new base A (35/24);  total (48/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTATCAACGTAATAGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAACA  >  NC_000913/1182563‑1183052
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
tttATCAACGTAATAGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATggg                                                                                                                                                                                                                                                  >  1:317163/1‑250 (MQ=255)
   aTCAACGTAATCGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGACATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCCAGCAGATGGCggg                                                                                                                                                                                                                                               >  2:411935/1‑249 (MQ=255)
   aTCAACGTAATAGGTTGAAGGACCCCCCAGATCATAGGCCCCTTTTTTGTATTTTTTCAGCTGGGCGTACATGGTTTCGTTCGACTCTTAATTCAAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCATGTACTCGGTCCAGTTG‑TAGAAATACAGCGTTTTTTTGTCATCGGCGTGAGCGG‑CGCTAATGCCCAGTGCCAAAGCACCCGCCGCGAGCAGGTGGCGTg                                                                                                                                                                                                                                               <  1:818982/251‑1 (MQ=255)
       aCGTAATAGGTTGAAGGACCCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTAATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGATGGCGTg                                                                                                                                                                                                                                               <  1:662465/247‑1 (MQ=255)
       aCGTAATAGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCCTGCCCAGTGCCAGAGCACCCGCCGCGAGCAGCTGGCGTGaggg                                                                                                                                                                                                                                           >  2:662465/1‑248 (MQ=255)
           aaTAGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCAtttt                                                                                                                                                                                                                                       >  1:200117/1‑251 (MQ=255)
                   gAAGGAAACACCAGATCATCGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCCGTGCCAGAGCACCCGCGGCGAGCAGGTGGCGTGACCATTTTTTctctt                                                                                                                                                                                                                                >  2:198875/1‑246 (MQ=255)
                     aGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTACGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGGGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAAc                                                                                                                                                                                                                             >  1:914597/1‑251 (MQ=255)
                       gAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTACTACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGTAATACAGCGTGTTGTTGTCATCGGCGTGCGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGAGGGCGTGACCATTTTTTCATTTAACGt                                                                                                                                                                                                                           >  2:741071/1‑251 (MQ=255)
                              caGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCtgt                                                                                                                                                                                                                    >  1:404153/1‑251 (MQ=255)
                                   cATAGGCTCCTTCTTTTTTTGTTTTCGTGTTGGCGTACTGGGTTTCTGTCGACTCATAAGTCAAATAGATAACCTTAATCCCGTTTTCTTGGGTGAAC‑TGTTAAAGCAGCCC‑GGGCGGCATGTCCTGGGTCCGGTTG‑TAAAAATCCAGCTTGTTGTTGTCATCGGGGTGAGCGG‑CGCTCATGCCCTGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTa                                                                                                                                                                                                                <  2:513059/250‑1 (MQ=255)
                                    aTAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACCTGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGATGGCGTGAACATTTTTTCATTTAAAGTCCCCTGTGTTAcc                                                                                                                                                                                                              >  2:173090/1‑251 (MQ=255)
                                         cACCGTCTTTGTATTTTTTCGGCTTGGCGTACATGTTTTCGTTCGACTCGTAAGTCGAATAAATACCCTTAATACCGTTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCCCGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTTTTGTCATCGGCGTGAGCGG‑CGCTCAGGCCCAGTCCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCg                                                                                                                                                                                                         <  1:584773/251‑1 (MQ=255)
                                              tCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCAACCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAAc                                                                                                                                                                                                                             >  2:805944/1‑226 (MQ=255)
                                                 gtgtagtttttCAGCTTGGCGTACATGGTTTCGTTCGATCCGTAAGACTAATAGTTATCCTTAAATCCGGTTTCTTTGGTGAAC‑TGTTCAAGCGGTCT‑GGGCGGCACGTACTCGGTCTAGTTG‑TAGAAATACATCGTTTTGTTGTCATGGGGTTGGGCGG‑CCCTCCTGCCCTTTGCCAGACCACCCGCCGCGAGCAGATTGCGTGACCATTTTTTCATTTAACGGCCCCTGTGTTACCTTTCGTTTTATc                                                                                                                                                                                                  <  2:899326/243‑1 (MQ=255)
                                                        ttttCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTACTACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGCGCAGATGGCGTGACCATTTTTTCCTTTAACGTCCCCTGTGTTTCCTTTCGTTTTATCACGTGCa                                                                                                                                                                                           >  2:925920/1‑250 (MQ=255)
                                                               cTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGAACGTACTCGGTCCAGTTG‑TAGAAATACAGCGCGTTGTTGTCATCGGCGTAAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGTGCAGATGGCGTGACCATTTTTTCATTTAACTTCCCCTGTTTTAACTTTCTTTTTTTCACGCGCAATAAg‑ccg                                                                                                                                                                                  >  2:543606/1‑249 (MQ=255)
                                                                   gcgtagttGGTTTCGTCCGACTGGAAAGTGGAAAAAATAACTTAAATCCCGTTTTTTTGGGTAAAC‑TGTACAAGC‑GTCCCGGGGGGCACGTACTGGGCCCAGTTG‑TAGAAATCCGGCGGGTTGTTGCCACCGGCGGGGGGGG‑CCCCCTTGCCCAGGGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAg‑ctggctg                                                                                                                                                                              <  2:142916/244‑1 (MQ=255)
                                                                        cATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACCCGGTCCAGTTG‑TAGCACTACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCCGTGCCAGAGCACCCGCCGCGGGCAGATGGCGGGACCATTTTTTCATTTCACGTCCCCTGTGTTACCTTTCGTTTTATCCCGAGCAATAAA‑Ccggctggcca                                                                                                                                                                          >  2:912792/1‑248 (MQ=255)
                                                                                cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGATGGCGGGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCCGGCAATAACCATc                                                                                                                                                                  >  2:23771/1‑250 (MQ=255)
                                                                                           gaaatCGAATAGATAACCTTAATCCCTGTTTTTTTGG‑GAACCTGTACAAGCAGCCC‑GGGCGGCACGTCCTCGGCCCAGTTG‑AAGAAATCCAGCGTTTTTTTTGCAGCGGGGGGAGCGG‑CGCTCATCCCCAGTGCAAGAGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGc                                                                                                                                                      <  1:912792/247‑1 (MQ=255)
                                                                                              gTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCCAGAGCACTAAG‑CCGGCTGGCAATAACCCTCACCAGCGACCGCAcc                                                                                                                                                   >  2:491584/1‑251 (MQ=255)
                                                                                                  aaTAGATCACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGGCCAGTTG‑GAGAAATACCGCATGTTGTTGTCATAGGCGTGCGCGGACGC‑CATGCCCAGTGACAAAGGAACCGCCGCGAGATGGTGGGGTGACCATTTTTTGAATTATAGTCCCCTGGTTTACCTTTCGCTTTATCACGAACAATAAGACTG‑CGGGCAATAACCAACAcc                                                                                                                                                               >  2:753977/1‑235 (MQ=255)
                                                                                                         acccTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCACCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGt                                                                                                                                        <  1:491584/249‑1 (MQ=255)
                                                                                                         aaCCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAAACATCACCAGCGACAGCACCAGTAATATGGt                                                                                                                                        >  1:96986/1‑251 (MQ=255)
                                                                                                              tAATACCGGTTTCTTGGGGGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACTTACTCGGTCCATTTG‑TAGAAATACAGCGTGTTGTTGCCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGcc                                                                                                                                    <  1:411935/250‑1 (MQ=255)
                                                                                                               aaTACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAg                                                                                                                                  <  1:23530/251‑1 (MQ=255)
                                                                                                                 tACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCCTTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAgcg                                                                                                                                >  1:233632/1‑251 (MQ=255)
                                                                                                                         tCTTTGGTGAAC‑TGTTAAAGCATGCC‑GGGCGGCGCGTACGCGGTCC‑GTTGGAAGAAATACGGGGTGTTGTTGGCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGAGTTAAc                                                                                                                         <  1:904732/250‑1 (MQ=255)
                                                                                                                          cTTTGGAGCAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCa                                                                                                                                                              >  1:244291/1‑212 (MQ=255)
                                                                                                                           tttGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCAGGTACTCGGTCCAGTTG‑TAGAATTACAGCGGGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGGGCAGATGGCGTGAc                                                                                                                                                                                                                                             >  2:290624/1‑133 (MQ=255)
                                                                                                                           tttGGTGAAC‑GGTTAAAGCAGTCC‑TGGCGGCACGTCCTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGCCATCGGCGGGAGCGG‑CGCTCATGCCCAGTGCAAGAGCACCCGCCGCGAGCAGATGGCGTGAc                                                                                                                                                                                                                                             <  1:290624/133‑1 (MQ=255)
                                                                                                                              ggTGAAC‑TGTTCACGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCCCCGCGGGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTCTTACCTTTCGTTTTATCCCGAGCAATAAG‑CTGGCTGGCAATAACCATCCCCAGCGCCAGACCCAGTAATATGGTTGCCAGCGCGTTAACTTCCgg                                                                                                                   >  2:793116/1‑251 (MQ=255)
                                                                                                                               gTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCCCCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGc                                                                                                                  >  2:293551/1‑251 (MQ=255)
                                                                                                                                 gAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCCGTAATATGGTTGCCAGCGCGTTAACTTCCGGCg                                                                                                                 >  2:270338/1‑250 (MQ=255)
                                                                                                                                  aaC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTGGGTCCAGTGG‑TAGAAATACAGCGTGTTGTTGCCATCGGCGTGAGCGG‑CCCCCATCCCCAGTCCCAAAGCTCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGAt                                                                                                               <  2:233632/251‑1 (MQ=255)
                                                                                                                                  aaC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCCGTAATATGGTTGCCCGCGCGTTAACTTCCGGCGAt                                                                                                               >  2:52630/1‑251 (MQ=255)
                                                                                                                                       gTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGc                                                                                                           <  1:293551/251‑1 (MQ=255)
                                                                                                                                        ttCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAAACATCACCCGCGACAGCACCAGTAATATGGTTGCCCGCGCGTTAACTTCCGGCGATACGcc                                                                                                          >  1:308412/1‑251 (MQ=255)
                                                                                                                                          cAAGCAGTCC‑TGGGGGCACGTCCTGGGTCAAGTTG‑TAGAAATACAGCGTTTTGTTTTCATGGGCGTGGCCGG‑CGCTAAGGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGc                                                                                                                  <  1:535770/241‑1 (MQ=255)
                                                                                                                                          cAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTTATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAAATTGGTTGCCAGCGCGTTAACTTCCGGc                                                                                                                  >  2:535770/1‑241 (MQ=255)
                                                                                                                                            atcggTCC‑TGGCGGCACGTCGTCGGTCCATTTG‑AGGAAATACAGCGTGTTTTTTTCATCGGCGTGGGCGG‑CGCTCAGGCCAAGTGCCAGATACCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGTCAATAACCATCACCAGCGACAGCACCAGTAATTTGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACt                                                                                                      <  1:600185/247‑1 (MQ=255)
                                                                                                                                             gCAGTCC‑TGGCGGCACGTATTCGGTCCGGTTG‑AAGAAATACAGGTTGTTGTTTTCATGGGCGGGAGCGG‑CGAGCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACtt                                                                                                     <  1:173090/251‑1 (MQ=255)
                                                                                                                                             gCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAAt                                                                                                                                                                                         >  1:679603/1‑168 (MQ=255)
                                                                                                                                              cAGTCC‑GGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACttt                                                                                                    <  1:52630/251‑1 (MQ=255)
                                                                                                                                               aGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGGGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCCGCGACAGCCCCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACtttt                                                                                                   >  2:2688/1‑250 (MQ=255)
                                                                                                                                                  cccgatcggcACGTACTCGGTCCATTTG‑AAGAAATACAGCGTTTTGTTGTCATGGGCGGGAGCGG‑CGCTCATGCCCAGGGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACGTCCGGCGATACGCCGACTTTGAc                                                                                                 <  1:804730/245‑1 (MQ=255)
                                                                                                                                                   c‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGCGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGAccc                                                                                               >  2:764149/1‑250 (MQ=255)
                                                                                                                                                         ggcACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAAATTCCGGCGCTACGCCGACTTTGACCATCGa                                                                                           >  1:912086/1‑250 (MQ=255)
                                                                                                                                                          gcACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAAt                                                                                         >  1:854662/1‑251 (MQ=255)
                                                                                                                                                          gcACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCCGCGACAGCACCCGGAATATGGTTGCCAGCGCTTTAACTTCCGGCGCTACGCCGACTTTGAccgttcgg                                                                                          >  2:804730/1‑244 (MQ=255)
                                                                                                                                                              gTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCCTGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGGGACCCTTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAA‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGAAACGCCGACTTTGAcccgcgcaaaaa                                                                                      >  2:30533/1‑240 (MQ=255)
                                                                                                                                                                ccTCGGTCCAGTTG‑TAAAAATACAGCGTGTTGTTGCCATCGGCGGGAGCGG‑CGCTCATGCCCAGTGCAAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAAttt                                                                                   <  1:154310/250‑1 (MQ=255)
                                                                                                                                                                aCTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGGGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGAc                                                                                                 >  1:120300/1‑237 (MQ=255)
                                                                                                                                                                   cGGTTCAGTTG‑TCGAAATACGGCGTTTTTTTGTCTCCGGCGGGGGCGG‑CGCTCATGCCGGGGGTCATTGCCCCCGCCGCGAGCAGATGGCGTGACCATTGTTTCATTTAACGTCCCCTGTGTTACCTTTCTTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCTCCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTGAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAAtttt                                                                                  <  1:861819/249‑1 (MQ=255)
                                                                                                                                                                         aGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGAAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGc                                                                           >  1:767699/1‑250 (MQ=255)
                                                                                                                                                                          gTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGGGAGCGG‑CCCTCATGCCCAGGGCCAGAGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACGTTGACCATCGAATAAAGTTTTAACGGCAg                                                                         <  1:793116/251‑1 (MQ=255)
                                                                                                                                                                                  aaaTACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTc                                                                   >  1:140059/1‑250 (MQ=255)
                                                                                                                                                                                   aaTACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGa                                                                                                                >  1:69132/1‑204 (MQ=255)
                                                                                                                                                                                         gCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGACGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCCGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACttt                                                           >  2:644020/1‑250 (MQ=255)
                                                                                                                                                                                           gtgtTGTTGTCATCGGCGTGAGCGG‑CGCTCCTGCCCAGGGCCAGAGCACCCGCCGCGAGCCGGTGGCGTGACCATTTTTTCATTTAACGCCCCCTGTGTTACCTTTAGTTTTATCACGAGCAATAAC‑CTCGCTGGCAATAACCCTCACCAGCGACAACACCCGTAATATGGTTGCCAGCGCGTTAACTTCCGGGGATACGCCGACTTTTAACCTCATAAAAAATTTTAACGGCAGAATTTCaacactttgg                                                         >  2:868523/1‑242 (MQ=255)
                                                                                                                                                                                             gttgttgtCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGCGCAGATGGCGTGACCATTTTTTCATTTCACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAAcc                                                                                                                                                                     >  2:830233/1‑141 (MQ=255)
                                                                                                                                                                                             gttgttgtCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGACCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAAcc                                                                                                                                                                     <  1:830233/141‑1 (MQ=255)
                                                                                                                                                                                               tgttgtCAGCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCg                                                      <  1:247922/250‑1 (MQ=255)
                                                                                                                                                                                                 ttgtCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAAAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGAGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCCGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGACACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTg                                                   >  2:183285/1‑251 (MQ=255)
                                                                                                                                                                                                      aTCGGCGTGAGCGG‑CGCTCAAGCCCAGTGCCAGAGAACCCGCCCCGGGCAGGGGGGGTGAACATTTTTTTATTTAACGGCCCCTGTGTTACCTTTCGTTTTTTCACGAGCAATAAG‑CTGGCTGGCAAAAACCATCCCCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCCTCGAATAAATTTTTAAAGGCAGGAATTCATAAATTTGGCCGGTggcgg                                               >  2:681191/1‑246 (MQ=255)
                                                                                                                                                                                                        cggcggGAGCGG‑CCCTCATGCCCGGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGacgaac                                             <  1:764149/250‑1 (MQ=255)
                                                                                                                                                                                                        cGGCGTGAGCGG‑CGCTCTTGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGacgaac                                             <  1:543606/250‑1 (MQ=255)
                                                                                                                                                                                                           gcgtgggcggccctcatcccgcgggCCAGCGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTGCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGAcgaacgaa                                          <  2:854662/227‑1 (MQ=255)
                                                                                                                                                                                                            gTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGacgaacga                                           >  1:738773/1‑248 (MQ=255)
                                                                                                                                                                                                              gAGCGG‑CCCCCAGGCCCGGTCCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGaaa                                      <  1:644020/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggcccccatgcccgtggCCAGAGCCCCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTAc                                  <  2:96986/235‑1 (MQ=255)
                                                                                                                                                                                                                     cgcTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTAcc                                 >  2:817191/1‑250 (MQ=255)
                                                                                                                                                                                                                      gcTCATGCCCAGGGCCAGAGCACCCGCCGCGAGCAGATGGCGTGCCCATTTTTTCATTTAACGTCCCCTGTGTTACATTTCGTTTTATCACGAGCCCTAAG‑CTGGCTGGCAATACCCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATGCGCCGACTTTGCGCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACca                                >  2:616784/1‑250 (MQ=255)
                                                                                                                                                                                                                       cTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTACCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCCCCAGCGACCGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACcac                               >  1:944790/1‑250 (MQ=255)
                                                                                                                                                                                                                              ccggtcccgcaccccccGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTTACTAACGAAGAAACTACCACATCGTCCa                       <  2:912086/237‑1 (MQ=255)
                                                                                                                                                                                                                               cAGTGCCAGAGCACCCGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCAt                      >  1:655176/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cccGCCGCGAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTaa           >  1:141582/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   cgAGCAGATGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTTCCACATCGTCCATCGACAGGGTAAAGCTTAaca  <  2:141582/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
TTTATCAACGTAATAGGTTGAAGGAACCACCAGATCATAGGCACCGTCTTTGTATGTTTTCAGCTTGGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAAC‑TGTTCAAGCAGTCC‑TGGCGGCACGTACTCGGTCCAGTTG‑TAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGG‑CGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAG‑CTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAATATGGTTGCCAGCGCGTTAACTTCCGGCGATACGCCGACTTTGACCATCGAATAAATTTTTAACGGCAGAATTTCATAACTTGGTCCGGTGACGAACGAAGAAACTACCACATCGTCCATCGACAGGGTAAAGCTTAACA  >  NC_000913/1182563‑1183052

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