Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 7.1% intergenic (+34/+26) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3360AC 7.1% 176.4 / 5.1 65intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base C (2/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

GAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACA  >  NC_000913/1159039‑1159629
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     
gAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                                                                                                                                                                                                                    >  2:195060/1‑301 (MQ=255)
       cAAAAGCGACAGGTTCCAGCGAAATGGAACCGGCTTTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                                                                                                                                                                                                                             <  2:182558/301‑1 (MQ=255)
              gACAGGTATCAGTGAAATGGCACCGGCTCTGGTTGCTGCATTTGGGGGTAAAGAAAACATTACTAACCGCGACGAATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGAACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttagtctcc                                                                                                                                                                                                                                                                                      <  1:294412/301‑8 (MQ=255)
               gccggtgccatttcgctGGTACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                                                                                                                                                                                                                     >  2:351683/17‑301 (MQ=255)
                   gTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTACACCGTAAGACGATGGGGAGACTAAGGCAGCAAGAAGGCTACATTTTTTACAggggg                                                                                                                                                                                                                                                                                 >  1:452990/1‑298 (MQ=255)
                            aaaTGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                                                                                                                                                                                                                        <  2:620225/301‑1 (MQ=255)
                              aTGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGGGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGAAGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGCCAGCCAAATGGCTGCCTTTTTTACAGGTGTTATACAAAAtt                                                                                                                                                                                                                                                                      >  1:582611/1‑301 (MQ=255)
                                 gtgtcgggtgtGGTTGCTGGAGTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCTTCTGCGCGTCAGCGTTGCTTCTGTGTCTAAAGTGGATCAGGGCGGCCTGAAGAAACTGGGCGCAGCGGGCGTATTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                                                                                                                                                                                                                                   <  1:351220/291‑1 (MQ=255)
                                  cACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATa                                                                                                                                                                                                                                                                  >  2:511042/1‑301 (MQ=255)
                                   tccGGCTCTGGTTGCTGCATTTGTTGTTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATAc                                                                                                                                                                                                                                                                 <  2:159475/300‑1 (MQ=255)
                                   aCCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATaa                                                                                                                                                                                                                                                                 >  2:430064/1‑300 (MQ=255)
                                     cGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGt                                                                                                                                                                                                                                                               >  2:226315/1‑301 (MQ=255)
                                          ttGGTTGCTTCATTTGGTGGTAAAGAAAACATTACTAACCTCTACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                                                                                                                                                                                                                          <  1:510370/300‑1 (MQ=255)
                                           tGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                                                                                                                                                                                                                         >  2:299918/1‑301 (MQ=255)
                                                gctgcATTTGGTGGTAAAGAAAACATTTCTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                                                                                                                                                                                                                    <  1:194817/301‑1 (MQ=255)
                                                  tgcATTTGGTTGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTg                                                                                                                                                                                                                                                  <  2:514096/301‑1 (MQ=255)
                                                           tggtAAAGAAAACATTACTAACCTCTAGGCATGTATTACCCGTTTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                                                                                                                                                                                                                         <  1:574162/301‑1 (MQ=255)
                                                                       cATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAGAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTaa                                                                                                                                                                                                                             >  2:978/1‑300 (MQ=255)
                                                                           aCTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGAGCCGGTAATGCTGAAATAACGCGGTATACCGTAGAcaa                                                                                                                                                                                                                         >  2:485784/1‑299 (MQ=255)
                                                                            ttttcttcgtctcATTTATTTCCCGTCTGCGCGTCAGCGTTGTTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCTTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAGGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGAt                                                                                                                                                                                                                        <  2:149519/289‑1 (MQ=255)
                                                                                 ctcgtggtaTGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                                                                                                                                                                                                                   <  1:242126/293‑1 (MQ=255)
                                                                                   ggacgtgtGTATTACCCGTCTTGGCGTATGCGTTGCTGATGTGGCTAATGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGTGGGCGTAGTGGTTGCTGGTTCTTGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt                                                                                                                                                                                                                 <  1:164909/294‑1 (MQ=255)
                                                                                       gCATGTATTACCCTGCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATGTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa                                                                                                                                                                                                             <  1:429500/301‑1 (MQ=255)
                                                                                                 cccGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATc                                                                                                                                                                                                   >  2:42248/1‑301 (MQ=255)
                                                                                                              gaGCGTTGTGGAGGGGTCTAAAGTGGATCAGGCCGGATGGAAGAAATGGGGCGCAGCGGGCGTAGTGGTGGCTGGTTGTGGTGGTCAGGCGATTTTCGGTACTAAAGCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAGTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCg                                                                                                                                                                                      <  1:497376/300‑1 (MQ=255)
                                                                                                                             gtCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGTTATCGTAGGTTTTGTCGAACAGGTTATTTAct                                                                                                                                                                       >  1:455214/1‑300 (MQ=255)
                                                                                                                             gtCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTCGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACg                                                                                                                                                                       >  2:205798/1‑301 (MQ=255)
                                                                                                                                   tgTGGATCAGGCCTGCCTGAAGAATCTGGGCGCTGCGGGCGTATTGTTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAAGTTATTGACGTTcccc                                                                                                                                                                 <  2:453554/300‑1 (MQ=255)
                                                                                                                                   aGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTcccc                                                                                                                                                                 >  2:569192/1‑301 (MQ=255)
                                                                                                                                     ttcttCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCAACGTTGGTATCGTAGGTTTTGTCGAAAAGGTTATTGACGTTCCCCTg                                                                                                                                                               >  2:642015/5‑301 (MQ=255)
                                                                                                                                          cAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTCGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACg                                                                                                                                                          <  1:205555/301‑1 (MQ=255)
                                                                                                                                            ggccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACgaa                                                                                                                                                        >  2:362568/1‑300 (MQ=255)
                                                                                                                                             gccggccTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACgaga                                                                                                                                                       >  2:554640/1‑301 (MQ=255)
                                                                                                                                                               ggCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGt                                                                                                                                     >  2:594856/1‑301 (MQ=255)
                                                                                                                                                                           gTAGTGGTTGCTGGTTCTGGTGTTCATGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAaga                                                                                                                         <  1:42195/301‑1 (MQ=255)
                                                                                                                                                                                tgTTGCTGGTTCTGGTGTTCAGGCTATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAGGTTTTTCGTCACCCGCTACCAGGTGAcgaaaaac                                                                                                                    <  2:455778/300‑6 (MQ=255)
                                                                                                                                                                                           cTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGCCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAACTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGCCGTTCCCCTGTAACGAGCAGTTTTTCGTCAACTAGTACAGGTGGAAGCGATCCACCAGCGCGtt                                                                                                         >  1:3139/1‑300 (MQ=255)
                                                                                                                                                                                             ggTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACGGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCAGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACATGGTAGCGGATGAAGAGATCAAACAGCGCGTAGc                                                                                                       >  2:338374/1‑301 (MQ=255)
                                                                                                                                                                                                   cAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCtt                                                                                                 >  2:360068/1‑301 (MQ=255)
                                                                                                                                                                                                                gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGaa                                                                                    <  1:568478/301‑1 (MQ=255)
                                                                                                                                                                                                                 aacataaTCCTATAACCTGTATACCGAGTTTGATGAGTACATCCTTAACCACTAATCCGTAAGACGTTTGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAg                                                                                   <  1:641138/296‑1 (MQ=255)
                                                                                                                                                                                                                  agtgaaggcggtAATCTGAAAATCTGGATTGTTTAGTACATCCGTATCCATTTAGCGGTAAGGTGTTTGGGGTACTAAGGCCGCCAGATTGCTGCCTGTTGTACAGGTGTTATTCAGAATTGATACGGGCTGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTAGCGTAGGTTTTGTCGAACAGGTTAGTGACGTGCCCCTGTAACGAGAAGGTTTTCGTCACCTGGTAGCGGGTGAAGAGATCAACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAgg                                                                                  <  1:368469/290‑1 (MQ=255)
                                                                                                                                                                                                                    tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTGTCTCGCCGCGGAAggag                                                                                >  1:88280/1‑299 (MQ=255)
                                                                                                                                                                                                                     aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                                                               >  2:211291/1‑301 (MQ=255)
                                                                                                                                                                                                                              aaCCTGAAAACCGAGATGGATGAGTACATCCGTAACCATTAATCCGTAAGTCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCg                                                                      <  1:299514/301‑1 (MQ=255)
                                                                                                                                                                                                                                ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTc                                                                    >  2:330449/1‑301 (MQ=255)
                                                                                                                                                                                                                                    aaaaCCGATATGTATGAGTACATCCGGAACCACGATTCCTGTAGTCGGTGGGTAGGATAGGGCAGCCAGATGGATGCCTGTTTGACAGGTGTTATTCAGAATAGATACGTGTCGGTAATGCTGAAATTACGCGGTGTGCCGTATACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGATGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                                                                <  1:337911/301‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGACACGgtg                                                              >  2:552029/1‑301 (MQ=255)
                                                                                                                                                                                                                                                 atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGCTAGCGGGTGACGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGCCACGTCGACGGTATacac                                                   >  1:509219/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          cttccgtttgttAATGCGTAAGACGTTGTGGAGATTAAGGCAGTCAGATGGCTGCTTTTTTTACAGGTGTTATTCAGAATTGATACTTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTACCCTGTAACGAGAAGTTTTTCGTAACCTGGTAGCGGGAGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTg                                          <  1:410068/290‑1 (MQ=255)
                                                                                                                                                                                                                                                            cGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCAAtttgtcc                                        >  2:229357/1‑297 (MQ=255)
                                                                                                                                                                                                                                                             gTAACCACTAAGCCGTAAGACGTTGGGGAGACTAAGGCATCCATATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCGGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCa                                       <  1:594095/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAg                                      <  2:583335/301‑1 (MQ=255)
                                                                                                                                                                                                                                                              aacctattAATCTGTAGGTCGTTGGGGGGACTGAGGCAGCCAGATGTCTGCCTTTTTTACTGGTGTTATTCAGAATTGATACGTGCCGGTAAGGCTGAATTTACGCGGTGTGCCGTAGACTATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAg                                      <  1:553967/294‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTaa                                  >  2:8966/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCACTGTAACGAGAAGTTTTTCGTCACCTGGTAGAGGGTGAAGAGATCAACCAGCGCGTACCTACATTGCACGGCGCGGAAGGTGCCATACGGCGTCACGGTGCCGGTATATACGCGATTTTGCCAGTTaa                                  >  1:450441/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                   ttcctttcacgacgctcttccgttctaaaGGCAGCCATCTTGCTGCTTTTTTTACAGGTTTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAac                                 <  1:198924/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTgagatcggaaa                                >  2:410631/1‑291 (MQ=255)
                                                                                                                                                                                                                                                                             aaGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACAAGCGATTTTGCCAGTTAACACCAccgcac                       >  1:83558/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                aCGTTGGGGAGACTAATGGAGCCTGTTGGCTGTGTTTTTTACAGGTTTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGAcc                    <  1:551356/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       tgtaaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCCTCACGGTGTCGGTCAACACCCGATTTTGCCAGTTAACACCCACCCCAACCGCCaaac             >  1:374570/8‑299 (MQ=255)
                                                                                                                                                                                                                                                                                             aaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAacagatcggaagagcacacgtctgaact       >  2:199167/2‑276 (MQ=255)
                                                                                                                                                                                                                                                                                               agccAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACAAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGAACGTCAACTATAGCATg     >  2:243076/4‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  cAGCCAGATGGTTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACa  <  2:3139/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  cAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACa  >  2:168339/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     
GAAGATGCAAAAGCGACAGGTACCAGCGAAATGGCACCGGCTCTGGTTGCTGCATTTGGTGGTAAAGAAAACATTACTAACCTCGACGCATGTATTACCCGTCTGCGCGTCAGCGTTGCTGATGTGTCTAAAGTGGATCAGGCCGGCCTGAAGAAACTGGGCGCAGCGGGCGTAGTGGTTGCTGGTTCTGGTGTTCAGGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGATTTTGCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACA  >  NC_000913/1159039‑1159629

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