Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,253 C→T S225S (TCG→TCA)  thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,2530CT94.2% 273.6 / 3.1 87S225S (TCG→TCA) thiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (53/29);  total (57/30)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

GCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACC  >  NC_000913/72007‑72496
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
gCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCGGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAg                                                                                                                                                                                                                                                   <  1:297154/251‑1 (MQ=255)
       ccGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGGGTTAATTCGCCACCCGGTGATGCCGGTGTATATCGGTATTTTGCGTAAAGTGGTTGGGGGTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTATTATCAGGCCTACAAAACGCACTGTCCGGAGGGTGGCAAAGGTGTTCAGGCCGCCTCTGACCGGGCATACAAACACGTAATCCCCCATAGTGCTGaaacccccgc                                                                                                                                                                                                                                            >  2:637950/1‑243 (MQ=255)
           tGATTTCTGCCAGCGGCGTGGTGGCGCTGGTGGTGTTAATGCGCCGCCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCtt                                                                                                                                                                                                                                         <  2:223891/250‑1 (MQ=255)
            gATTGCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTGGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTa                                                                                                                                                                                                                                        <  1:708057/250‑1 (MQ=255)
              tttCTGCCTGCGGCTTGGTGGCGCTGTTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATTGTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTAcc                                                                                                                                                                                                                                      <  2:188197/250‑1 (MQ=255)
                   gCCATCGGCGTGGTTGCGCGGTTGGTTTTAAGTCGCAACCCGTTGATGCCGGTGTATATCGATATTTTGTGTAAAGGGGTGGGGGTTGAGGGGATTGTCTTG‑TGCGGGGTGCCTGATGCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAACCGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCt                                                                                                                                                                                                                                   <  2:38484/248‑1 (MQ=255)
                                  gtgttggtggGGTTAATTCGCCACCCGGTGATGCCGGGGTTTATCGATATTTTGGGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTctt                                                                                                                                                                                                                 <  2:183566/247‑3 (MQ=255)
                                      tggtggTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAccc                                                                                                                                                                                                             >  1:566792/1‑251 (MQ=255)
                                                aaTTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGAGGTTGGGGTTTAGGGGATTGTCTTG‑TGAGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCCTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCg                                                                                                                                                                                                    >  2:82102/1‑250 (MQ=255)
                                                        aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCACTAGTGCTGAAGCACTCTCCTTACCGCTCACCAACTCATTCGTCTTACCCTGCCAGGCGCTGCTcccc                                                                                                                                                                                           >  2:419050/1‑250 (MQ=255)
                                                           cGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCg                                                                                                                                                                                        <  1:112324/251‑1 (MQ=255)
                                                                   ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCgcccac                                                                                                                                                                                >  1:215203/1‑246 (MQ=255)
                                                                   ccGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGaccac                                                                                                                                                                                >  2:146419/1‑248 (MQ=255)
                                                                          atatCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGATCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGACCGTCGGCGACCACCACCg                                                                                                                                                                           >  1:420059/1‑249 (MQ=255)
                                                                               gATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGt                                                                                                                                                                    >  1:246578/1‑251 (MQ=255)
                                                                                 tatTTTGCTGAAAGTGTTTGGGGTTTAGGGTATTGTCTTG‑GCCGGGTTGCTGGATCCGCCGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                                                                                                                                   <  2:12912/250‑1 (MQ=255)
                                                                                  atTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTg                                                                                                                                                                   >  1:12912/1‑249 (MQ=255)
                                                                                   ttttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCCCGTGGCg                                                                                                                                                                >  2:649887/1‑251 (MQ=255)
                                                                                     ttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGAt                                                                                                                                                              >  1:364894/1‑251 (MQ=255)
                                                                                           aaaGGGGTTGGGGGTTAGGGGATTGTCTGGGTG‑GGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGGTAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGGCTTTCCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCgcgc                                                                                                                                                        <  2:441971/251‑1 (MQ=255)
                                                                                                   tGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCtt                                                                                                                                                >  1:659108/1‑251 (MQ=255)
                                                                                                    ggggTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCCTTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGATCGCGTGGCGCGCCGCGAGGCATCTTc                                                                                                                                               >  2:139350/1‑251 (MQ=255)
                                                                                                     gggTTTAGGGGATTGTCTTG‑TGGGGTTTGCTTGATGCGCCG‑TGGACCGCGTCTTATCAGGCCTTAAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                                                                                                                              <  2:543732/251‑1 (MQ=255)
                                                                                                     gggTTTAGGGGATTGTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTcc                                                                                                                                              >  1:374851/1‑251 (MQ=255)
                                                                                                         ttAGGTGATTGTCTTG‑TGCGGGTTGCCTGATGGGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAca                                                                                                                                           <  1:663953/250‑2 (MQ=255)
                                                                                                           aggggattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGGCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGTGGCGGTCCGCGCGGCATCTTCCACACtg                                                                                                                                        >  2:614622/9‑251 (MQ=255)
                                                                                                            ggggattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgt                                                                                                                                       >  1:207465/8‑251 (MQ=255)
                                                                                                               gattgTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCCCACtgtg                                                                                                                                      >  2:348239/5‑249 (MQ=255)
                                                                                                                   gTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGCGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                >  1:130820/1‑251 (MQ=255)
                                                                                                                   gTCTTG‑TGCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGacac                                                                                                                                >  1:348505/1‑251 (MQ=255)
                                                                                                                       tg‑tgCGGGTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAcc                                                                                                                               <  1:496612/248‑1 (MQ=255)
                                                                                                                              ggTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                                                                                                                                      >  2:117255/1‑155 (MQ=255)
                                                                                                                              ggTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                                                                                                                                      <  1:117255/155‑1 (MQ=255)
                                                                                                                               gTTGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGACGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATaacg                                                                                                                     >  2:422700/1‑248 (MQ=255)
                                                                                                                                 tGCCTGATGCGACGCTGA‑CCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTc                                                                                                                   >  1:402160/1‑251 (MQ=255)
                                                                                                                                       aTGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCAttttt                                                                                                             >  1:23102/1‑251 (MQ=255)
                                                                                                                                         gCGACGCTCA‑CCGCGCCTTATCCGGCCTACAAAACGCACTGCCCGTAGGTCGGATAAGTTTTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACATCCCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATTCGCGCGGCATCTTCCACACTGTGTGACACCAGCAATAACGTCATTTTTtg                                                                                                           <  1:207806/251‑1 (MQ=255)
                                                                                                                                              cgctgaccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTtgctgc                                                                                                       <  1:82102/245‑1 (MQ=255)
                                                                                                                                                     cgcgTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGc                                                                                                <  1:391518/251‑1 (MQ=255)
                                                                                                                                                            tATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑c                                                                                        <  1:649887/250‑1 (MQ=255)
                                                                                                                                                             aTCAGGCCTACAAAACGCACTACCCGTAGGTCGGCTAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCACCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCCTTTTTTGCTGCTGGCAGCTC‑GGG‑CTc                                                                                      >  2:325396/1‑251 (MQ=255)
                                                                                                                                                                aGGCCTACAACACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTACCCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGACGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGCGTGACACCATAAATAACGTCATTTTTTGCTGCTGGCACCTC‑GGG‑CTAAcc                                                                                   >  2:533882/1‑251 (MQ=255)
                                                                                                                                                                aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCAcc                                                                                   >  1:253871/1‑251 (MQ=255)
                                                                                                                                                                aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGGCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCAcc                                                                                   >  2:172565/1‑251 (MQ=255)
                                                                                                                                                                    cTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCg                                                                               >  1:552869/1‑251 (MQ=255)
                                                                                                                                                                      aCAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTc                                                                             <  1:419050/251‑1 (MQ=255)
                                                                                                                                                                       cAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCa                                                                            <  1:538629/251‑1 (MQ=255)
                                                                                                                                                                        caaaCGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGCCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTc                                                                             >  2:87633/2‑249 (MQ=255)
                                                                                                                                                                             gcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCAtct                                                                      >  1:508501/8‑251 (MQ=255)
                                                                                                                                                                              cactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCCTCTTCCACACTGTGTGACACCATCCATAACGTCCTTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCAtctc                                                                     >  2:44122/7‑251 (MQ=255)
                                                                                                                                                                                 tgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCCTCAATAACGCAATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCACCTg                                                                  >  2:149741/4‑251 (MQ=255)
                                                                                                                                                                                 tACCCGTAGGTCAGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAACATCTCCTg                                                                  >  1:55711/1‑251 (MQ=255)
                                                                                                                                                                                  gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCt                                                                   <  1:228813/247‑1 (MQ=255)
                                                                                                                                                                                  gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCt                                                                   <  1:144894/247‑1 (MQ=255)
                                                                                                                                                                                  gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTg                                                                  >  1:664746/3‑250 (MQ=255)
                                                                                                                                                                                    ccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACg                                                               >  2:522392/1‑251 (MQ=255)
                                                                                                                                                                                     cGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGt                                                              >  1:290148/1‑251 (MQ=255)
                                                                                                                                                                                      gTAGGTCGGATAAGGTGTTCACGCCGAATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCTGCGACCACCACCGAGCGCTTGGCGATCCGCGCGGCATCTTCCTCACTGTGTGCCACATTCCATAACGTCCTTTTTTGCTACTGGCAGCTC‑‑TGTCTCCCCAGCGTCAGCATCTCCTGACGt                                                              >  2:630387/1‑250 (MQ=255)
                                                                                                                                                                                       tAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTa                                                             >  1:215659/1‑250 (MQ=255)
                                                                                                                                                                                         ggTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACg                                                          >  1:588556/1‑251 (MQ=255)
                                                                                                                                                                                             ggATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCCGCATCTCCTGACGTAACGcagg                                                      >  2:485785/1‑248 (MQ=255)
                                                                                                                                                                                              gATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAca                                                                                                                                           >  1:122417/1‑166 (MQ=255)
                                                                                                                                                                                              gATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAca                                                                                                                                           <  2:122417/167‑2 (MQ=255)
                                                                                                                                                                                               atcaggtgTTCACGCCGCCTCTGACCGTGCATACTAACCCGTAATACCCAATAGCGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACACTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGAGTGGCGATCCGCGCGGCCTCTTCCACACTGTGGAACACCATCAATAACTTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCATCAGCATCAGCATCTCCTGACGCAACGCCGGac                                                    >  2:219768/8‑250 (MQ=255)
                                                                                                                                                                                               aTAAGGCTTTCACGCCGCACCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAACATCTCCTGACGTAACGCCGGAt                                                    <  2:659108/251‑1 (MQ=255)
                                                                                                                                                                                                  aggtgTTCACGCCGCATCTGACAGTGCATACTAACCAGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCCCCAGCGTCAGCATCTCCTGACGTACCGCCGGATCg                                                  >  2:389903/5‑250 (MQ=255)
                                                                                                                                                                                                    gtgTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTGACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAgc                                               <  1:515805/249‑1 (MQ=255)
                                                                                                                                                                                                            gccgcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAgaga                                        >  1:77195/1‑250 (MQ=255)
                                                                                                                                                                                                              cgcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACg                                     <  1:485785/251‑1 (MQ=255)
                                                                                                                                                                                                               gcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACgg                                    >  2:236088/1‑251 (MQ=255)
                                                                                                                                                                                                               gcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACg                                     <  1:236088/250‑1 (MQ=255)
                                                                                                                                                                                                                       cAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCACAGTCAGAATCCCCTGACGTAACGCCGGATCAAGCGAAGCGAACGTTTCATCg                             >  1:491846/1‑250 (MQ=255)
                                                                                                                                                                                                                        agagCATACTAACCCGTACTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGCTGCGCCCGTCGGCGACCACCACCGAGCGCGTGTCGATCCGCGCGACATCTTCCCCACAGTGTGACACCCTCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGAAGTAACTCCGGCTCGAGCGCCGAGAACGGTTCATCg                             >  2:122261/4‑249 (MQ=255)
                                                                                                                                                                                                                        aGTGCATACTAACCCGTACTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGAAGAAAACGGTTCATCgcg                           >  2:68298/1‑249 (MQ=255)
                                                                                                                                                                                                                        aGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATCACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCgcg                           >  1:392351/1‑249 (MQ=255)
                                                                                                                                                                                                                           gCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCa                         <  1:146419/250‑1 (MQ=255)
                                                                                                                                                                                                                             aTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATa                      <  1:207359/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaCCCGTAATCCCAAATAGAGCTGAAGCACTCGCCTTACCGCTCACCAACTCCTTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGCGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATACCGTCATTTTTTGCTGCTTGCAG‑TCAGTG‑CTCACCCGCGTCAGCATCTCCTGACGTCACGCCGCATCGCGCGCAgagacagcttcatctagcactcaaatc                 >  2:207164/1‑228 (MQ=255)
                                                                                                                                                                                                                                    cccGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTAATCGAGCAATacattcgg               >  2:707109/1‑244 (MQ=255)
                                                                                                                                                                                                                                          aTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTc         >  2:272069/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ccccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcg        <  1:172565/250‑1 (MQ=255)
                                                                                                                                                                                                                                             cccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcc      >  2:166643/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAAACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGAATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg      >  1:348263/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCAcc  >  1:446767/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCCCCACCGAGCGCGt                                                                                                                                                                    <  2:463212/90‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGt                                                                                                                                                                    >  1:463212/1‑90 (MQ=255)
                                                                                                                                                                                                                                                    gTGCTGAAGCACACGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCCGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGAATCTTCCACACTGTGTGACACAATCAAAAACGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CGCACCAGCATCAACATCTCCTGACGTAACGCCGGATCGATGGCAGCAACCGGTTCATCGAACAATAGAATCGGCtatgctcgcaca  >  2:274140/1‑237 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTG‑TGCGGGTTGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTG‑CTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACC  >  NC_000913/72007‑72496

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