Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,790,473 A→G 8.2% N135D (AAC→GAC)  dxr → 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,790,4730AG 8.2% 233.7 / 5.9 69N135D (AAC→GAC) dxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base G (4/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGCGTCATCGTCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCAAGCCT  >  AE015451/1790241‑1790707
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
cgtgcgtCATCGTCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCcaaca                                                                                                                                                                                                                                          >  1:845459/1‑236 (MQ=255)
 gtgcgtCATCGTCCGGCCTTCGCTGTGGTGCCGAGCCCCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCcaacaa                                                                                                                                                                                                                                         <  1:292266/236‑1 (MQ=255)
  tgcgtCATCGTCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGGGTGGCTGGCCAACAAg                                                                                                                                                                                                                                        >  1:260831/1‑236 (MQ=255)
        aTCGTCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCa                                                                                                                                                                                                                                  >  1:905250/1‑236 (MQ=255)
            tCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACtt                                                                                                                                                                                                                               <  1:1103028/235‑1 (MQ=255)
                       cTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGGGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  1:996698/236‑1 (MQ=255)
                       cTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  1:1087077/236‑1 (MQ=255)
                        tgtgGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  1:537708/235‑1 (MQ=255)
                          tggtgCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCcaacaa                                                                                                                                                                                                                                         <  2:192554/211‑1 (MQ=255)
                          tggtgCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  2:1235536/233‑1 (MQ=255)
                             tgCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCAGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  2:768818/230‑1 (MQ=255)
                              gCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTg                                                                                                                                                                                                             >  1:776289/1‑235 (MQ=255)
                                  aGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTc                                                                                                                                                                                                          >  2:1191199/1‑234 (MQ=255)
                                     aCCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATgg                                                                                                                                                                                                      >  1:699820/1‑235 (MQ=255)
                                         aGGCGGCGGCGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGACGTGCTGGAAGGTGAGGCCGTGCTTTGCCAAGTGGCGTCGGCCGCTTAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCTACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTtgt                                                                                                                                                                                                                             <  2:1210349/208‑1 (MQ=255)
                                          ggcggcGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACt                                                                                                                                                                                                              >  1:175342/1‑222 (MQ=255)
                                          ggcggcGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCg                                                                                                                                                                                                 >  1:546751/1‑235 (MQ=255)
                                            cggcggTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCg                                                                                                                                                                                                 <  2:450394/233‑1 (MQ=255)
                                                 gTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCa                                                                                                                                                                                        >  1:845746/1‑236 (MQ=255)
                                                 gTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCa                                                                                                                                                                                        <  1:28899/236‑1 (MQ=255)
                                                  tGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGc                                                                                                                                                                                          <  1:1116285/234‑1 (MQ=255)
                                                            tgcgAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGCAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGCGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGgcagca                                                                                                                                                                                                                                  >  1:565159/1‑184 (MQ=255)
                                                             gcgAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGt                                                                                                                                                                            >  1:1038277/1‑236 (MQ=255)
                                                              cgAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGt                                                                                                                                                                            >  1:642567/1‑235 (MQ=255)
                                                                aGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAgg                                                                                                                                                                                                   <  1:94166/211‑1 (MQ=255)
                                                                      gCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGcc                                                                                                                                                                   <  1:877627/236‑1 (MQ=255)
                                                                        tGCGGCCGGTTGCGCCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGa                                                                                                                                                                 <  1:260745/236‑1 (MQ=255)
                                                                              cGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   <  2:722322/181‑1 (MQ=255)
                                                                              cGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTACCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGgc                                                                                                                                                                                                                   >  1:722325/1‑181 (MQ=255)
                                                                               ggTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGaa                                                                                                                                                                                                                                     >  2:66008/1‑162 (MQ=255)
                                                                                           gAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAAc                                                                                                                                                <  1:1142132/234‑1 (MQ=255)
                                                                                            aGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACg                                                                                                                                               <  2:1038274/234‑1 (MQ=255)
                                                                                                 cTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTAt                                                                                                                                           >  1:446578/1‑234 (MQ=255)
                                                                                                 cTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGGCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGGGCCGggc                                                                                                                                                                          >  1:259010/1‑202 (MQ=255)
                                                                                                      aGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGt                                                                                                                                   >  1:311/1‑236 (MQ=255)
                                                                                                        gTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCa                                                                                                                                     <  2:812492/232‑1 (MQ=255)
                                                                                                           aGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATg                                                                                                                              >  1:440530/1‑236 (MQ=255)
                                                                                                              ccGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGt                                                                                                                                                                                 >  1:1565/1‑182 (MQ=255)
                                                                                                                 ggCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGCAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGGTGTCCGgcgcgc                                                                                                                                                                                                               >  1:220409/1‑148 (MQ=255)
                                                                                                                    gttGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAAc                                                                                                                                                <  1:847739/208‑1 (MQ=255)
                                                                                                                        ccAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACg                                                                                                                 <  1:263219/236‑1 (MQ=255)
                                                                                                                          aaGTGGCGTCGGCCGCTGAAGTGGATGCAGTTATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACg                                                                                                                                               <  1:345179/204‑1 (MQ=255)
                                                                                                                              ggCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCg                                                                                                                                                              <  1:956479/185‑1 (MQ=255)
                                                                                                                               gCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGt                                                                                                                                                                                                            <  2:1135781/139‑1 (MQ=255)
                                                                                                                                  tCGGCCGCTGAGGTGGAGGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCGGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCgg                                                                                                                                                                                                                    <  2:970245/128‑1 (MQ=255)
                                                                                                                                  tCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCt                                                                                                       <  1:66008/236‑1 (MQ=255)
                                                                                                                                   cGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGa                                                                                                                      >  1:123060/1‑220 (MQ=255)
                                                                                                                                       cGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTcc                                                                                                                                      >  1:551389/1‑200 (MQ=255)
                                                                                                                                                    gCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACt                                                                                                                    >  2:798866/1‑205 (MQ=255)
                                                                                                                                                     cAGGGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGcc                                                                                                               <  2:261853/209‑1 (MQ=255)
                                                                                                                                                              cAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCa                                                                                                                           <  2:311/188‑1 (MQ=255)
                                                                                                                                                               aGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGccc                                                                                                              >  2:836731/1‑200 (MQ=255)
                                                                                                                                                                   aTTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCaaa                                                                                                                                                                                      <  2:858524/124‑1 (MQ=255)
                                                                                                                                                                   aTTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCaaa                                                                                                                                                                                      >  1:858527/1‑124 (MQ=255)
                                                                                                                                                                    ttGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGcc                                                                                                               <  2:699817/194‑1 (MQ=255)
                                                                                                                                                                       tCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACcgcg                                                                  >  1:1120222/1‑236 (MQ=255)
                                                                                                                                                                                ccGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGc                                                                                                 <  2:123060/196‑1 (MQ=255)
                                                                                                                                                                                   gCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGCCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGAc                                                                                                                                                            >  2:52079/1‑134 (MQ=255)
                                                                                                                                                                                           cccACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGcc                                                       >  1:1096779/1‑227 (MQ=255)
                                                                                                                                                                                            ccACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCt                                                                                                       <  2:381565/178‑1 (MQ=255)
                                                                                                                                                                                               cgggggcggcAGTCGAGGCAGGCAAGAAGGTGTTGGTGGCCAACAAGGAAGCACTAGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATc                                                                                                                                          <  2:95327/136‑1 (MQ=255)
                                                                                                                                                                                                gCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGc                                         >  1:392327/1‑236 (MQ=255)
                                                                                                                                                                                                             gAGGCAGGCAAGAAGGTGTTGCTGGCCGACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGgcg                              >  1:640262/1‑234 (MQ=255)
                                                                                                                                                                                                                       agaagGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGt                   >  1:64111/1‑235 (MQ=255)
                                                                                                                                                                                                                        gaagGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGa                                                                                                     >  2:812086/1‑152 (MQ=255)
                                                                                                                                                                                                                        gaagGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGa                                                                                                     <  1:812089/152‑1 (MQ=255)
                                                                                                                                                                                                                        gaagGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGa                 <  1:1084873/236‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAAc                                                                                                                                                <  2:1196072/100‑1 (MQ=255)
                                                                                                                                                                                                                                 gCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAAc                                                                                                                                                >  1:1196075/1‑100 (MQ=255)
                                                                                                                                                                                                                                 gCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGAGCGGCAAAAGGGTGCCGCGCTGCTGCCGATCGACAGCGAACACAAc                                                                                                                                                >  2:175293/1‑100 (MQ=255)
                                                                                                                                                                                                                                   tGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGGGCCGGCGgt                                                                                          >  2:115204/1‑152 (MQ=255)
                                                                                                                                                                                                                                    ggCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCaa      <  1:1191202/235‑1 (MQ=255)
                                                                                                                                                                                                                                       aaacaaGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCAAGCCt  <  2:632384/235‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCAAGCCt  >  1:707745/1‑235 (MQ=255)
                                                                                                                                                                                                                                        aacaaGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCAAGCCt  <  2:845743/235‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           
CGTGCGTCATCGTCCGGCCTTCGCTGTGGTGCCGAGCACCGAGGCGGCGGTGAGGTTGCGTGCGAGTCTGGCTGCGGCCGGTTGCACCACCGAGGTGCTGGAAGGTGAGGCCGGGCTTTGCCAAGTGGCGTCGGCCGCTGAAGTGGATGCAGTGATGGCAGCCATTGTCGGTGCTGCCGGCCTGCGCCCCACGCTGGCGGCAGTCGAGGCAGGCAAGAAGGTGTTGCTGGCCAACAAGGAAGCACTTGTGATGTCCGGCGCACTGTTCATGGAGGCGGTGCGGCAAAAGGGTGCCGTGCTGCTGCCGATCGACAGCGAACACAACGCTATCTTCCAGTGCATGCCAGGCGACTACGCCCGTGGCCTGAGTGCCGTCGGTGTGCGCCGGATCCTGCTTACCGCGTCCGGCGGGCCGTTCCGCGAAACGCCGGTCGAGGCGCTGCTGGACGTGACTCCGGAGCAAGCCT  >  AE015451/1790241‑1790707

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