Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,215 A→G 5.5% intergenic (+21/‑209) 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,2150AG 5.5% 1913.4 / 62.8 569intergenic (+21/‑209)mtlA/mtlDmannitol‑specific PTS enzyme: IIA, IIB and IIC components/mannitol‑1‑phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (259/274);  new base G (14/18);  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 = 1.42e-06

Only 100 of 596 total aligned reads displayed.

ACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGG  >  NC_000913/3773941‑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)
                                  aaCGGGCGTCGTGTTCGGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCCGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTa                                                                                                                                                                                                                                              <  4:101254/274‑1 (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)
                                                                                                                            ccGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATCCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACTTTCTTGTATGCCGCATAAGGCGCTTACGCCCAAAGCGGCGCTGc                                                                                                                                                   >  4:176560/1‑275 (MQ=255)
                                                                                                                                aacaacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCct                                                                                                                                               >  2:111665/1‑275 (MQ=255)
                                                                                                                                     cGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTctcct                                                                                                                                          <  4:159014/275‑1 (MQ=255)
                                                                                                                                         cacaTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCcaca                                                                                                                                                                                                                                                                               >  5:1062299/1‑138 (MQ=255)
                                                                                                                                              tCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      >  6:987444/1‑142 (MQ=255)
                                                                                                                                                  gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGCCGTAAGTAATCCAACCCCACCTTCTCCATGTGGAGAGGGTGGGAttgg                                                                                                                                                                                                                                                             <  5:672174/147‑3 (MQ=38)
                                                                                                                                                  gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGCCGTAAGTAATCCAACCCCACCTTCTCCATGTGGAGAGGGTGGGAttgg                                                                                                                                                                                                                                                             >  6:672174/1‑145 (MQ=38)
                                                                                                                                                     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)
                                                                                                                                                                          cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   >  6:641905/1‑117 (MQ=255)
                                                                                                                                                                          cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTgg                                                                                                                                                                                                                                                                   <  5:641905/117‑1 (MQ=255)
                                                                                                                                                                          cTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTG                                                                                                                                                                                                                                                      <  2:557039‑M1/130‑1 (MQ=255)
                                                                                                                                                                          cTGGATGATCAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTCCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     <  5:915659/115‑1 (MQ=255)
                                                                                                                                                                           tGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTatcagg                                                                                                                                                                                                                                   >  6:897892‑M2/1‑142 (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)
                                                                                                                                                                                 atgaGTCCGTCATCGGGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                              <  1:30482/99‑1 (MQ=255)
                                                                                                                                                                                 atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      >  5:48478/1‑107 (MQ=255)
                                                                                                                                                                                 atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      <  6:48478/107‑1 (MQ=255)
                                                                                                                                                                                   gaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAgg                                                                                                                                                                                                                           <  7:270243/148‑1 (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)
                                                                                                                                                                                      tCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAt                                                                                                                                                                                                                                                                              >  7:1068231/1‑94 (MQ=255)
                                                                                                                                                                                       ccGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                                                                                                                                                            <  6:219327/111‑1 (MQ=255)
                                                                                                                                                                                         gTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                                                                                                                                                                                                               >  6:871620/1‑138 (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)
                                                                                                                                                                                                aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTa                                                                                                                                                                                                                                                              >  5:996699/1‑100 (MQ=255)
                                                                                                                                                                                                aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                                                                                                                                                            >  8:495898/1‑102 (MQ=255)
                                                                                                                                                                                                aGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAAtt                                                                                                                                                                                                              <  5:562429/148‑1 (MQ=255)
                                                                                                                                                                                                   gTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                                                                                                                                                           <  6:350012/148‑1 (MQ=255)
                                                                                                                                                                                                     cTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGGAAGTAATCCAATCCCACCCTATCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCGCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATtgatg                                                                           <  4:123915/274‑1 (MQ=255)
                                                                                                                                                                                                           cacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                                                                                                                                                     <  8:417780/82‑1 (MQ=255)
                                                                                                                                                                                                             caccaaCAGCGTGGGTGAAGTGCTTGAACTGCTGTCACGTCGTAAGTAATCCAATCCCGCCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACTGGATGCATCACAAGTTGTTGAATTTGCACGTTCTTGTAGGCCGGATATGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCGCCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAAt                                                                  <  4:484061/275‑1 (MQ=255)
                                                                                                                                                                                                               ccaccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAttt                                                                                                                                                                                                    >  7:1128804/1‑143 (MQ=255)
                                                                                                                                                                                                               ccaccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTG                                                                                                                                                                                                                                                          >  5:403502‑M1/1‑89 (MQ=255)
                                                                                                                                                                                                                 accaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCC                                                                                                                                                                                                                                                        <  5:885645‑M1/89‑1 (MQ=255)
                                                                                                                                                                                                                  ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTggg                                                                                                                                                                                                                                                                  <  6:588839/78‑1 (MQ=37)
                                                                                                                                                                                                                  ccaGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAgg                                                                                                                                                                                                                                                                      >  7:1022594/1‑74 (MQ=255)
                                                                                                                                                                                                                    aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTg                                                                                                                                                                                                                                                                    >  6:440088/1‑74 (MQ=255)
                                                                                                                                                                                                                    aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTg                                                                                                                                                                                                                                                                    <  5:440088/74‑1 (MQ=255)
                                                                                                                                                                                                                    aGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGa                                                           <  2:77314/275‑1 (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)
                                                                                                                                                                                                                                            gctgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcg                                                                                                                                                                   <  6:758561/148‑1 (MQ=255)
                                                                                                                                                                                                                                              tgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTaa                                  >  2:455694/1‑275 (MQ=255)
                                                                                                                                                                                                                                               ggCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAtt                                                                                                                                                                                                     <  7:87221/111‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttg                                                                                                                                                                                                            <  8:725222/97‑1 (MQ=255)
                                                                                                                                                                                                                                                            gaatccaaaccCCCCCTCCCCAAATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGc                     <  2:178437/273‑1 (MQ=255)
                                                                                                                                                                                                                                                               tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATgcg                                                                                                                                                                                                                                                 <  7:501148/51‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          <  7:537946/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCt                                                                                                                                                                                                                                          >  8:537946/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                     cccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGaa                                                                                                                                                                                                       >  8:361288/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                     cccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGaa                                                                                                                                                                                                       <  7:361288/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAAttt                                                                                                                                                                                                    <  5:589804/88‑1 (MQ=25)
                                                                                                                                                                                                                                                                         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                                                                                                                                                                                           >  7:598530/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                                  caTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCt                                                                                                                                                                                           <  8:598530/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                             
ACGGTACGGTTGAAGCGAAAGATCGCGTACTGAAAACGGGCGTCGTGTTCTGCCAGTACCCGGAAGGCGTGCGCTTCGGTGAAGAAGAAGATGACATTGCCCGTCTGGTGATTGGTATTGCTGCCCGTAACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACGGAGAGGGTTTGGGTGAGGGAAAAGCCTCACCCCAGCCCTCTCGGGTAAAAACATTGATGAAGGTTAATACTATGAAAGCATTACATTTTGGCGCAGGTAATATCGGTCGTGGCTTTATCGGTAAACTGCTGG  >  NC_000913/3773941‑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: