Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,214 C→T 5.4% intergenic (+20/‑210) mtlA → / → mtlD mannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,774,2140CT 5.4% 1905.8 / 63.7 569intergenic (+20/‑210)mtlA/mtlDmannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (259/273);  new base T (14/17);  total (276/293)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-06

Only 100 of 599 total aligned reads displayed.

CACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGG  >  NC_000913/3773940‑3774484
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                              
     tacggtTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGga                                                                                                                                                                                                                                                                           >  4:203171/1‑275 (MQ=255)
          ttGAAGCGTAAAATCGCGGACTGAAAACGTGCGTCTGGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGGGATTGGTATGGCTGCCCGGAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                                                                                                                                                       <  4:244100/274‑1 (MQ=255)
             aaGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGGGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   <  2:129926/275‑1 (MQ=255)
               ggaaaaGATTGCGGACGGACAACGTGNGNGGTGTTGTGCCAGTACCCGGAAGGTGTGCTCTTCGGTGAAGAAGANGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCNTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   <  3:48875/270‑1 (MQ=255)
               gCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTggg                                                                                                                                                                                                                                                                  >  1:14749/1‑274 (MQ=255)
                                         cgtcgtGTGCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTAGTGCGGCCCGTAACAACGGGCACATTCAGGTGATCACCAGCCTGACCAATGCAATGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACGACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTat                                                                                                                                                                                                                                       >  2:332943‑M2/1‑273 (MQ=255)
                                                        gTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACCACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGTGTGGATGAAGTCCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAAAAGGTGGGGTTAATTGCCTGATGCGCTCCGCTTatcagcactaccggaag                                                                                                                                                                                                                        <  4:438932‑M2/275‑18 (MQ=255)
                                                                               gTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTTCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGAATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAGTCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCa                                                                                                                                                                                                 <  2:498300/275‑1 (MQ=255)
                                                                                  aaaaagagtaTGCCATTGCCCCTCTGGTGATTGGTATTGTCGCCCGTAACAACGAGAACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGGGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACg                                                                                                                                                                                               <  4:504393/265‑1 (MQ=255)
                                                                                        aagaGGACAGTGGCCGTCTTGGGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCAATGGATGAGGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTg                                                                                                                                                                                         <  4:369165/274‑1 (MQ=255)
                                                                                          gaTGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGt                                                                                                                                                                                        >  1:101679/1‑273 (MQ=255)
                                                                                          gaTGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCAGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGt                                                                                                                                                                                        >  1:117185/1‑273 (MQ=255)
                                                                                                  tGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATACCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCGCAATTTGTTGAATTTTCTCGTTCTTCTGAGCCGGAt                                                                                                                                                                               >  1:62610/1‑274 (MQ=255)
                                                                                                              gATTGGTATTGCTGCCCATAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcg                                                                                                                                                                   >  2:159015/1‑274 (MQ=255)
                                                                                                                       tgctgcCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCgg                                                                                                                                                         <  2:357640/275‑1 (MQ=255)
                                                                                                                       tgctgcCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCgg                                                                                                                                                         <  2:410082/275‑1 (MQ=255)
                                                                                                                                 aacaacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCct                                                                                                                                               >  2:111665/1‑275 (MQ=255)
                                                                                                                                          cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCcaca                                                                                                                                                                                                                                                                               >  5:1062299/1‑138 (MQ=255)
                                                                                                                                                   gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGCCGTAAGTAATCCAACCCCACCTTCTCCATGTGGAGAGGGTGGGAttgg                                                                                                                                                                                                                                                             <  5:672174/147‑3 (MQ=38)
                                                                                                                                                    ttATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                              >  6:97800/1‑129 (MQ=255)
                                                                                                                                                      aTCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGC                                                                                                                                                                                                                                                         <  8:453476‑M1/148‑1 (MQ=255)
                                                                                                                                                         aCCAGCCGGACAAACGTACTGGAGGATTAGTCCGTCATCGATCGTCTGCCACACACCACCAGCGGGGATGAAGTGCTGGAACTGATGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGTGTTAATTGCCTGATGCGCTACGCTTAGCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTg                                                                                                                       <  1:86492/275‑1 (MQ=255)
                                                                                                                                                          ccAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGA                                                                                                                                                                                                                                                     >  6:1068834‑M1/1‑148 (MQ=255)
                                                                                                                                                            aGCGTGACACATGCACTGGATGAGGGATCCGTCATCGAGCGGCTGGCACAACCCACCCGCGTGGATGAAGTGCTGGAACTGCTGACAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGGGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTg                                                                                                                       <  2:434206/272‑1 (MQ=38)
                                                                                                                                                            aGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATg                                                                                                                                                                                                                                                   >  8:270243/1‑148 (MQ=255)
                                                                                                                                                             gCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccc                                                                                                                                                >  2:27191/1‑246 (MQ=255)
                                                                                                                                                             gCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAgg                                                                                                                    >  1:4466/1‑274 (MQ=255)
                                                                                                                                                                 gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGgagagg                                                                                                                                                                                                                                                                       >  7:787831/1‑121 (MQ=255)
                                                                                                                                                                 gACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCtctc                                                                                                                                            <  4:27190/246‑1 (MQ=255)
                                                                                                                                                                     aatgaacgtgtggatgaGTCCGGCATCGCGCGCCTGGCACACACCACCAGGTGTGATGAAGGGCTGGAACTGCTGGCNGGTCGTAAGTAATCCCATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAGTGCCAGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCt                                                                                                           <  3:117185/263‑1 (MQ=255)
                                                                                                                                                                     aaTGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCt                                                                                                           <  2:236956/275‑1 (MQ=255)
                                                                                                                                                                        gCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCAACCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCATAATTTGTTGAATTTGCCCGTTCTTGTAGGCAGGATAAGGCGCTTAAGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAGGCCTCAc                                                                                                        >  4:300843/1‑275 (MQ=255)
                                                                                                                                                                         cACTGGATGATGAGTCCGTCATCGAGCGTCTGTCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTatc                                                                                                                                                                                                                                      >  7:269123‑M2/1‑145 (MQ=255)
                                                                                                                                                                           cTGGATGATGATTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCcac                                                                                                                                                                                                                                                                                <  5:546062/104‑1 (MQ=255)
                                                                                                                                                                           cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   <  5:641905/117‑1 (MQ=255)
                                                                                                                                                                           cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTG                                                                                                                                                                                                                                                      <  2:557039‑M1/130‑1 (MQ=255)
                                                                                                                                                                           cTGGATGATCAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTCCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     <  5:915659/115‑1 (MQ=255)
                                                                                                                                                                            tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTTCCCCTCTCCTCACGGAGCGGGTTTGGGTGAGGGAAAAGCCTCACCCCa                                                                                                    >  2:22141/1‑275 (MQ=255)
                                                                                                                                                                             ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTg                                                                                                                                                                                                                                                                    >  7:720915/1‑114 (MQ=255)
                                                                                                                                                                             ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTatca                                                                                                                                                                                                                                     >  8:914439‑M2/1‑141 (MQ=255)
                                                                                                                                                                             ggATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGC                                                                                                                                                                                                                                                         <  8:131655‑M1/125‑1 (MQ=255)
                                                                                                                                                                              gatgatGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                                 >  7:641268/1‑148 (MQ=255)
                                                                                                                                                                               atgatgAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                                                                                                                                                            >  6:542499/1‑120 (MQ=255)
                                                                                                                                                                               atgatgAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTa                                                                                                                                                                                                                                              >  6:758301/1‑134 (MQ=255)
                                                                                                                                                                                 gatgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCC                                                                                                                                                                                                                                                        >  7:1115367‑M1/1‑122 (MQ=255)
                                                                                                                                                                                  atgaGTCCGTCATCGGGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                              <  1:30482/99‑1 (MQ=255)
                                                                                                                                                                                  atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      <  6:48478/107‑1 (MQ=255)
                                                                                                                                                                                  atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      >  5:48478/1‑107 (MQ=255)
                                                                                                                                                                                     aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   >  7:969801/1‑107 (MQ=255)
                                                                                                                                                                                     aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTatcagg                                                                                                                                                                                                                                   >  8:1036438‑M2/1‑133 (MQ=255)
                                                                                                                                                                                     aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGTAAAAGCCTCCCCCCAGccccccag                                                                                           >  4:41814/1‑270 (MQ=255)
                                                                                                                                                                                     aGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGC                                                                                                                                                                                                                                                         <  7:455295‑M1/117‑1 (MQ=255)
                                                                                                                                                                                      gTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          >  7:1062572/1‑131 (MQ=255)
                                                                                                                                                                                       tCCGTCATCGAGCGTCTGGCACACACCGCCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     <  6:79231/103‑1 (MQ=255)
                                                                                                                                                                                        ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                                                                                                                                                            <  6:219327/111‑1 (MQ=255)
                                                                                                                                                                                        ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGc                                                                                                                                                                                                                                           >  4:165332/1‑128 (MQ=255)
                                                                                                                                                                                           tCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATc                                                                                                                                                                                                                    >  5:758561/1‑148 (MQ=255)
                                                                                                                                                                                            cATGGAGCTGCTGGGACACACCACACGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCGCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCGTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTaaa                                                                                     <  4:84812/274‑1 (MQ=255)
                                                                                                                                                                                               cgggCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGGCGTAAGTAAGCCAATCCCCCCCTCGCCACATGGAGAAGGTGGGGTTAATTGCCTGGTGCGCTACGCTTATCAGGCCGACAGGAGGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAGAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAAc                                                                                  <  4:62499/271‑1 (MQ=255)
                                                                                                                                                                                                gAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATTCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGAT                                                                                                                                                                                                                                                    <  8:826063‑M1/111‑1 (MQ=255)
                                                                                                                                                                                                 aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTggg                                                                                                                                                                                                                                                                  >  6:669362/1‑96 (MQ=38)
                                                                                                                                                                                                 aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTa                                                                                                                                                                                                                                                              >  5:996699/1‑100 (MQ=255)
                                                                                                                                                                                                 aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                                                                                                                                                            >  8:495898/1‑102 (MQ=255)
                                                                                                                                                                                                 aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTa                                                                                                                                                                                                                                        <  7:360832‑M2/122‑2 (MQ=255)
                                                                                                                                                                                                    gTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                           <  6:350012/148‑1 (MQ=255)
                                                                                                                                                                                                        ggCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          >  6:63704/1‑113 (MQ=255)
                                                                                                                                                                                                            cacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     <  8:417780/82‑1 (MQ=255)
                                                                                                                                                                                                            cacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                                                                                                                                                             <  6:895945/122‑1 (MQ=255)
                                                                                                                                                                                                              caccaaCAGCGTGGGTGAAGTGCTTGAACTGCTGTCACGTCGTAAGTAATCCAATCCCGCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACTGGATGCATCACAAGTTGTTGAATTTGCACGTTCTTGTAGGCCGGATATGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCGCCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAAt                                                                  <  4:484061/275‑1 (MQ=255)
                                                                                                                                                                                                                ccaccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTG                                                                                                                                                                                                                                                          >  5:403502‑M1/1‑89 (MQ=255)
                                                                                                                                                                                                                  accaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCC                                                                                                                                                                                                                                                        <  5:885645‑M1/89‑1 (MQ=255)
                                                                                                                                                                                                                   ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTggg                                                                                                                                                                                                                                                                  >  5:588839/1‑78 (MQ=37)
                                                                                                                                                                                                                     aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTg                                                                                                                                                                                                                                                                    >  6:440088/1‑74 (MQ=255)
                                                                                                                                                                                                                     aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTg                                                                                                                                                                                                                                                                    <  5:440088/74‑1 (MQ=255)
                                                                                                                                                                                                                     aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGa                                                           <  2:77314/275‑1 (MQ=255)
                                                                                                                                                                                                                               gAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGAGCTACGCTTATCAGGCCTACAGGATGCATCACAAt                                                                                                                                                                                                               >  6:584507/1‑118 (MQ=255)
                                                                                                                                                                                                                                  gTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     >  8:643221/1‑61 (MQ=255)
                                                                                                                                                                                                                                      tGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCt                                                                                                           <  2:538568/211‑1 (MQ=255)
                                                                                                                                                                                                                                        gAACTGCTGGCAGGTTGTAAGTCATCCAATCCCACCCTCTCCCCATGGGGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGGTGAATTTGCACGGTCTTGTAGGCAGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGgc                                          <  4:328650/274‑1 (MQ=255)
                                                                                                                                                                                                                                         tgataagcgtagcgcatcaggcaattAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   <  8:1069646‑M2/30‑1 (MQ=255)
                                                                                                                                                                                                                                         aaCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGgcg                                         <  4:128791/274‑1 (MQ=255)
                                                                                                                                                                                                                                             gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGCAAAGCATCACCCAAGCCCTCCCGGGTAAAAACATTGATGAAGGGTAATACACTGAAAGCATAACATTTTGGAGCAGGt                                    >  4:178436/1‑275 (MQ=255)
                                                                                                                                                                                                                                             gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCCCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGt                                    >  4:101481/1‑275 (MQ=255)
                                                                                                                                                                                                                                               tgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTaa                                  >  2:455694/1‑275 (MQ=255)
                                                                                                                                                                                                                                                ggCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                     <  7:87221/111‑1 (MQ=255)
                                                                                                                                                                                                                                                  cAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATaa                                                                                                                                                                             >  7:1021264/1‑133 (MQ=255)
                                                                                                                                                                                                                                                       cGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttg                                                                                                                                                                                                            <  8:725222/97‑1 (MQ=255)
                                                                                                                                                                                                                                                             gaatccaaaccCCCCCTCCCCAAATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGc                     <  2:178437/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcg                                                                                                                                                                                                                                                 <  7:501148/51‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          >  8:537946/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                 ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          <  7:537946/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGaa                                                                                                                                                                                                       >  8:361288/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                          cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATaagcgtagcgcatcagg                                                                                                                                                                                                                                   <  5:193411‑M1/55‑18 (MQ=255)
                                                                                                                                                                                                                                                                          cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTatca                                                                                                                                                                                                                                     <  6:407194‑M2/53‑5 (MQ=255)
                                                                                                                                                                                                                                                                          cccTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAAct       >  2:307806/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACAtt                                                                               >  4:551559/1‑198 (MQ=255)
                                                                                                                                                                                                                                                                                  acaTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAGTTTGCACGTTCTTGTAGGCCGGATAAg                                                                                                                                                                            >  5:313074/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                                   caTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGc                                                                                                                                                                                                  >  7:409801/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                                   caTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCt                                                                                                                                                                                           <  8:598530/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   caTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCt                                                                                                                                                                                           >  7:598530/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                              
CACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGG  >  NC_000913/3773940‑3774484

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