Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,028 T→C G74G (GGT→GGC)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,0280TC98.8% 265.8 / NA 84G74G (GGT→GGC) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/53);  total (30/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TGCAAAAAACAGGAATAACCATGCAACCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGT‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGC  >  NC_000913/45787‑46270
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
tGCAAAAAGCAAGAATACCCCTGCAACCGT‑CCAGAAACTGTGACGCTCTCAATTTCTCCTCTATTCACCGCCGCATTTTGCTGTGGTGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGTTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCCAAAACTGGTCGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGTC‑TATATTTc                                                                                                                                                                                                                                            <  2:143646/250‑1 (MQ=255)
   aaaaaaCAGGAATCCCCATGCA‑CCGTGCCAGAAACTTTGACGATCTCAAATTCTCCTCTTTTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGTTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGAt                                                                                                                                                                                                                                        <  2:449561/251‑1 (MQ=255)
   aaaaaaCAGGAATACCCATGCAACCGT‑CCAGAACCTTGGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGAt                                                                                                                                                                                                                                        <  2:477478/251‑1 (MQ=255)
         caagAATAACCATGCAACCGT‑ACAAAAACTTTGACGATCTCAAATTCTCCTGTTTGCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTc                                                                                                                                                                                                                                  <  1:13268/248‑1 (MQ=255)
           ggAATAACCATGCAACCGT‑CCAGAACCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCgg                                                                                                                                                                                                                                <  2:484254/251‑1 (MQ=255)
                aaCCATGCAACCTT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgc                                                                                                                                                                                                                           <  1:540728/251‑1 (MQ=255)
                aaCCATGCAACCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgc                                                                                                                                                                                                                           >  2:564517/1‑251 (MQ=255)
                aaCCATGCAACCGT‑CCAGAAACGTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGCAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCGGGCGCTAAAACTGGACGCTGACTGGATTGGCTCGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGCTAAAGTCGGACggg                                                                                                                                                                                                                           >  2:535620/1‑250 (MQ=255)
                 cccATGCAACTGC‑ACAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGTAGTTTGCTGTGGGACAGCGGTTGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgcg                                                                                                                                                                                                                          <  1:108355/250‑1 (MQ=255)
                 cccATGCAACCGT‑CCAGAACCTTTGCCGACCGCAATGTCTCCTCGATTCCCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTTGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAATCTGGACGCTGTCTGGATTGGCTTGTTGGGCTCGGGAACGCTCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgcg                                                                                                                                                                                                                          <  2:201116/250‑1 (MQ=255)
                   caggCAACCGT‑CCAGAACCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCa                                                                                                                                                                                                                        <  1:564517/248‑1 (MQ=255)
                   caggCAACCGT‑CCAGAAACTTTGACGATCTAAAATTCCCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGGCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCa                                                                                                                                                                                                                        <  2:52506/248‑1 (MQ=255)
                      gCAACCGT‑CGAGAAACTTTGACGATCTCAAATTCCCCCCTATTCCCCGCCGCATTTTGCTGTGGGGAAGCGGTTGTCTGTTTCTGGATGGTTAGGTACTGGTAATGATGGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGTGCGCCGGGCTGTTCGTTGGCTCATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGTCGCaaaa                                                                                                                                                                                                                     <  1:404963/251‑1 (MQ=255)
                       cAACCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAAt                                                                                                                                                                                                                    >  2:506313/1‑251 (MQ=255)
                        aaCCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCCTTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCATAATg                                                                                                                                                                                                                   >  1:257115/1‑251 (MQ=255)
                        aaCCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATg                                                                                                                                                                                                                   <  2:574844/251‑1 (MQ=255)
                            cgtacagaaattttaactacctcAAATTCTACTCTATGCCCCGCCGCATTTTGCTGTGGTGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGTGCAACTGACGCCGGCGCTAAAACTGGATGCGGACTGGATTGGCTGGGTGGGCGCGGGAACTCTCGCCGGTCGGTTCGTGGGCACATCGCTGTGTGGC‑TATATTTCCGATAAAGTCGGGCGACGCAGAATGttt                                                                                                                                                                                                                <  2:38582/231‑2 (MQ=255)
                                cAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGCGGTCCGTTTCTGAATGGTTATGTACTGGTAATGATTGACGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGctcctc                                                                                                                                                                                                             >  1:127087/1‑250 (MQ=255)
                                 aGAAACTTTGCCGACCTTAAATGCTCCTCTATTCCCCGCCGCATTTTGCTGTGGGGAAGCGGTGTTCCGTTTCTGGATGGTTATGTACGGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGATGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACTGCGCAAAATGTTctt                                                                                                                                                                                                              <  2:514474/248‑3 (MQ=255)
                                   aaaCTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTCTGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATACAGTAGGTCGGCGCAAAATGTTCCTCAttt                                                                                                                                                                                                         >  2:30579/1‑250 (MQ=255)
                                         tGACGATCTCAAATTCTCCTTTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATatc                                                                                                                                                                                                    <  2:346876/250‑1 (MQ=255)
                                                      ttCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGt                                                                                                                                                                                      >  2:504521/1‑251 (MQ=255)
                                                          cctcTATTCCCCGCCGGATTTTGCTGTGTGGAAGCGGTGGTCGGTTTCTGTATGTTTATGTACTGGTAATGATTGGCGTGGCGCTGGTGCAACTGTCGCCGGCGCTAAAACTGGACGCTGACTGGATTTGCTTTTTGGGCGCGGGAACGCTCGCCGGGCTTTTTGTTGGCACAGCGCTGTTTTGC‑TATATTTCCGACAGTGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCTGCGTGa                                                                                                                                                                                    <  1:584726/249‑1 (MQ=255)
                                                             ctATTCACCGCCTCAGTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATGGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGAGTGGCTTGGTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGTCATCGGCGTGatat                                                                                                                                                                                 <  2:284198/249‑1 (MQ=255)
                                                             cgaTTCGCCGCCGCATTTTGCTGTGGGGAAGCGGTGGGCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGTACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTT‑GGTCTATATTGCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGatat                                                                                                                                                                                 <  1:68499/247‑1 (MQ=255)
                                                               aTTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGt                                                                                                                                                                             <  2:550686/251‑1 (MQ=255)
                                                                ttCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAAGGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTg                                                                                                                                                                            >  1:344934/1‑251 (MQ=255)
                                                                  ccccgccgCATTTTGCTGTGGGGACGCGGTGGTCCTTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTTATATCGGTgt                                                                                                                                                                           <  1:334663/248‑2 (MQ=255)
                                                                  cACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTgg                                                                                                                                                                           <  1:71622/250‑1 (MQ=255)
                                                                  cACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGc                                                                                                                                                                          <  1:74968/251‑1 (MQ=255)
                                                                   aCCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAATGTCGGACGGCGCAAAATGTTCCTCAGTGATATCATCGCCATCGGCGTGATATCGGTGGcg                                                                                                                                                                         >  2:42336/1‑251 (MQ=255)
                                                                           aTTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGAtgtt                                                                                                                                                                 <  1:506313/251‑1 (MQ=255)
                                                                                       ggAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATcccc                                                                                                                                                     >  1:289469/1‑251 (MQ=255)
                                                                                         aaGCGGGGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCCTTGATATCATCGCCATCGGCCTGGTATCGGTGGCGACGATGTTTGTTTCATCCCCCg                                                                                                                                                   >  2:549486/1‑251 (MQ=255)
                                                                                           gCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgc                                                                                                                                                                                                                           >  2:371407/1‑177 (MQ=255)
                                                                                           gCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGgc                                                                                                                                                                                                                           <  1:371407/177‑1 (MQ=255)
                                                                                           gCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTc                                                                                                                                                 <  1:228904/251‑1 (MQ=255)
                                                                                                   ccGTTTCTGGATGGTTATGTACTGTTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGGCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGt                                                                                                                                          <  2:342165/250‑1 (MQ=255)
                                                                                                    cGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTg                                                                                                                                        <  1:25166/251‑1 (MQ=255)
                                                                                                         cTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGAGCTAAAACTGGACGCTGACTGGATTGGCTTGGTGGGCGCGGGGACGCTCGACCGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGCTTTCTTCGCCATCGGCGTGATATCGGTGGCGAAGCAGGTTGTTTCATCCCCCGTCGTACTGTTGGTGAt                                                                                                                                   >  2:487455/1‑251 (MQ=255)
                                                                                                          tGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATg                                                                                                                                  >  1:359104/1‑251 (MQ=255)
                                                                                                           ggATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGc                                                                                                                                 <  1:539334/251‑1 (MQ=255)
                                                                                                                gTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTa                                                                                                                            <  1:148010/251‑1 (MQ=255)
                                                                                                                gTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTa                                                                                                                            <  1:67285/251‑1 (MQ=255)
                                                                                                                     gTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTAt                                                                                                                       >  2:550889/1‑251 (MQ=255)
                                                                                                                     gTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTAt                                                                                                                       >  1:479154/1‑251 (MQ=255)
                                                                                                                       aCTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATc                                                                                                                      <  1:576697/250‑1 (MQ=255)
                                                                                                                         tGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATc                                                                                                                                                                                >  1:276054/1‑190 (MQ=255)
                                                                                                                         tGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTGCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATc                                                                                                                                                                                <  2:276054/190‑1 (MQ=255)
                                                                                                                         tGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGAGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCCCATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATCTCGGTGGCGACGCTGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTTTCGGc                                                                                                                   >  2:549944/1‑251 (MQ=255)
                                                                                                                          ggTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGGACTTATCGGcc                                                                                                                  >  2:482580/1‑250 (MQ=255)
                                                                                                                              aTGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGt                                                                                                              <  1:319687/251‑1 (MQ=255)
                                                                                                                                   tggcgtggcgCTAGAGCAACTGGCGCCGGCGCTAAAACTGGACGCGGACTTGATTGGCCTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATc                                                                                                          <  1:27060/250‑1 (MQ=255)
                                                                                                                                    ggcgtggcgCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCgg                                                                                                        >  1:33931/1‑251 (MQ=255)
                                                                                                                                      cgtggcgCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGGCGGACGGCGCAAAATGTTCCTCATTGATATCATCTCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTGCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTg                                                                                                      <  2:258219/251‑1 (MQ=255)
                                                                                                                                         ggcgCTGGAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCa                                                                                                    <  2:119492/250‑1 (MQ=255)
                                                                                                                                               ggAGCAACTGACGCCGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTa                                                                                               <  1:131335/249‑1 (MQ=255)
                                                                                                                                               ggAGCAACTGAAGCCGGCGCTAAAACTGGACGCTGACTGGCTTGGATTGTTGGGCGCGGGAAAGCTCGCCGGGCTGTTCGTTGGCACATCGCTGATTGGC‑TATATTTCCGATAAAGTCGGCCGGCGCAACATGATCCTCATTGATATCATCGCCATCGGCGTGATATCGGAGGCGGCGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCGGATTATc                                                                                             >  2:300925/1‑251 (MQ=255)
                                                                                                                                                        gACGCCGGCGCTAACACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCCCATCGCTGTTTGGC‑TATATCTCCGATAAAGGCGGACGGCGCCAAATGTTCCTCATTTATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGct                                                                                     <  2:127087/250‑2 (MQ=255)
                                                                                                                                                            ccGGCGCTAAAACTGGACGCTGCCTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTACGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTc                                                                                >  1:96151/1‑251 (MQ=255)
                                                                                                                                                            ccGGCGCTAAAACTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTc                                                                                <  2:263780/251‑1 (MQ=255)
                                                                                                                                                                       aCTGGACGCTGACTGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCCTCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCCGTTTATCCCATCGCCACCTCAATGATCACCg                                                                     >  2:15839/1‑251 (MQ=255)
                                                                                                                                                                                    tGGATTGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTAc                                                        >  1:130341/1‑251 (MQ=255)
                                                                                                                                                                                       aTTGGCTTGGTTGGCGCGGTAACGCTCGCCGGGCTGTGCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAAGGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGGGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGGTGGTGATGCGGGGACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCATTGATCACAGAGTTCTCAAGTTCCCg                                                     <  2:503587/251‑1 (MQ=255)
                                                                                                                                                                                         tGGCTTGTTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTc                                                   >  2:134906/1‑251 (MQ=255)
                                                                                                                                                                                              tgttGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTTTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTAGCGGCATTGTCAGCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCg                                               <  2:59434/247‑1 (MQ=255)
                                                                                                                                                                                                 tGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCGGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAAGGTTCCTCATTGATATCATCGCCATCGGCGTGATAGCGGTGGCGACGATGTTTGTTGCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCAGCTCAATGATCACCGAGTTCGCCAGTACCCGTCAGCGGGc                                            <  2:359104/250‑1 (MQ=255)
                                                                                                                                                                                                  gggCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGt                                          <  1:34652/251‑1 (MQ=255)
                                                                                                                                                                                                        gggAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCa                                    <  1:581768/251‑1 (MQ=255)
                                                                                                                                                                                                         ggAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCAt                                   >  2:106016/1‑251 (MQ=255)
                                                                                                                                                                                                          gAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATc                                  <  1:141060/251‑1 (MQ=255)
                                                                                                                                                                                                              gCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGGGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCt                              >  1:103916/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGCCGTGCTGTTCGTTGGCACATCGCTGGTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTGCCATCAGCttt                            <  2:344934/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGCCGGGCTGTTCGTTGGCACAGCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGGGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCtt                             <  2:96151/249‑1 (MQ=255)
                                                                                                                                                                                                                   ccGGGCTGTTCGTTGGCACATCGCTGTTTGGC‑TAGATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTGTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCTGTGCAGATTAGCCCAGCGCCAGCTCAGTGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCTTGTTCCATCAGCTTTAtc                         <  2:298230/251‑2 (MQ=255)
                                                                                                                                                                                                                          gTTCGTTGGCACATCGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATAGCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATcgccgc                    <  2:479154/249‑1 (MQ=255)
                                                                                                                                                                                                                           tgcGTTGGCTCATCGCTGTTTGGC‑GATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATGGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCAGCCCCCGTCGAACTGGTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGGCAGCGGGCGTTTTCCATCAGCGTTATcgccgcca                  <  2:554195/248‑3 (MQ=255)
                                                                                                                                                                                                                                        cGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTTCCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgc    <  1:30579/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGTTATGTCGgc    <  2:289469/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgc    <  1:550889/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgc    >  1:45149/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCTGTTTTGC‑TATATGTCCGATCAAGTCGGGCGGCGCACACTGTTCCTTATGGATATCCGCTCCATCGGCGTGATCTCGATGGCGACAATGTTTGTTTCATCCCCCGTCGAACTGTTGGTAATGCGTGGACTTATCGTCATTATCATCTGTGCAGATTATCCCATCGCCACATCAATACTCACCCAGTTCTCCAGTACCCATCAGCGGACTCTTTCCCTCAGCTTGATCGCCGCCATGTGGTATgacgacg   >  2:177218/1‑245 (MQ=255)
                                                                                                                                                                                                                                         gCTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgcg   <  1:504521/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cTGTTTGGC‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATCGCCGCCATGTGGTATGTCGgcgc  >  1:40527/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
TGCAAAAAACAGGAATAACCATGCAACCGT‑CCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGGTGGTCCGTTTCTGGATGGTTATGTACTGGTAATGATTGGCGTGGCGCTGGAGCAACTGACGCCGGCGCTGAAACTGGACGCTGACTGGATTGGCTTGCTGGGCGCGGGAACGCTCGCCGGGCTGTTCGTTGGCACATCGCTGTTTGGT‑TATATTTCCGATAAAGTCGGACGGCGCAAAATGTTCCTCATTGATATCATCGCCATCGGCGTGATATCGGTGGCGACGATGTTTGTTTCATCCCCCGTCGAACTGTTGGTGATGCGGGTACTTATCGGCATTGTCATCGGTGCAGATTATCCCATCGCCACCTCAATGATCACCGAGTTCTCCAGTACCCGTCAGCGGGCGTTTTCCATCAGCTTTATTGCCGCGATGTGGTATGTCGGCGC  >  NC_000913/45787‑46270

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