Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,984 2 bp→CT coding (1315‑1316/1656 nt) sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,9840TC93.3% 226.6 / ‑0.5 76D439S (GAC→AGC) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base C (36/34);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_00091375,9850CT92.1% 221.5 / 2.6 76D439S (GAC→AGC) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base T (36/34);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACAT  >  NC_000913/75749‑76227
                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                    
gtggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGAttt                                                                                                                                                                                                                                      <  1:17773/241‑1 (MQ=255)
 tggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCTTTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTCA‑GCCAGATATCGCTCTCGATCTCTCCTTAATGCCACTGATCGTAGCT‑GT‑TCTCTTTGAtttc                                                                                                                                                                                                                                     >  1:125343/10‑250 (MQ=255)
          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                                                                                                                                                        <  1:170534/244‑1 (MQ=255)
                                                                                    ggcGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAA‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATaaa                                                                                                                                                        >  2:170534/2‑245 (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:513627/249‑1 (MQ=255)
                                                                                                                                                                                             aa‑GCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCT‑GA‑TCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTCAGGATGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGt                                          <  1:560336/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)
                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                    
GTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTA‑GCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTC‑GA‑TCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTGATAAAAGGTTAGGGTGAGGCTTTCCAGGCCAGCCGGTTTTTCGCTCTTTATGGTGCGGGCATGGTGCCAACGGGGGAGCAGTCCATAAGCCGGGAACCACAGTTGCTGGTACTGTTCCTCAGCGAAATAGACCAGATTAGTTGGAGAAAGCACAT  >  NC_000913/75749‑76227

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: ‑