Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 951,297 Δ1 bp 5.6% coding (128/2019 nt) cyoB → cytochrome bo terminal oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451951,2970A. 5.6% 252.4 / 7.0 71coding (128/2019 nt)cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base . (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GCGGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGTGATCATGATCATCTTCATGGCGATGCCGTTCT  >  AE015451/951072‑951524
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   
gcgGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGc                                                                                                                                                                                                                           <  1:482245/236‑1 (MQ=255)
                               aGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGa                                                                                                                                                                                            <  1:280967/236‑1 (MQ=255)
                                   cTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCg                                                                                                                                                                                         <  1:929120/235‑1 (MQ=255)
                                         acaACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGAt                                                                                                                                                                                  >  1:772436/1‑236 (MQ=255)
                                                   gtgtGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTgg                                                                                                                                                                                                                     >  2:271570/1‑191 (MQ=255)
                                                         aaTCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAACAAGATCGGCGTGATGTACATCAtcgtc                                                                                                                                                                     >  1:146899/1‑233 (MQ=255)
                                                               gTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGactac                                                                                                                                                                                                             >  1:825409/1‑187 (MQ=255)
                                                               gTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGactac                                                                                                                                                                                                             <  2:825406/187‑1 (MQ=255)
                                                                tATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTtgt                                                                                                                                                                                                                                >  2:630049/1‑167 (MQ=255)
                                                                 aTGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCg                                                                                                                                                                                         >  2:679714/1‑205 (MQ=255)
                                                                 aTGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCg                                                                                                                                                                                         <  1:679717/205‑1 (MQ=255)
                                                                     aaCCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGAc                                                                                                                                                                                                      >  2:787604/1‑188 (MQ=255)
                                                                     aaCCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAATGTCGGCGTGATGTACATCATCGTCGCGATGGTCAtg                                                                                                                                                         >  1:748085/1‑233 (MQ=255)
                                                                                 gTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCttt                                                                                                                                          <  1:393015/236‑1 (MQ=255)
                                                                                    cAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTg                                                                                                                                                                                                           <  2:148553/168‑1 (MQ=255)
                                                                                            tAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCcat                                                                                                                               <  1:50836/236‑1 (MQ=255)
                                                                                               gagaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGCGCTACTTGTGGAGCGAGTGGCTGACTACTg                                                                                                                                                                                                           >  2:255824/1‑157 (MQ=255)
                                                                                               gagaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTATTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAACAAGATCGGCGTGATGTACATCATCGTCGCGCTGGTCAtgctgct                                                                                                                                                    >  1:522563/1‑212 (MQ=255)
                                                                                               gagaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTAcatcat                                                                                                                                                                         >  2:417388/1‑191 (MQ=255)
                                                                                               gagaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTCTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGAt                                                                                                                                                                                  >  2:204280/1‑181 (MQ=255)
                                                                                               gagaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTCTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTc                                                                                                                                                            >  2:432329/1‑203 (MQ=255)
                                                                                                 gaTGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCatga                                                                                                                          >  1:475811/1‑236 (MQ=255)
                                                                                                       cGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGAGGTACATCAtcgtc                                                                                                                                                                     >  1:3146/1‑187 (MQ=255)
                                                                                                        ggTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTAc                                                                                                                   >  1:646121/1‑236 (MQ=255)
                                                                                                                     tcgAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCGTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCAtgctgct                                                                                                                                                    <  1:76461/188‑1 (MQ=255)
                                                                                                                       gAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTAc                                                                                                                   >  2:801835/1‑221 (MQ=255)
                                                                                                                          gCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTa                                                                                                                    >  1:452301/1‑217 (MQ=255)
                                                                                                                          gCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAAc                                                                                                 >  1:495232/1‑236 (MQ=255)
                                                                                                                          gCGATACCCTAGCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCGTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCg                                                                                                                                                                                         <  2:29101/148‑1 (MQ=255)
                                                                                                                               gcccTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGc                                                                                                                                        <  1:84862/191‑1 (MQ=255)
                                                                                                                                cccTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATg                                                                                                                                                                                 >  1:855656/1‑150 (MQ=255)
                                                                                                                                cccTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACcggcgg                                                                                           >  1:855145/1‑236 (MQ=255)
                                                                                                                                 ccTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAAcc                                                                                                <  1:238301/230‑1 (MQ=255)
                                                                                                                                           gCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGc                                                                                 >  1:270520/1‑235 (MQ=255)
                                                                                                                                           gCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGACGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGCGGCTGCCTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCttt                                                                                                                                          <  2:353704/178‑1 (MQ=255)
                                                                                                                                             cGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCaa                                                                                                  >  2:815959/1‑216 (MQ=255)
                                                                                                                                               aTAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTAcc                                                                             <  1:426782/235‑1 (MQ=255)
                                                                                                                                                  gTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCatg                                                                                                                           >  2:464694/1‑186 (MQ=255)
                                                                                                                                                    aaTGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgc                                                                       <  2:646118/235‑1 (MQ=255)
                                                                                                                                                      tGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgcc                                                                      <  1:911974/235‑1 (MQ=255)
                                                                                                                                                         tGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTAc                                                                              >  2:486348/1‑224 (MQ=255)
                                                                                                                                                               ttGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTAt                                                            <  1:255827/236‑1 (MQ=255)
                                                                                                                                                                tGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTAt                                                            >  1:541990/1‑235 (MQ=255)
                                                                                                                                                                       aTCGCGCTAGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGAACTACTTGTGGAGCGAGTGGCTGACTACTGTCg                                                                                                                                                                                                        <  2:662253/88‑1 (MQ=255)
                                                                                                                                                                         cgcgCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTCTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCtt                                                 <  1:432332/236‑1 (MQ=255)
                                                                                                                                                                          gcgcTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTCTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCtt                                                 <  1:204282/235‑1 (MQ=255)
                                                                                                                                                                             cTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCtt                                                 <  2:823832/233‑1 (MQ=255)
                                                                                                                                                                              tCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTATTTGTGGAGCGCTTGGCTGACTACTGTCGACCACATAAAGATCGGCGTGAt                                                                                                                                                                                  <  2:407730/103‑1 (MQ=255)
                                                                                                                                                                              tCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCg                                                                                                                                      <  2:758432/147‑1 (MQ=255)
                                                                                                                                                                              tCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCg                                                                                                                                      >  1:758435/1‑147 (MQ=255)
                                                                                                                                                                              tCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCAc                                              >  1:383929/1‑235 (MQ=255)
                                                                                                                                                                              tCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGGACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCCACCGGCGGCTCCg                                                                                      >  2:163720/1‑195 (MQ=255)
                                                                                                                                                                                ggcggTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGc                                                                                                                                                                                       >  1:383327/1‑96 (MQ=255)
                                                                                                                                                                                ggcggTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGc                                                                                                                                                                                       <  2:383324/96‑1 (MQ=255)
                                                                                                                                                                                ggcggTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCg                                            <  1:541953/235‑1 (MQ=255)
                                                                                                                                                                                   ggTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCCAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACCAAAAGATCg                                                                                                                                                                                         <  2:930916/91‑1 (MQ=255)
                                                                                                                                                                                          cTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCAcc                                             <  1:775286/224‑1 (MQ=255)
                                                                                                                                                                                           tGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGtga                                <  1:384578/236‑1 (MQ=255)
                                                                                                                                                                                             tcgtctGTCTTATCACCGATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgc                                                                       <  2:600962/195‑1 (MQ=255)
                                                                                                                                                                                             tcgtcgGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGgtgt                                  <  2:382978/232‑1 (MQ=255)
                                                                                                                                                                                             tcgtcgGGCTTATCACCTATTTCCGCAAGTGGACCTACTTGGGGGGCGAGGGGCTGACTACTGTCCACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGa                                                                                                                                                                 >  2:157078/1‑105 (MQ=255)
                                                                                                                                                                                                tcgGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTAt                                                            <  2:340672/203‑1 (MQ=255)
                                                                                                                                                                                                tcgGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGtgatcatg                           <  1:193746/236‑1 (MQ=255)
                                                                                                                                                                                                tcgGTCTTATCACCGATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGaa                                                                                    <  2:141153/179‑1 (MQ=255)
                                                                                                                                                                                                   gTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGc                                           <  2:382798/217‑1 (MQ=255)
                                                                                                                                                                                                            ccTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGTGATCATGATCATCTTCATg               >  1:492715/1‑236 (MQ=255)
                                                                                                                                                                                                              tATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGTCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGTGATCATGATCATCTTCATgg              <  1:877852/235‑1 (MQ=255)
                                                                                                                                                                                                                ttgccGCAAGTGGACCTACGTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGCTGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgc                                                                       <  2:239093/173‑1 (MQ=255)
                                                                                                                                                                                                                      cAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgcc                                                                      >  2:400609/1‑171 (MQ=255)
                                                                                                                                                                                                                      cAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTgccgcc                                                                      <  1:400612/171‑1 (MQ=255)
                                                                                                                                                                                                                         gTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGTGATCATGATCATCTTCATGGCGATGCCGttct  <  1:699036/236‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   
GCGGTCGTAGCCACGAAGAAGCAGGCAGCAAAGATCTGGCCACAACCAAGGGTGTGGAATCGAGTATGCAACCAGCTGCCGGTGCAGAGGAGTAAGAGATGTTCGGTAAATTAAGCCTGGAGGCGATACCCTATCACGAGCCGATAGTCATGGTGACGCTTGCCATGATCGCGCTCGGCGGTATCGCTGTCGTCGGTCTTATCACCTATTTCCGCAAGTGGACCTACTTGTGGAGCGAGTGGCTGACTACTGTCGACCACAAAAAGATCGGCGTGATGTACATCATCGTCGCGATGGTCATGCTGCTGCGCGGCTTTGCCGACGCCATCATGATGCGTACCCAGCTGGCTGCCGCAACCGGCGGCTCCGAAGGCTACCTGCCGCCTGAACACTATGACCAGATCTTCACCGCTCACGGTGTGATCATGATCATCTTCATGGCGATGCCGTTCT  >  AE015451/951072‑951524

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