Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,739,789 A→C 6.1% E240A (GAA→GCA)  hutF → probable formiminoglutamate deiminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,739,7890AC 6.1% 227.1 / 7.2 66E240A (GAA→GCA) hutFprobable formiminoglutamate deiminase
Reads supporting (aligned to +/- strand):  ref base A (36/26);  new base C (2/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-02

GCTTTGGCGGCCAGGCGCCGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTCCCGGCAGTGGATTTCCTGGCACAGGGC  >  AE015451/5739558‑5740009
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            
gCTTTGGCGGCCAGGCGCCGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAgaaggaa                                                                                                                                                                                                                             >  2:1062762/1‑233 (MQ=255)
       cggcCAGGCGCCGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGCAGTCgtcgac                                                                                                                                                                                                                    >  1:530959/1‑235 (MQ=255)
        ggcCAGGCGCCGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCgacgac                                                                                                                                                                                                                    >  2:937594/1‑234 (MQ=255)
                  cGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGCAGTCGTCGACTtccttctgct                                                                                                                                                                                                         <  2:530958/235‑7 (MQ=255)
                   gAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTgcctggcctg                                                                                                                                                                                                         <  1:185460/234‑1 (MQ=255)
                       gacgGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCgg                                                                                                                                                                                                   <  1:613418/236‑1 (MQ=255)
                                    tcTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGcctggcctgg                                                                                                                                                                                                        <  1:232085/217‑1 (MQ=255)
                                        catcaACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGCGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTgcctg                                                                                                                                                                                                              >  1:135200/1‑208 (MQ=255)
                                                 cccggAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGcctggcctgg                                                                                                                                                                                                        <  1:897008/201‑1 (MQ=255)
                                                      aaCAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCac                                                                                                                                                                    >  1:338506/1‑236 (MQ=255)
                                                           tACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGc                                                                                                                                                                      <  1:907219/229‑1 (MQ=255)
                                                            aCCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtgg                                                                                                                                                                >  2:1119658/1‑234 (MQ=255)
                                                              ctgcaactgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCGCTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGAGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTgcc                                                                                                                                                                                             >  1:258132/1‑202 (MQ=255)
                                                              ctgcaactgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTcgacga                                                                                                                                                                                                                     >  1:541011/1‑179 (MQ=255)
                                                              ctgcaactgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTgcctggcctg                                                                                                                                                                                                         >  1:408968/1‑191 (MQ=255)
                                                              ctgcaactgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCac                                                                                                                                                                    >  1:989333/1‑228 (MQ=255)
                                                               tgcaactgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtggacgtg                                                                                                                                                           <  1:641112/236‑1 (MQ=255)
                                                                   actgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtg                                                                                                                                                                 >  2:911776/1‑226 (MQ=255)
                                                                     tgcaaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAg                                                                                                                                                                                   <  2:526534/206‑1 (MQ=255)
                                                                        aaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGCAGTCGTCGACTtcctt                                                                                                                                                                                                              <  2:551311/176‑2 (MQ=255)
                                                                        aaGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGCAGTCGTCGACTtcctt                                                                                                                                                                                                              >  1:551312/1‑175 (MQ=255)
                                                                          gcgcAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtggacgtg                                                                                                                                                           <  1:542596/225‑1 (MQ=255)
                                                                            gcAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGt                                                                                                                                              >  1:492237/1‑236 (MQ=255)
                                                                              aaCTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGc                                                                                                                                                                                                  <  1:829831/182‑1 (MQ=255)
                                                                                 tctTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTg                                                                                                                                         <  1:495381/236‑1 (MQ=255)
                                                                                  ctTCCCGCTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTgcctggcctg                                                                                                                                                                                                         <  2:492702/171‑1 (MQ=255)
                                                                                          tGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGc                                                                                                                                 <  1:1119659/235‑1 (MQ=255)
                                                                                               gCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTaca                                                                                                                                    <  1:587032/227‑1 (MQ=255)
                                                                                                ccGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTgg                                                                                                                                        >  2:649275/1‑222 (MQ=255)
                                                                                                  gcgcAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCCGAAGGAAGCCGACGACCGCCTGGCCTGGAg                                                                                                                                                                                                      >  2:163611/1‑158 (MQ=255)
                                                                                                  gcgcAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTacc                                                                                                                                    >  2:396747/1‑223 (MQ=255)
                                                                                                              cagcaCCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTtcctgg                                                                                                                                                                                                             >  2:493277/1‑139 (MQ=255)
                                                                                                                        gCCTGGGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCt                                                                                                                                          <  1:95295/196‑1 (MQ=255)
                                                                                                                         ccTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCg                                                                                                                 >  1:872985/1‑220 (MQ=255)
                                                                                                                          cTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtg                                                                                                                                                                 >  1:248744/1‑171 (MQ=255)
                                                                                                                            gtgCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGGACACGCTACCCATGCCGAACCgta                                                                                                               >  1:386278/1‑217 (MQ=255)
                                                                                                                              gCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGcc                                                                                            >  1:716577/1‑236 (MQ=255)
                                                                                                                                 tCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAg                                                                                                          <  1:433773/219‑1 (MQ=255)
                                                                                                                                  ccACTCTCGGCGGGCGGTGCGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACGGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTg                                                                                                                                                                                               <  2:263585/133‑1 (MQ=255)
                                                                                                                                         ctGCGGGCGGTGAGGCCGGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTg                                                                                                                                                <  1:941865/173‑1 (MQ=255)
                                                                                                                                          tGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAgcggcgg                                                                                 <  1:502280/235‑2 (MQ=255)
                                                                                                                                          tGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGcgcg                                                                                >  1:1092122/1‑236 (MQ=255)
                                                                                                                                          tGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGcgcg                                                                                >  1:657828/1‑236 (MQ=255)
                                                                                                                                                cggTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTAccc                                                                                                                            <  1:359419/186‑1 (MQ=255)
                                                                                                                                                cggTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAgcggcggtccacgt                                                                          <  1:9201/236‑9 (MQ=255)
                                                                                                                                                  gTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCg                                                                         <  1:219870/235‑1 (MQ=255)
                                                                                                                                                     aGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGtt                                                                     >  1:14892/1‑236 (MQ=255)
                                                                                                                                                     aGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGCGGCt                                                                                                                                                                              >  2:345441/1‑131 (MQ=255)
                                                                                                                                                               aGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGa                                                                                                               >  2:407598/1‑184 (MQ=255)
                                                                                                                                                                gATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTggc                                                            >  1:812444/1‑232 (MQ=255)
                                                                                                                                                                    gCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAAc                                                                                                                   >  1:709993/1‑175 (MQ=255)
                                                                                                                                                                    gCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAg                                                                                       >  2:834096/1‑203 (MQ=255)
                                                                                                                                                                    gCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGaccac                                                         >  2:947020/1‑233 (MQ=255)
                                                                                                                                                                    gCCGAAGGGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCCCATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCt                                                                                                                                          <  1:163890/152‑1 (MQ=255)
                                                                                                                                                                     ccGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCacgtgga                                                                                                                                                               >  2:89982/1‑130 (MQ=255)
                                                                                                                                                                       gcgggttgcgcggcCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGAACGAGCACGCGGACGAGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAgcggc                                                                                   >  1:244117/10‑204 (MQ=255)
                                                                                                                                                                         aGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTCCGAGCccgtgga                                                                                                                                                               >  2:137082/1‑126 (MQ=255)
                                                                                                                                                                         aGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGGGCATATCCAAATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTgc                                                                                                                                                                                              >  2:274746/1‑95 (MQ=255)
                                                                                                                                                                                gcggcCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTgg                                           >  1:704831/1‑235 (MQ=255)
                                                                                                                                                                                  ggcCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTc                                                                             >  1:120305/1‑199 (MQ=255)
                                                                                                                                                                                   gcCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTgg                                           <  2:492236/232‑1 (MQ=255)
                                                                                                                                                                                               cAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTccc                           <  1:498408/236‑1 (MQ=255)
                                                                                                                                                                                                                   caTCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGc                                                                                                                                                                                     <  2:166912/62‑1 (MQ=255)
                                                                                                                                                                                                                      cGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCCGCGGCGCGGTCGCCGGGTTGTGcc                                                                >  2:362246/1‑176 (MQ=255)
                                                                                                                                                                                                                      cGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTCCCGGCAgg                     >  1:297414/1‑218 (MQ=255)
                                                                                                                                                                                                                         cGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTCCCGGCAGTGGATTTCCTGGCACAgggc  <  1:1062763/235‑1 (MQ=255)
                                                                                                                                                                                                                              agcagAAGGAAGTCGACGACTGCCTGGCATGGAGCGGCCTGCGCCCGATTCAGTGGCTGTAGGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCg                                                                                                                 >  2:799317/1‑119 (MQ=255)
                                                                                                                                                                                                                                   aaggaagTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTCCCGGCAGTGGATTTCCTGGca        >  1:250871/1‑219 (MQ=255)
                                                                                                                                                                                                                                      gaagTCGCCGGCGGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAg                                                                                                          <  2:309522/118‑1 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                            
GCTTTGGCGGCCAGGCGCCGAACGACGGGCAGCGGCGCTTCATCAACTCCACTGAACAGTACCTGCAACTGCAAGCGCAACTCTTCCCGTTGCTGGCCGCGCAACCCGCGCAGCACCTGGGCCTGTGCTTCCACTCTCTGCGGGCGGTGAGGCCTGGGCAGATCGCCGAAGTGCTGGCGGCCAGCGACAGGCAGTGCCCGGTGCATATCCACATCGCCGAACAGCAGAAGGAAGTCGACGACTGCCTGGCCTGGAGCGGCCTGCGCCCGCTGCAGTGGCTGTACGAGCACGTGGACGTGGACCGCCGCTGGTGCCTGGTACACGCTACCCATGCCGAACCGGACGAAGTGACTGCCATGGCCCGCAGCGGCGCGGTCGCCGGGTTGTGCCTGACCACCGAAGCCAACCTGGGCGATGGGATCTTCCCGGCAGTGGATTTCCTGGCACAGGGC  >  AE015451/5739558‑5740009

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