Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,158,383 G→A 10.8% W172* (TGG→TAG)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,3830GA10.8% 324.0 / 27.7 122W172* (TGG→TAG) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (54/54);  new base A (7/6);  total (61/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 123 total aligned reads displayed.

TATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTA  >  NC_000913/1158098‑1158670
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               
gatcgccgcTGTCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCGGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACAGGTTTAACCGGTTCTACCGTATTAAGCTGCCTGAGTATCGTCGCTTCTTTGCCGGTAACCGCTTTGTGACGATCATTTCTGGCCTGGCTGCCATCTTTACGGGCGTTGTGCTGTCCTTACATTTGgc                                                                                                                                                                                                                                                                                               <  1:715718/283‑1 (MQ=255)
    cgcgCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGTCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGccgccg                                                                                                                                                                                                                                                                                          <  2:131674/289‑1 (MQ=255)
         tagCCGCAGTTGTTGCCTATGGCATCATGGTTAAACCATGGCCGTGGTTGCGCCCCTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTTGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTgg                                                                                                                                                                                                                                                                                     <  1:807807/286‑1 (MQ=255)
              gCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTg                                                                                                                                                                                                                                                                                <  2:507849/289‑1 (MQ=255)
                    gttgCCTATGGCATCATGGTTAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCa                                                                                                                                                                                                                                                                         <  2:160090/289‑1 (MQ=255)
                         cTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGAcc                                                                                                                                                                                                                                                                     <  1:906129/289‑1 (MQ=255)
                          tATGGCATCATGGTTGAAACCATGGCCGGGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACGCCTGGCGGATACTAGCGTACTCGGGGGGATTATCTCCGGTGCGATCGCTGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTTCCGGTAAACGCGTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCt                                                                                                                                                                                                                                                                    <  2:975880/289‑1 (MQ=255)
                            tGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCtt                                                                                                                                                                                                                                                                   >  1:356256/1‑288 (MQ=255)
                                  catGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTCCTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCTGAGGGATTATCCCCGGTTCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAATCTGCCTGAGTATCTTGGCTTCTTTCCCGGTCAACGCTTTGTTCCGATCATTTCTGGCCTGGCTGCCATCTTTACCGGCGTTGTGCTGCCCTTCATTTGGCCTCCGATGGCTTCTGCATTCCAGACCTTCGCCCAg                                                                                                                                                                                                                                                            >  1:856431/1‑289 (MQ=255)
                                    tGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTg                                                                                                                                                                                                                                                          <  1:466800/289‑1 (MQ=255)
                                      ggttaaaCCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTggg                                                                                                                                                                                                                                                        <  2:483106/285‑1 (MQ=255)
                                       ttAAAATCATGGCCGTGGTTGCGGCACTGTTACTGCATTTATCTGCTGTAGAAAGCGCCTCTTAACACCTGGCGGATACTGGCGTACTCTGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTggg                                                                                                                                                                                                                                                        <  2:505923/288‑1 (MQ=255)
                                       ttAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACCTGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTggg                                                                                                                                                                                                                                                        <  2:716826/288‑1 (MQ=255)
                                       ttAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCt                                                                                                                                                                                                                                                                                 <  1:980251/263‑1 (MQ=255)
                                         aaaaCCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGACATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGgct                                                                                                                                                                                                                                                      >  2:572337/1‑288 (MQ=255)
                                                tGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCa                                                                                                                                                                                                                                              >  2:579390/1‑289 (MQ=255)
                                                 ggCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACAGCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTCAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCTATCTTTAGTGGCGTTGTGCTGTCCTTACATTTGGCCGCCGATTGGTTCTGCAATCGAGAACTTCGCTTAGTGGGCCGCTTACCa                                                                                                                                                                                                                                              >  2:825752/1‑289 (MQ=255)
                                                 ggCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGa                                                                                                                                                                                                                                                                       >  2:952474/1‑263 (MQ=255)
                                                 ggCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGa                                                                                                                                                                                                                                                                       <  1:952474/263‑1 (MQ=255)
                                                      tGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAAccc                                                                                                                                                                                                                                        >  1:304011/1‑289 (MQ=255)
                                                            cgcCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCa                                                                                                                                                                                                                                              <  2:721976/277‑1 (MQ=255)
                                                            cgcCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCa                                                                                                                                                                                                                                              >  1:721976/1‑277 (MQ=255)
                                                            cgcCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGt                                                                                                                                                                                                                                  >  1:288690/1‑289 (MQ=255)
                                                             gcCACGGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACCCCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGt                                                                                                                                                                                                                                  <  1:448134/288‑1 (MQ=255)
                                                                         gCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGCGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTACATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCAt                                                                                                                                                                                                                      >  1:538930/1‑289 (MQ=255)
                                                                          cATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTACATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCAtt                                                                                                                                                                                                                     >  2:686976/1‑289 (MQ=255)
                                                                          cATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTg                                                                                                                                                                                                                                                                                <  1:379650/229‑1 (MQ=255)
                                                                          cATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTg                                                                                                                                                                                                                                                                                >  2:379650/1‑229 (MQ=255)
                                                                            tttACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGGACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATTATTTCTGGCATGGCTGCGATCTTTACTGGCGTTGTGCTGTCCTTCTTTTGGCCGCCGATTGGTTCTGCAATCTAGACTTTCTCTCAGTGGGCTGCTTACCAGAACGCGGTAGTTGCGTTTGGCTTTTaa                                                                                                                                                                                                                  >  1:270230/1‑288 (MQ=255)
                                                                               aCCTGCTGCAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTCAACGCTTTGTGCCGACCATCTCTGGCCTGGCTGCCATCTTTACCGGCGTTGTGCTGTCCTTCATTTACCCGCCGATTGGTTCTGCACTCCAGCCCTTCCCTCCGTGGGCTCCTTCCCCGAACCCGGTACTTGCGTTTGGCATTTCCg                                                                                                                                                                                                                 >  1:495408/1‑287 (MQ=255)
                                                                               aCCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGGTCTGCAATCCAGACCTTCCCTCAGCGGGCTGCTTACCAGAACCCGGTAGTTGCGGTTGGCATTTAAGGt                                                                                                                                                                                                               >  2:190850/1‑289 (MQ=255)
                                                                               aCCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACCCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTCCCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTCCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGTCTGCTTACAAGAACCCGGTTGATGCGTTTGGCATTTACGGt                                                                                                                                                                                                               >  2:910630/1‑289 (MQ=255)
                                                                                   gctgAGGAAATCGCCCGTCACTACCTGGGGTATACGGGCGTACTCGACGGGATTATCTGCGGTGCGGTCGTAGCGTACATGTTTAACCGTGTCTACCGTATTAAACTGCCTGAGTATCTTGGCTTCTGTGCCGGTAAACGCTTTGTGCCGATCATTGCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCa                                                                                                                                                                                                           <  1:757698/289‑1 (MQ=255)
                                                                                    ctgAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGAGCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTTGGCTGCTTACCAGAACGCGGTAGTTGCGTTTGGCATTTACGGTTTCAt                                                                                                                                                                                                          >  1:566846/1‑289 (MQ=255)
                                                                                    ctgAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCGTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCATTCATTTGGCTGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAATGGGCTGCTTACCAGAACCCGGTAGTTCCGTTTGGCATTTACGGTTTCa                                                                                                                                                                                                           >  2:818106/1‑288 (MQ=255)
                                                                                      gaatacaTCGCCTCTAGACACCTGGCGGATACTGGCGTAGTCGGAGGGATTATCTCCGGTGCGTTCGCAGCGTACATGTTTAACCGTGTCTACCGTATTAAGCTGCCTGAGTATCCTGGCTTCTTTGCCGGTAAACGCTTGTTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTACATTGGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGCGCTGCTTACCAGAACCCGGTAGTTGCGTTCGGCATTTACGGTTTCATc                                                                                                                                                                                                         <  1:825752/283‑1 (MQ=255)
                                                                                          aaaTCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACg                                                                                                                                                                                                    <  1:579390/289‑1 (MQ=255)
                                                                                          aaaTCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACg                                                                                                                                                                                                    <  1:179390/289‑1 (MQ=255)
                                                                                               gCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGcc                                                                                                                                                                                               <  2:144223/289‑1 (MQ=255)
                                                                                                         acCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCTGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGtt                                                                                                                                                                                     <  2:280283/289‑1 (MQ=255)
                                                                                                         acCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGtt                                                                                                                                                                                     <  1:751436/289‑1 (MQ=255)
                                                                                                         acCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGtt                                                                                                                                                                                     <  2:356256/289‑1 (MQ=255)
                                                                                                           cTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACTGCTTg                                                                                                                                                                                   >  2:972660/1‑289 (MQ=255)
                                                                                                           cTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTg                                                                                                                                                                                   >  2:749958/1‑289 (MQ=255)
                                                                                                             ggCGGGTACTGTGGAACGCGGATGGATTATCTCCGGTGCGATTGAGACGTACATGTTTAACCGTTGCTACCGTATTAAGCTGCCTGAGTAGCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGt                                                                                                                                                                                 <  1:217437/289‑1 (MQ=255)
                                                                                                                   ttcTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTAGCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGcac                                                                                                                                                                           <  1:586882/287‑1 (MQ=255)
                                                                                                                   tACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGcac                                                                                                                                                                           <  2:658175/289‑1 (MQ=255)
                                                                                                                   tACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGcac                                                                                                                                                                           <  2:266396/289‑1 (MQ=255)
                                                                                                                        gCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTTGTACCGTTTGGTCTGCACcaca                                                                                                                                                                       >  2:869448/1‑288 (MQ=255)
                                                                                                                            attCGGAGGGATGTTCTCCGGTGCGATCGAGGCTTAATTGTTTAACCGTTTCTACCGTATTAAGCTGCCTTAGTATCTTGGCTTCTCTGCCGGTCAACGCTTTGTGCCGAGCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCAGCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTg                                                                                                                                                                   <  1:142542/286‑1 (MQ=255)
                                                                                                                             cTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGa                                                                                                                                                                 >  1:384085/1‑289 (MQ=255)
                                                                                                                             cTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGa                                                                                                                                                                 >  1:719689/1‑289 (MQ=255)
                                                                                                                             cTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGa                                                                                                                                                                 >  2:18422/1‑289 (MQ=255)
                                                                                                                                      aTTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTTGTACCGTTTGGTCTGCACCACATCTGGAACGTACCtt                                                                                                                                                        >  2:615735/1‑289 (MQ=255)
                                                                                                                                        tATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCCCATCTGGAACGTACCttt                                                                                                                                                       >  2:497739/1‑288 (MQ=255)
                                                                                                                                                      tagagtcttACAGGTTTGACCTTTACTACCGCATTGAGCTGCGGGTATATGTTGGCTTCTTTGCCGGGAATCGCTTATTGCCGCTCATGTCCGAGCTGTGTTCCATCTTTACAGCCGTTGTTCTGTCCCTCCTTTGGCTGCCGATTGATTCTGCAATCCAGACCTTCTCTCCGTGGGCGGCTTACCAGAACCCGAAAGTTGCGTTGAGCATTTACGGTTTCATCTAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTg                                                                                                                                         <  1:172140/280‑1 (MQ=255)
                                                                                                                                                      tCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTg                                                                                                                                         <  1:366858/288‑1 (MQ=255)
                                                                                                                                                      tCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTg                                                                                                                                         <  2:360954/288‑1 (MQ=255)
                                                                                                                                                          aGCGTACGTGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCACTCCAGACCTTCTCTCAGTGGGCTGCTTACAAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGACACACATCTGGAACGTACCTTTCCAGCTGCAGATTGGTGaa                                                                                                                                    >  2:388432/1‑289 (MQ=255)
                                                                                                                                                                      tAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgca                                                                                                                        >  1:536350/1‑289 (MQ=255)
                                                                                                                                                                      tAACCGTTTCTACCGTATTAAGCTGCCTGACTATCTTGGCTTCTTTGCCGGTAACCGCTTTGTGCCGACCATTTCTGGCCTGGCTGCCATCCTTACTGGCCTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCACACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTATTTGCGTTTGGCATTTACGGTTTCATCGTCCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGGACCTTTCCCGATGCCGATTGGTGAATACACCAACgcc                                                                                                                        >  1:252173/1‑288 (MQ=255)
                                                                                                                                                                        aCCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTTCCTGGTACCGTTTGGTCTGCACCCCATCTGGAACCTACCTTTCCAGATGCAGACTGGTGAATCCACCAACgcagc                                                                                                                      >  1:964812/1‑289 (MQ=255)
                                                                                                                                                                        aCCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgcagc                                                                                                                      >  2:685379/1‑289 (MQ=255)
                                                                                                                                                                             ttCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggtc                                                                                                                 >  1:953541/1‑289 (MQ=255)
                                                                                                                                                                             ttCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCATGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTCGTCTGCACAACATCTGGAACGTACCTTTCCAGATGCCGATTGGTGAATACACCAACGCAGcaggtt                                                                                                                 >  2:384038/1‑288 (MQ=255)
                                                                                                                                                                              tCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCCCAGTGGGCTGCTTACCAGAGCACGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTAACTTTCCAGATGCAGATTGGTGCATACCACAACGCAGtcggtca                                                                                                                >  1:545166/1‑289 (MQ=255)
                                                                                                                                                                              tCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTg                                                                                                                                                                                                                                                                                <  1:429761/129‑1 (MQ=255)
                                                                                                                                                                              tCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTg                                                                                                                                                                                                                                                                                >  2:429761/1‑129 (MQ=255)
                                                                                                                                                                                     tATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGTTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGAACCACATCTGGAACGTACCTTTCCAGTTGCAGTTTGGTGAATACACCAACGCAGCAGATCAGGtttt                                                                                                          >  2:339008/1‑288 (MQ=255)
                                                                                                                                                                                       ttAAGCTGCCTGAGTATCTTTGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCGGGCCTGGCTGCCATCTTTACTGGCGTTGGGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCa                                                                                                       <  1:304174/289‑1 (MQ=255)
                                                                                                                                                                                       ttAAGCTGCCTGAGTATCTTGGCTTCTTTGCCCGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCa                                                                                                       <  2:545166/289‑1 (MQ=255)
                                                                                                                                                                                       ttAAGCTGCCGGAGTATCTTGGCTTCTTTGCCGGTAATCGCTGTGTGCCGATCATTTCTGGACTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCGTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTGGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACAGCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCa                                                                                                       <  2:495408/289‑1 (MQ=255)
                                                                                                                                                                                               ccTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACa                                                                                               >  2:85722/1‑289 (MQ=255)
                                                                                                                                                                                               ccTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTAGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACCCCAACGCAGCAGGTCAGGTTTTCCACGGCGACa                                                                                               >  2:303985/1‑289 (MQ=255)
                                                                                                                                                                                               ccTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGAACATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGGCCGGTTTTCCACGTCGACa                                                                                               >  2:979874/1‑289 (MQ=255)
                                                                                                                                                                                                  gAGTATCTTGGTTTCTTTTCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTAGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACTCCAACGCAGCAGGTCAGGTTTTCCACGGCGACAtt                                                                                             <  1:979874/288‑1 (MQ=255)
                                                                                                                                                                                                    ggaTCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTcc                                                                                           <  1:893245/286‑1 (MQ=255)
                                                                                                                                                                                                    gTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTcc                                                                                           <  2:304011/288‑1 (MQ=255)
                                                                                                                                                                                                         ttGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCCAGGCCACATTCCGCGTta                                                                                     >  1:533570/1‑289 (MQ=255)
                                                                                                                                                                                                         ttGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTta                                                                                     >  1:336588/1‑289 (MQ=255)
                                                                                                                                                                                                                  ttGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCggg                                                                            <  1:85722/289‑1 (MQ=255)
                                                                                                                                                                                                                         tcagcgcgTTGTGCCAATCATTGCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTcc                                                                                           <  2:191177/263‑1 (MQ=255)
                                                                                                                                                                                                                           acgtcgttGTGCCGATCATTTTTTGCCTGGTTCCCATCTTTACTGGAGTTGTGATGTCCTTCATTTGGACGCCGATTGGTTCTGCCATCCAGACCTTCTCTCGGTGGGCTGCTTACCTGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTTGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTTAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGat                                                                   <  1:386834/283‑2 (MQ=255)
                                                                                                                                                                                                                                     gCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                         <  1:18422/289‑1 (MQ=255)
                                                                                                                                                                                                                                      ccGATCCTTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCGGTCCTTCAATTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCACCGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                         <  1:928861/288‑1 (MQ=255)
                                                                                                                                                                                                                                       cGATAATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTACATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                                         <  2:899828/288‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcctggcTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCt                                            >  2:519763/1‑288 (MQ=255)
                                                                                                                                                                                                                                                   gcctggcTGCCATCTTTACTGGCGTTGTGCTGTCCTTCACTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTCCCAAAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAATGTTGCCTGGTACCGTTTGGTCTGCACAACATCTGGAACGTACCTTTCCAGATTCAGATTGGTGAATACACAAATGCAGAAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTCCGGGGACACTGTCTGGTGGCtt                                           >  2:489446/1‑289 (MQ=255)
                                                                                                                                                                                                                                                     cgtgcTGCCATCTTTACGGCTGTTGTGCCGTTCTTCATTTGGCCGGGGTGTGGTTCTGCATGCCAGACCTTCTCTCAGTGGGCTCCTTACCCGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTTGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTcc                                         <  1:757536/286‑1 (MQ=255)
                                                                                                                                                                                                                                                        gcTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCGGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTCCGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGGAGATTGGTGTATACACCAACGCAGCCGGTCAGTTTTTCCACGGCGAACTTCCGCGTTATAGGGCGGGTGACCGGACTGCGGGGAAACGGTTgggtggatcacttg                                      >  1:583400/1‑281 (MQ=255)
                                                                                                                                                                                                                                                        gcTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGt                                      >  2:278424/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                  ttACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGcagcag                                                                                                                    <  2:725172/201‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGcagcag                                                                                                                    >  1:725172/1‑201 (MQ=255)
                                                                                                                                                                                                                                                                           ttGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGcc                   >  2:119086/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                           ttGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGc                    >  1:579843/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             gtgCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCa                  >  2:242965/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                               gCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGc                <  1:240307/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  gTCCTTACATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGc             <  2:538930/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          aTTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTg     <  1:151635/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ggCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTa  >  2:602190/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               
TATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAAGCTGCCTGAGTATCTTGGCTTCTTTGCCGGTAAACGCTTTGTGCCGATCATTTCTGGCCTGGCTGCCATCTTTACTGGCGTTGTGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTA  >  NC_000913/1158098‑1158670

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: