Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,620 G→A 6.0% intergenic (+47/‑45) yfjP → / → yfjQ CP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,768,6200GA 6.0% 320.0 / 4.9 99intergenic (+47/‑45)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base A (3/3);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

GTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACGGGTA  >  NC_000913/2768347‑2768880
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                     
gtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGGTTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAATTGCGGc                                                                                                                                                                                                                                                                      >  2:495741/1‑275 (MQ=255)
 tcgtgcgcgTGAGGCCAGCATCCCGGTGGTGGCGCTACTGTAAACTCGCGTTCGACCAACGGTAGCCCGGGAGCAGGAACGTTACGATTTCGGGTAAACCGGAGGGGCCATACTGGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTTGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGc                                                                                                                                                                                                                                                                      <  4:312815/269‑1 (MQ=255)
                        cgggattggcgcTACTGCAAACTACCTTGCGCACCAAGGTAGCCCGGGAGCAGGCACGGGACGATTTCGGTGAAACCGGAGGTGCAATACTCGTTACGGTAAGCACCTTTCCCCTCATGCCCGCCCCGGGGCGGACCATTATTCATGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                                                                                                                                               <  4:495736/266‑1 (MQ=255)
                        cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                                                                                                                                               <  2:237360/274‑1 (MQ=255)
                          gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                                                                                                                                                                                             <  2:56466/274‑1 (MQ=255)
                                  gcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGGAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                                                                                                                                                     <  2:212435/274‑1 (MQ=255)
                                     acggCAACAGCCCTGTCGAGCAACGGTAGCACGGGTGCAGGCAAGTGACGATTTCGGGCAAAACGTAGGGGCAATACTCGATTCGGTAAGCAACTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGAGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                                                                                                                                                                                                 <  4:241463/272‑1 (MQ=255)
                                                               tAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTTCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                                                                                                                                                                                                       >  2:357442/1‑243 (MQ=255)
                                                                      ggAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATTACCCGTCTGGCTTCGCGCTTTGgc                                                                                                                                                                                                >  2:348792/1‑275 (MQ=255)
                                                                                                aaaCCGTAGGGGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc                                                                                                                                                                       <  2:388869/274‑1 (MQ=255)
                                                                                                        gctGCCATACTCGTTACGGTAAGCACCTTTCCCCTTATCCCCGCCCCGTTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGTGTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGTGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTaa                                                                                                                                                              <  2:171471/273‑1 (MQ=255)
                                                                                                                    gATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTACTGAGTGTTTAATTCTGCACTAGTATTATCCGAGCACTGTCGCAAGTACGGCAGGGTTTTCTTTATTTGTTTTACCCATACTGAGGTGTCTGCTTATGACACGGCCGGCTTCGCGCTTTGGCTCTTCAACTCTTCTTCGCCGTGACCGTCTGTTacaccttggaggtt                                                                                                                                                  >  4:251962/1‑262 (MQ=255)
                                                                                                                              gCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGTTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt                                                                                                                                        <  2:109410/275‑1 (MQ=255)
                                                                                                                              gCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgt                                                                                                                                        <  2:251966/275‑1 (MQ=255)
                                                                                                                                  cTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAggg                                                                                                                                                                                                                                                                  <  8:461058/149‑1 (MQ=25)
                                                                                                                                     tCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctt                                                                                                                                                                                                                                                               <  5:273970/149‑1 (MQ=255)
                                                                                                                                        ccTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                                                                                                                                                                                            >  7:128744/1‑149 (MQ=255)
                                                                                                                                         cTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttctt                                                                                                                                                                                                                                                           >  8:459248/1‑149 (MQ=255)
                                                                                                                                          ttATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttcttt                                                                                                                                                                                                                                                          >  5:50129/1‑149 (MQ=255)
                                                                                                                                           tATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTa                                                                                                                                                                                                                                                         <  8:363270/149‑1 (MQ=255)
                                                                                                                                            aTTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                                                                                                                                                                                     >  4:435878/1‑168 (MQ=255)
                                                                                                                                                ccGCCCGGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                                                                                                                                                    >  8:176972/1‑149 (MQ=255)
                                                                                                                                                ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                                                                                                                                                    <  6:50129/149‑1 (MQ=255)
                                                                                                                                                ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                                                                                                                                                                                    >  7:268686/1‑149 (MQ=255)
                                                                                                                                                ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                                                                                                                                                                                                 <  2:435883/168‑1 (MQ=255)
                                                                                                                                                ccGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTgg                                                                                                                                                                                                 <  2:474262/200‑1 (MQ=255)
                                                                                                                                                  gCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGCTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCCTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACACGTCTGGCTTCGAGCTTTGGCGCTGCAAATCTTATTCGCCGTAACCGTCCGTAAACCAGTGAAGGGTTATCTCTCGTGGTGCCCAGCGTAtgtaggt                                                                                                                    >  4:16660/1‑269 (MQ=255)
                                                                                                                                                  gCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                                                                                                                                                                                                <  2:309783/167‑1 (MQ=255)
                                                                                                                                                     ccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtttt                                                                                                                                                                                                                                                 <  7:355226/147‑1 (MQ=255)
                                                                                                                                                     ccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtttt                                                                                                                                                                                                                                                 >  8:355226/1‑147 (MQ=255)
                                                                                                                                                     ccGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAc                                                                                                                                                                                                                                               <  7:459248/149‑1 (MQ=255)
                                                                                                                                                      cGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                                                                                                                                                                                              <  7:361469/149‑1 (MQ=255)
                                                                                                                                                                cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAg                                                                                                                                                                                                                                                                    >  7:302268/1‑117 (MQ=25)
                                                                                                                                                                cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAg                                                                                                                                                                                                                                                                    >  7:294066/1‑117 (MQ=25)
                                                                                                                                                                cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAg                                                                                                                                                                                                                                                                    <  8:294066/117‑1 (MQ=25)
                                                                                                                                                                cATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGACCTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGTCAg                                                                                                                                                                                                                                                                    <  8:302268/117‑1 (MQ=25)
                                                                                                                                                                  ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGttt                                                                                                                                                                                                                                                  >  5:350614/1‑133 (MQ=255)
                                                                                                                                                                  ttattCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGttt                                                                                                                                                                                                                                                  <  6:350614/133‑1 (MQ=255)
                                                                                                                                                                            cGTGCGTTCCTCTGTGGTGTCAGTGGCCCGCGCTGTCTGTGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATGTTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGaa                                                                                          <  1:18048/275‑1 (MQ=255)
                                                                                                                                                                            cGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATg                                                                                                                                                                                                                        <  8:345990/149‑1 (MQ=255)
                                                                                                                                                                             gTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAg                                                                                                                                                                                                                                                                    >  2:211331/1‑104 (MQ=255)
                                                                                                                                                                               gCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCgtgg                                                                                                                                      >  2:235946/1‑228 (MQ=255)
                                                                                                                                                                                 gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGct                                                                                                                                                                                                                                                                <  4:211328/104‑1 (MQ=255)
                                                                                                                                                                                 gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCg                                                                                                                                                                                                                   >  7:374698/1‑149 (MQ=255)
                                                                                                                                                                                 cttacTCGGTGGTTTAAGTGGCCCGCGCTGTCTTGGCTTTTTTCTTCGTGTTGTGTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCGAGTGCGGCATGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCGTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTa                                                                                     <  4:348787/271‑1 (MQ=25)
                                                                                                                                                                                  ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTac                                                                                    <  2:16660/275‑1 (MQ=255)
                                                                                                                                                                                   tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTc                                                                                                                                                                                                                 <  7:129890/149‑1 (MQ=255)
                                                                                                                                                                                   tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc                                                                                                                                  <  4:235942/228‑1 (MQ=255)
                                                                                                                                                                                    ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTTACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTccac                                                                                  >  2:382714/1‑275 (MQ=255)
                                                                                                                                                                                         gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCtt                                                                                                                                                                                                           >  8:413709/1‑149 (MQ=255)
                                                                                                                                                                                                 aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTc                                                                                                                                                                                                          <  7:242685/142‑1 (MQ=255)
                                                                                                                                                                                                 aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTc                                                                                                                                                                                                          >  8:242685/1‑142 (MQ=255)
                                                                                                                                                                                                 aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttt                                                                                                                                                                                                   >  8:45192/1‑149 (MQ=255)
                                                                                                                                                                                                   tGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTgg                                                                                                                                                                                                 >  8:60004/1‑149 (MQ=255)
                                                                                                                                                                                                       ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCt                                                                                                                                                                                             >  7:346328/1‑149 (MQ=255)
                                                                                                                                                                                                            cTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaaa                                                                                                                                                                                        <  7:60004/149‑1 (MQ=255)
                                                                                                                                                                                                             tGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAAt                                                                                                                                                                                       <  7:45192/149‑1 (MQ=255)
                                                                                                                                                                                                             tGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAAt                                                                                                                                                                                       >  5:152776/1‑149 (MQ=255)
                                                                                                                                                                                                              gTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATc                                                                                                                                                                                      <  7:413709/149‑1 (MQ=255)
                                                                                                                                                                                                               tCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                                                                                                                                                                                     <  6:100171/149‑1 (MQ=255)
                                                                                                                                                                                                                 tGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                                                                                                                                                   <  6:152776/149‑1 (MQ=255)
                                                                                                                                                                                                                 tGGGATCTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                                                                                                                                                                                   >  8:297586/1‑149 (MQ=255)
                                                                                                                                                                                                                        tcttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCACAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCctg                                                            >  4:221073/1‑261 (MQ=255)
                                                                                                                                                                                                                         cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGt                                                                                                                                                                           >  2:368340/1‑149 (MQ=255)
                                                                                                                                                                                                                         cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGt                                                                                                                                                                           >  4:522932/1‑149 (MQ=255)
                                                                                                                                                                                                                         cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCa                                              >  2:70715/1‑274 (MQ=255)
                                                                                                                                                                                                                         cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCGGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCCCCGTGACCGTCCGTTACCCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTAAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCCAAATTTCCCTTCTGGACAGCCTGCa                                              >  4:472756/1‑274 (MQ=255)
                                                                                                                                                                                                                            cttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTgg                                                        <  2:221077/261‑1 (MQ=255)
                                                                                                                                                                                                                             ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc                                                                                                                                                                       <  4:368335/149‑1 (MQ=255)
                                                                                                                                                                                                                             ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc                                                                                                                                                                       <  2:522937/149‑1 (MQ=255)
                                                                                                                                                                                                                               ctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgt                                                                                                                                                                     >  7:104092/1‑149 (MQ=255)
                                                                                                                                                                                                                                    tgtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTCTTCGCCGTGACCGTCCGtt                                                                                                                                                                <  7:281125/149‑1 (MQ=255)
                                                                                                                                                                                                                                    tgtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCg                                                                                                 >  2:214689/1‑212 (MQ=255)
                                                                                                                                                                                                                                      tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTa                                                                                                                                                                                                                                                <  7:152556/67‑1 (MQ=255)
                                                                                                                                                                                                                                      tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTa                                                                                                                                                                                                                                                >  8:152556/1‑67 (MQ=255)
                                                                                                                                                                                                                                      tttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCtt                                  >  2:444920/1‑273 (MQ=255)
                                                                                                                                                                                                                                       ttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATAc                                                                                                                                                                                                                                         >  8:128103/1‑73 (MQ=255)
                                                                                                                                                                                                                                       ttAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATAc                                                                                                                                                                                                                                         <  7:128103/73‑1 (MQ=255)
                                                                                                                                                                                                                                        tAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGt                                                                                             <  4:214686/212‑1 (MQ=255)
                                                                                                                                                                                                                                                cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAg                                                                                                                                                                                                                                                                    >  2:349039/1‑37 (MQ=255)
                                                                                                                                                                                                                                                cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                                                                                                                                                    <  7:176972/149‑1 (MQ=255)
                                                                                                                                                                                                                                                cccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAaga                                                                                                                                                    <  8:268686/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                                                                                                                                                   >  7:448519/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAagag                                                                                                                                                   >  5:424567/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   tattatTTTCCGTGCCCTGTGTCAAGTGGGGAAGGGCTGTCTTTATTTGTTTTTCCCTTACTTCGGAGTCTGCTTATGACCCGTCTGGCGTCGCCCTTGTGCCCTGCAAACCCTCTTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCTTGGTTGCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGc                   <  4:444915/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   tattatTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGtt                                                                                                                                                 >  7:72883/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTc                                                                                                                                            <  5:202969/148‑1 (MQ=255)
                                                                                                                                                                                                                                                         tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTCACCGTCCGTTAACCCGTGAAGAGTTATTTc                                                                                                                                            >  6:202969/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTcgcg                                                                                                                                         >  6:127388/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             cGAGCCCTGTCGCAAGTTCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGac         >  2:467338/1‑274 (MQ=255)
                                                                                                                                                                                                                                                             cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAAc                                                                                         >  4:288721/1‑195 (MQ=255)
                                                                                                                                                                                                                                                               aGCCCTGTCGCAAGTGCGGCAGGGCTCTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGATCTATTTCGCGTGGTACACAGCGTATTCAGTGAGGAAAAACACGCATCCCGCAGTGAACACTACACGTATATACACACAAATTACCGGCTGGGCAGCCTGCACCGGGAAGGCTTCCAGCCATTTCTTGCCTGTGCAGacc        >  4:148559/1‑273 (MQ=21)
                                                                                                                                                                                                                                                                 cccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTa                                                                                     <  2:288725/195‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGcc                                                                                                                                  >  6:168541/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                                                                                                                                                  >  6:4307/1‑85 (MQ=21)
                                                                                                                                                                                                                                                                  ccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTg                                                                                                                                                                                                  <  5:4307/85‑1 (MQ=21)
                                                                                                                                                                                                                                                                   cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGccc                                                                                                                                 >  7:250753/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGac         <  2:148560/268‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACACTTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACGGGTa  >  2:481029/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                                 gCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGa                                                                                                                   <  8:250753/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                     
GTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTCCCGCCCCGGTGCGGACCATTATTCAGGCCGTGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTCCGTTAACCCGTGAAGAGTTATTTCGCGTGGTGCCCAGCGTATTCAGTGAGGATAAACACGAATCCCGCAGTGAACGCTACACGTATATACCCACAATTTCCCTGCTGGACAGCCTGCAGCGGGAAGGCTTCCAGCCATTCTTTGCCTGTCAGACACGGGTA  >  NC_000913/2768347‑2768880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: