Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 76,000 T→C K434E (AAA→GAA)  sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091376,0000TC94.9% 243.0 / 2.9 79K434E (AAA→GAA) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (38/36);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTG  >  NC_000913/75759‑76244
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ggCACCATTGCTTGGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGCTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGAGCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGAc                                                                                                                                                                                                                                             <  2:121576/251‑1 (MQ=255)
 gCACCAACGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGGGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTGA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGAcc                                                                                                                                                                                                                                            <  2:249012/251‑1 (MQ=255)
  cACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCt                                                                                                                                                                                                                                           <  1:217257/251‑1 (MQ=255)
      attgCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATAGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTg                                                                                                                                                                                                                                       >  2:513627/4‑251 (MQ=255)
               cTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGAGGCACCTCGCACAAGTGCGCGAACAGCGCAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATAACGCTCTCGATCCCTCCTTCATGCCACTGATCGTAGCT‑AA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGcccg                                                                                                                                                                                                                              >  2:541533/1‑249 (MQ=255)
                       cAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATctgct                                                                                                                                                                                                                      <  1:519518/251‑1 (MQ=255)
                           tGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                                                                                  >  2:560336/1‑251 (MQ=255)
                           ggcgggCACCAGTTCGCCAGGTTCATTTCGCCAGTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                                                                                                                                                                                                  <  2:635438/247‑1 (MQ=255)
                                        tcgccagggtcattgcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCTCCTCGCACAAGCGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAAt                                                                                                                                                                                                      <  2:246722/235‑1 (MQ=255)
                                            cgggttcattccGCCATTGCGCCAGCGGGCGGCGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAAGCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCAccc                                                                                                                                                                                                 <  2:277345/240‑1 (MQ=255)
                                             aggttCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg                                                                                                                                                                                                >  1:350258/4‑251 (MQ=255)
                                                   atttcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTc                                                                                                                                                                                          >  1:703806/4‑251 (MQ=255)
                                                       cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTccatg                                                                                                                                                                                      >  2:143041/1‑247 (MQ=255)
                                                       cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtg                                                                                                                                                                                      >  1:58436/1‑251 (MQ=255)
                                                           aTTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGAt                                                                                                                                                                                  <  1:317737/251‑1 (MQ=255)
                                                             tGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATc                                                                                                                                                                                 <  1:455182/250‑1 (MQ=255)
                                                             ggcgcCAGCGGGCGGGGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAATTTGGCGCTGGTAA‑GCCAGATATCGCTCTCTATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCTCCCGATGCTCACTGTGATcc                                                                                                                                                                                <  2:13816/250‑1 (MQ=255)
                                                               cgcCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCTGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTg                                                                                                                                                                              <  1:210514/251‑1 (MQ=255)
                                                                      gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                                                                         >  1:36657/6‑249 (MQ=255)
                                                                      gggcggcGTCGGCTTGCCAGTCGATCGGAATACAGGGTGGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCACTCTCGATCTCTCCTTCAGGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTaaagg                                                                                                                                                                       >  1:181635/6‑249 (MQ=255)
                                                                       ggcggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGGGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAgg                                                                                                                                                                      <  2:667825/247‑1 (MQ=255)
                                                                          ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                                                                         >  2:170534/2‑245 (MQ=255)
                                                                          ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                                                                         <  1:170534/244‑1 (MQ=255)
                                                                          ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGtc                                                                                                                                                                    >  1:301539/2‑249 (MQ=255)
                                                                             ggcGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACGGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATTCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                                                                                <  2:620400/249‑2 (MQ=255)
                                                                             gTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAgga                                                                                                                                                                <  1:551284/251‑2 (MQ=255)
                                                                                      ccAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCttt                                                                                                                                                       >  1:702439/1‑251 (MQ=255)
                                                                                      ccAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCttt                                                                                                                                                       >  1:268914/1‑251 (MQ=255)
                                                                                       cagtcgatcggcatacagtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTACAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTACATCGTGACCTGGTGACTTGCCAGCATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCCGGATGAGGCTTTc                                                                                                                                                      >  2:518493/19‑251 (MQ=255)
                                                                                        agtcgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTcc                                                                                                                                                     <  1:12260/243‑1 (MQ=255)
                                                                                           cgatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAgg                                                                                                                                                  >  1:494151/6‑251 (MQ=255)
                                                                                            gatcggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGCGGCTTTcccggg                                                                                                                                                 >  2:541974/5‑247 (MQ=255)
                                                                                               cggAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCGTTCCAGGCCAg                                                                                                                                              <  2:242798/250‑1 (MQ=255)
                                                                                                        gtgttggAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGCTAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGttttt                                                                                                                                     >  1:397623/2‑251 (MQ=255)
                                                                                                                           tCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGATGTTTA‑GCCAGATATCATTTTCGATCGCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCgg                                                                                                                    <  1:249174/249‑1 (MQ=255)
                                                                                                                                    tgcgcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGcc                                                                                                          >  2:249124/19‑250 (MQ=255)
                                                                                                                                       gcgaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGATGCCAACg                                                                                                      >  1:610071/16‑251 (MQ=255)
                                                                                                                                         gaacagcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAC‑GCCAGATATCGCTCTCGATATCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGCTGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCACTTTATGGGGCGGGCCTGGTGCAAACggg                                                                                                    >  2:41020/14‑251 (MQ=255)
                                                                                                                                             agcgaaaaatCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAg                                                                                                <  1:9836/242‑1 (MQ=255)
                                                                                                                                                           cGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCg                                                                                  <  1:541119/251‑1 (MQ=255)
                                                                                                                                                             gCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                                                                                <  1:41020/251‑1 (MQ=255)
                                                                                                                                                             gCAGCGTAAAGTTGGCGC‑GTTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTTACCTGGTGACTTGCCAGAATCCGCTGCATGAGCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCggg                                                                                <  2:449298/251‑1 (MQ=255)
                                                                                                                                                                 cgTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCATTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGCCCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTACAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAAcc                                                                            >  2:609615/2‑251 (MQ=255)
                                                                                                                                                                   tAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACcac                                                                          >  1:121123/1‑251 (MQ=255)
                                                                                                                                                                   tAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGAATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGACCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGACCATAAGCCGGGAACcac                                                                          >  1:625164/1‑251 (MQ=255)
                                                                                                                                                                         tGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGc                                                                    >  1:26759/1‑251 (MQ=255)
                                                                                                                                                                           gcgcTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGATCCACAGTTGCTg                                                                  <  1:577362/251‑1 (MQ=255)
                                                                                                                                                                           gcgcTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTg                                                                  <  1:167223/251‑1 (MQ=255)
                                                                                                                                                                               tGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGACAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACTGTTGCTGGTAc                                                              >  1:701329/1‑251 (MQ=255)
                                                                                                                                                                                gttaa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACt                                                             >  2:595019/5‑251 (MQ=255)
                                                                                                                                                                                   aa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt                                                           <  1:560336/249‑1 (MQ=255)
                                                                                                                                                                                   aa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt                                                           <  1:513627/249‑1 (MQ=255)
                                                                                                                                                                                          gATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAg                                                    <  1:550940/251‑1 (MQ=255)
                                                                                                                                                                                          gATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCATTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTATTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCGGGTACTGTTCCta                                                      >  1:583831/1‑248 (MQ=255)
                                                                                                                                                                                           atatCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGc                                                   <  1:333883/251‑1 (MQ=255)
                                                                                                                                                                                             atcgctctCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCg                                                  >  1:407483/8‑250 (MQ=255)
                                                                                                                                                                                                gctctCGATCTCTCCTACATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGCGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCAGGGAACCCCAGTTGCAGGTACTGATCCTCAGAGaaa                                               >  2:507126/5‑250 (MQ=255)
                                                                                                                                                                                                     cGATCTCTCCTTCATGCCACTTATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGa                                           <  1:510839/249‑1 (MQ=255)
                                                                                                                                                                                                       attcctctttcaTGCCATTGATCGTAG‑C‑TACTCTCTTTGAATTCCAGGGTGACCTGGTGATTTGACCGAATCTGCTGCAAGATCCAGGCAATCACCCGATGTGCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGACGGTTTTTCGCTCTTTATGGTGCTGGCATGGTGCCAACGGGGGAGCAGTCAATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAg                                       <  2:526431/240‑1 (MQ=255)
                                                                                                                                                                                                       atccctccttcaTGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAg                                       <  1:536089/240‑1 (MQ=255)
                                                                                                                                                                                                              cttcaTGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCCTAAGCCGGGAACCCCAGTTGCTGGTACTGTTCCTCAGCGACATAGCCCAGATTAGtt                                >  2:263731/5‑251 (MQ=255)
                                                                                                                                                                                                                tcaggCCACTGATCGTA‑TCTGA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGtt                                <  2:301539/245‑1 (MQ=255)
                                                                                                                                                                                                                tcaTGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTg                               <  2:36657/248‑1 (MQ=255)
                                                                                                                                                                                                                gcaTGCCTCTGATCGTAGCT‑GA‑TCCCTTTGATTGCCAGCGTGACCTGGGGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCATCGGGGGAGCAGTCCATAAGCCGGGAATCAGAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAggtg                               <  2:19690/248‑4 (MQ=255)
                                                                                                                                                                                                                 caTGCCGCTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGGTGCTCACTGTGATCCTGATAAAAGGACAGGATGAGGCTTTACAAGCCAGCCGGTTTTTCGCTATTTATGGTGCGGGCCAGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACCGTTGATGGTACTGTTCCACAGCGGAATAGACAAGATTAGTTgg                              >  2:549031/2‑250 (MQ=255)
                                                                                                                                                                                                                  aTGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGga                             <  2:268914/250‑1 (MQ=255)
                                                                                                                                                                                                                  aTGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATACGACGGGAACCACAGTTGCTGGTCCTGTTCCTCAGCGAATTAGACCAGATTattatggg                            >  2:98335/1‑244 (MQ=255)
                                                                                                                                                                                                                    gCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGaa                          <  1:249124/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGaa                          >  1:504524/1‑250 (MQ=255)
                                                                                                                                                                                                                     ccACTGAGCGTAGCT‑GA‑TCTCGTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCGGCATTATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAgg                                                                                                                                                           <  1:101632/121‑1 (MQ=255)
                                                                                                                                                                                                                      cACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTa                                                               <  1:274950/212‑1 (MQ=255)
                                                                                                                                                                                                                      cACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACAACAGTTGCTGGTa                                                               >  2:274950/1‑212 (MQ=255)
                                                                                                                                                                                                                          gATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGcac                     >  1:359475/1‑250 (MQ=255)
                                                                                                                                                                                                                           aTCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACAt                   >  1:447936/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tgA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCAAAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGCTTAGTTGGAGAAAGAACATAGCTTAcc           >  1:487205/2‑251 (MQ=255)
                                                                                                                                                                                                                                         ctctTTGATTTCCAGCGTGACCTGGCGACTTGCCCGAATCTGCTGCATGATCCAGGCAATCACCCGATGCTAACTGTGATCCTGATAAAAGGTAAGGCTGAGGATTTCCAGGCCAGCCGGTTTTTCGATCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGAACAGATTAGTTGGAGAAAGCACCTAGCTTAccc          >  1:657599/1‑248 (MQ=255)
                                                                                                                                                                                                                                           ctTTGATTTCCAGCGTGACCTGGTGACTTGCCAGATTCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCCGATTAGTTGGAGAAAGCACATAGCTTACCCAGTcc     >  2:347638/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGTTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTg  <  2:359475/250‑1 (MQ=255)
                                                                                                                                                                                                                                               gATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTg  >  1:340032/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
GGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACATAGCTTACCCAGTCCCTG  >  NC_000913/75759‑76244

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