Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 339,769 G→A 6.3% intergenic (+26/‑396) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,7690GA 6.3% 380.8 / 11.9 124intergenic (+26/‑396)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (56/57);  new base A (3/4);  total (62/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 131 total aligned reads displayed.

GTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGC  >  NC_000913/339499‑340037
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                            
                             aaCCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAg                                                                                                                                                                                                                                             <  4:300740/275‑1 (MQ=255)
                                   aaCGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATCAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGACCGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCCTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                                                                                                                                                                                                       >  2:15934/1‑275 (MQ=255)
                                   aaCGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCCTCTGACCTACATGTTTATGCCGGATGCGGCCTGAAAGCCTTACCCGGCCTACCAAAGGGGCaaagaccgt                                                                                                                                                                                                                                       >  4:364417/1‑268 (MQ=255)
                                         aCCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                                                                                                                                                                  <  2:364418/274‑1 (MQ=255)
                                         aCCATTCCGCTGCTCGGGTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGGGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                                                                                                                                                                  <  1:160884/274‑1 (MQ=255)
                                          ccATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTATCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGGTGCGGGCGTCGGCATGATTGGCGCCGGCATTGTGCATCCACCTTTTGCCTTCTTCGAGGAAGCCATTCTTGGCTGGTGCGAACGTGACGGCGTCCGACCTACATGTTCATTCCGGAAGCGGCGTGGCCCCCTTGTCCCCCCTACCATGGGCGCACGAAACTGCAggacc                                                                                                                                                                                                                                >  4:146278/1‑272 (MQ=255)
                                              tCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCCCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTCCGGGGTCTGGCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTTCCCGGCCTACGAAAGTGGCAAGAAACTGTAGNGCTgca                                                                                                                                                                                                                             >  4:512200/1‑272 (MQ=255)
                                                     ctCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCACTGCGCATCCACCTTTTGCCTTCTTCGAGACAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTTGTACCTACCTGTTCATGCCGGCTGCGGCCGCAACCCCTTATCCGGCCTACGAATGGCGCCAGAATCTCGACGCCTGATAAGCGt                                                                                                                                                                                                                       >  4:200178/1‑273 (MQ=255)
                                                          ttttATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAt                                                                                                                                                                                                                <  2:512201/275‑1 (MQ=255)
                                                           tttATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTGGCAGCGGCAGTCTGCCGGGGATCAGCACCGCCATCGCCCATAAAGAGGCGTGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCTTGCTTCGAGAAAGCCAGTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATc                                                                                                                                                                                                               <  2:309329/275‑1 (MQ=255)
                                                               tGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGAAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                                                                                                                                            >  2:371062/1‑274 (MQ=255)
                                                                    gTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGTGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTg                                                                                                                                                                                                      >  2:7416/1‑275 (MQ=255)
                                                                            ggCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAATGGCGCGAGAATTGTAGGCCTTATAAGCGTAGTGCATCAGGCagtctggcgttgg                                                                                                                                                                                              >  2:512127/1‑262 (MQ=255)
                                                                                 tCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          >  2:471327/1‑210 (MQ=255)
                                                                                     cATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGaa                                                                                                                                                                                                                                                      <  4:471326/210‑1 (MQ=255)
                                                                                      aTCACCCGCGTGGGGAGCAGCGGCATTCTGCCGGTGATCAACCCCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATGGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGTCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgcg                                                                                                                                                                                    <  2:212681/275‑1 (MQ=255)
                                                                                            cgcgGGGGTAGGAGCGGCATTTTGCCGGGGATCAACACCGCCATCGCCCATAAAGATGGGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTAAATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                                                                                                               <  4:349996/274‑1 (MQ=255)
                                                                                            agcgTTGGTAGCAGCGGCAGTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                                                                                                               <  1:148988/273‑1 (MQ=255)
                                                                                               gTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGGTTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAAc                                                                                                                                                                            >  2:99173/1‑274 (MQ=255)
                                                                                                       cagcGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                                   >  4:368846/1‑275 (MQ=255)
                                                                                                            gCATTCTGCCGGTGATCAACACCGCCATCGCCCATAATGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCCTTGATAAAGCCACTCTTGGCTGGTGCGAACGTTACGGCGGCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCATGCACCGGTAGGCCTGATAAGCGTACCCCATCCGGCTTTTGCCCCACTTTCCCGATCCCGCGTGTACGGCTTATCCGGCCt                                                                                                                                                              >  4:166234/1‑275 (MQ=37)
                                                                                                            gCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgcgg                                                                                                                                                                                   >  2:427131/1‑254 (MQ=255)
                                                                                                                tttGCCCGTGATCAACCCCGCCCTCGTCCCTAAAGAGTCGGGCGTCGGCCTGATTGGCGCGGGCATTGTTCAGCCACCGTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                                           <  4:414619/272‑1 (MQ=25)
                                                                                                                tCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                                                                                                               <  4:427130/254‑1 (MQ=255)
                                                                                                                  tGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAACCCGCAGGCCTGATAAGCGTAACGCATCGAGCATTTGTCACCATTGCCGGGTGGGGCGTGAACGCCTTTACCGGCCTaacagt                                                                                                                                                        >  4:236834/1‑270 (MQ=38)
                                                                                                                         gATCCACTCCGCCATCGCCCATGAAGATTCGGGCGTCGGGATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGCATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACGAATGGCGCa                                                                                                                                                  <  1:58894/274‑1 (MQ=255)
                                                                                                                         gATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcg                                                                                                                                                    <  2:15285/272‑1 (MQ=255)
                                                                                                                          aTCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGt                                                                                                                                                                                                                       >  4:571905/1‑204 (MQ=255)
                                                                                                                           tCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGGTTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                                                                                                                                               <  2:200845/275‑1 (MQ=255)
                                                                                                                            cAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTGGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                                                                                                                                               <  2:217075/274‑1 (MQ=255)
                                                                                                                               cacCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                                                                                                                                          <  5:531834/148‑1 (MQ=255)
                                                                                                                                acCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGAAc                                                                                                                                                                                                                                                                         <  8:205750/148‑1 (MQ=255)
                                                                                                                                      accGCCCATAAAGATGCTGGCGTCGGCATGATTGGCGCGTGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGTCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                                                                     <  2:483941/272‑1 (MQ=255)
                                                                                                                                       tCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGACAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTCGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCATTTTCCGGCCTCCGAACGGCGCAAGAATCTGAACGccg                                                                                                                                   >  2:210004/1‑274 (MQ=255)
                                                                                                                                          cccATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGATAAAGCCATTCTTGCCTGGCGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAACGATCAGTAGGCCTgct                                                                                                                                >  1:95292/1‑273 (MQ=255)
                                                                                                                                            cATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATg                                                                                                                                                                                                                                                                            >  6:657692/1‑133 (MQ=255)
                                                                                                                                            cATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGaa                                                                                                                                                                                                                                                                          >  6:267484/1‑135 (MQ=255)
                                                                                                                                            cATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGaa                                                                                                                                                                                                                                                                          <  5:267484/135‑1 (MQ=255)
                                                                                                                                              tAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                                   >  6:384288/1‑140 (MQ=255)
                                                                                                                                              tAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                                                                                                                                                   <  5:384288/140‑1 (MQ=255)
                                                                                                                                                     gCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACATTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          >  2:307227/1‑142 (MQ=255)
                                                                                                                                                         gCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACATTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGaa                                                                                                                                                                                                                                                      <  4:307226/142‑1 (MQ=255)
                                                                                                                                                            tCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                                                                                                                                                                     <  3:26745/156‑1 (MQ=255)
                                                                                                                                                              ggCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACCAATGGCGCCAGAATCTGTAGGCCTGATAAGCGTACCCCATCCGGCAt                                                                                                            >  4:109276/1‑275 (MQ=255)
                                                                                                                                                               gCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                                                                                                                                      >  5:425114/1‑120 (MQ=255)
                                                                                                                                                                      tGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                                                                                                                                        >  7:1020912/1‑111 (MQ=255)
                                                                                                                                                                      tGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                                                                                                                                        <  8:1020912/111‑1 (MQ=255)
                                                                                                                                                                      gggggcgGGAATTGTGCACCCACCTTTTGCCTGCTTCGGGAAAGCTATTCTTGGATGGGGCGAACGTTACGGGGTCTGACGTACATGTTCATGCCGGATGCGGGTTGAACGCCTTATCCGGCCTACGAATGGGGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccac                                                                                                     <  4:204985/270‑1 (MQ=18)
                                                                                                                                                                                  tgtgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          >  5:896628/1‑113 (MQ=255)
                                                                                                                                                                                  tgtgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          <  6:896628/113‑1 (MQ=255)
                                                                                                                                                                                   gtgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcg                                                                                       >  2:88400/1‑275 (MQ=255)
                                                                                                                                                                                    tgCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCTGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGCTgcgg                                                                                      >  1:22220/1‑275 (MQ=255)
                                                                                                                                                                                     gCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATg                                                                                                                                                                                                                                                    <  8:1047304/116‑1 (MQ=255)
                                                                                                                                                                                     gCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATg                                                                                                                                                                                                                                                    >  7:1047304/1‑116 (MQ=255)
                                                                                                                                                                                       aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          >  7:793681/1‑108 (MQ=255)
                                                                                                                                                                                       aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                                                                                                                          <  8:793681/108‑1 (MQ=255)
                                                                                                                                                                                       aTCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCACGCCGCATCCGGCATGa                                                                                                                                                                                                                                                                           <  7:621426/91‑1 (MQ=255)
                                                                                                                                                                                        tCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGt                                                                                   >  2:168284/1‑274 (MQ=255)
                                                                                                                                                                                         ccACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              <  8:1024447/86‑1 (MQ=255)
                                                                                                                                                                                         ccACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              >  7:1024447/1‑86 (MQ=255)
                                                                                                                                                                                         ccACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTg                                                                                  <  2:368847/274‑1 (MQ=255)
                                                                                                                                                                                         ccACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCCGCGTGAACGCCTTATCCGGCCTACGCATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCTCATCAAGCCTTTCCCACCACTGCCGGATGCGGCGtt                                                                                  >  2:420544/1‑273 (MQ=255)
                                                                                                                                                                                           aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                                                                                                                              <  5:439136/100‑1 (MQ=255)
                                                                                                                                                                                           aCCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                                                                                                                              >  6:439136/1‑100 (MQ=255)
                                                                                                                                                                                                ttGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                                                                                                                                         <  7:67107/148‑1 (MQ=255)
                                                                                                                                                                                                     tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                                  <  5:214636/86‑1 (MQ=255)
                                                                                                                                                                                                     tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                                                                                                                                                                  >  6:214636/1‑86 (MQ=255)
                                                                                                                                                                                                     tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCa                                                                                                                                                                                                                                               <  8:743984/105‑1 (MQ=255)
                                                                                                                                                                                                     tGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGt                                                                                                                                         >  4:121188/1‑207 (MQ=255)
                                                                                                                                                                                                        ttCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              <  5:915692/71‑1 (MQ=255)
                                                                                                                                                                                                        ttCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              >  6:915692/1‑71 (MQ=255)
                                                                                                                                                                                                         tCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                                                                                                                                                                     >  8:506611/1‑111 (MQ=255)
                                                                                                                                                                                                         tCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                                                                                                                                                                     <  7:506611/111‑1 (MQ=255)
                                                                                                                                                                                                         tCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTa                                                                                                                                                                                                                      <  5:611576/126‑1 (MQ=255)
                                                                                                                                                                                                         tCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGc                                                                                                                                     <  2:121189/207‑1 (MQ=255)
                                                                                                                                                                                                           gagaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTa                                                                                                                                                                                                                      >  8:828883/1‑124 (MQ=255)
                                                                                                                                                                                                           gagaAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTa                                                                                                                                                                                                                      <  7:828883/124‑1 (MQ=255)
                                                                                                                                                                                                              aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              <  7:254883/65‑1 (MQ=255)
                                                                                                                                                                                                              aaaGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATgcggcg                                                                                                                                                                                                                                                                              >  8:254883/1‑65 (MQ=255)
                                                                                                                                                                                                                 gCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                                                                                                                                          <  6:463266/66‑1 (MQ=255)
                                                                                                                                                                                                                   cATTCTTGGCTGGTGCGAACGTTACGGCGTCGGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGgcg                                                       <  1:121088/275‑1 (MQ=255)
                                                                                                                                                                                                                        ttGGCTGGTGCGAACGTTACGGCGTCTGACGTACATGTTCAGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCgaga                                                  <  2:520507/275‑1 (MQ=255)
                                                                                                                                                                                                                        ttGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCgaga                                                  <  2:129116/275‑1 (MQ=255)
                                                                                                                                                                                                                          ggCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGACTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCCTTTGCCACCACTGCCg                                                                                             >  2:256970/1‑230 (MQ=255)
                                                                                                                                                                                                                           gCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGt                                                                                                                                                                                                                                                                             <  6:543586/53‑1 (MQ=255)
                                                                                                                                                                                                                           gCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGt                                                                                                                                                                                                                                                                             >  5:543586/1‑53 (MQ=255)
                                                                                                                                                                                                                            cTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTTTAGGCCTGATAAGGGTAGCGCAACAGGCATTTGCCACCCATTCCGGAATCGGACTGGACGCCTTATCCGGCCTACCAGTGGCCCGAGAATc                                               >  4:92626/1‑274 (MQ=17)
                                                                                                                                                                                                                              gttggcaaCGTTACGTCGTCTGTCCTACAGGTTCCTGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATg                                                                                         <  4:256969/224‑1 (MQ=255)
                                                                                                                                                                                                                               gTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTa                                           <  2:129456/275‑1 (MQ=255)
                                                                                                                                                                                                                                ttcGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTGTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTa                                           <  2:410629/272‑1 (MQ=255)
                                                                                                                                                                                                                                tGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATTGCGCAAGGATCTGTAGGCCTGCTAAGGGTAGCGCATCAGGACTTTTGCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCCAGAGAATATGTaa                                          >  4:509696/1‑274 (MQ=17)
                                                                                                                                                                                                                                         tACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCa                                                                                                                                                                                                              <  7:855548/102‑1 (MQ=255)
                                                                                                                                                                                                                                              cGTATGACCTATCTGTTCATTCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGAAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGc                                        <  1:4385/263‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTTCCGGATTCGGCGTGGACGCCTTAACCGGCCTACGGGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAgcgc                        >  2:292436/1‑275 (MQ=255)
                                                                                                                                                                                                                                                        aCATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                                 <  6:332342/132‑1 (MQ=255)
                                                                                                                                                                                                                                                        aCAGGTTCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAAGCTGTAGGCATGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAgg                  <  2:227104/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccactGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCgaga                                                  <  4:217979/234‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACg                                                                                                                                                                           >  5:925353/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                       aTGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGCTGCGGCGTGGACCCCTTATCCGGCCTACGAGTGGCGCGAGGATCTTGAGGCCTGATAAGCGGAGCGCAGACGGCATTTGTCACCaatt   >  4:45238/1‑272 (MQ=14)
                                                                                                                                                                                                                                                                         ggggcgTAAACGCCTGTGCCGGCCCACGAACGGCGCAAGGATCTGTAGGCCTGGTGAGCGGCGGGCATCCGGCATGTGGCGCCATGGCCGGACGCGGCTGTAACGCCTTGGCCGGCCTACGAGTGGCGCAAGAATCTGTAGGCCTGGTAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGGGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGc  <  4:303131/272‑1 (MQ=14)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                            
GTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGGCGTCGGCATGATTGGCGCGGGCATTGTGCATCCACCTTTTGCCTGCTTCGAGAAAGCCATTCTTGGCTGGTGCGAACGTTACGGCGTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGC  >  NC_000913/339499‑340037

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