Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,421,183 A→C 25.2% T402P (ACC→CCC)  ulaA → L‑ascorbate‑specific enzyme IIC permease component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,421,1830AC25.2% 172.4 / 9.9 92T402P (ACC→CCC) ulaAL‑ascorbate‑specific enzyme IIC permease component of PTS
Reads supporting (aligned to +/- strand):  ref base A (43/19);  new base C (5/21);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTG  >  NC_000913/4420949‑4421420
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
aaCGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTTCGGCGCCGAAGATTTGTCTGGTGGTGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGCGGAACTCACCgg                                                                                                                                                                                                                                            >  1:786067/1‑239 (MQ=255)
aaCGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGGTTTGTCTGGTGATGGGGATGCTCGAAATCTTTGGTTGCGGCTGGGCGGGGAAACTCCCCgg                                                                                                                                                                                                                                            >  1:854131/1‑239 (MQ=255)
   gacGTGGTCTGGGGCGTGATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTAt                                                                                                                                                                                                                                         <  1:456196/237‑1 (MQ=255)
   cgcggggTCTGGGAATTGATGTGGGGCACCGTCGGTCAGCTGGTCGGTGTTGTCCGCCGGGTCGCCTGCGGCGCCTCGATCCTGATTATTCCTGGCTTTAGCCCGAGGTTCTTCTCTAACGCCACTATCGGCGTGTTCGCTAACCGCTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAGTCTTTGGTGGCGTCTTGGCGGTGAAACTCACCGGTAt                                                                                                                                                                                                                                         <  2:1024601/234‑1 (MQ=255)
    gacatgTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTAt                                                                                                                                                                                                                                         <  1:575391/233‑1 (MQ=255)
     cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATg                                                                                                                                                                                                                                       <  1:175928/238‑1 (MQ=255)
     cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                      >  1:565961/1‑239 (MQ=255)
     cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGa                                                                                                                                                                                                                                      <  2:281493/239‑1 (MQ=255)
     cGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCGGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATg                                                                                                                                                                                                                                       <  2:877060/238‑1 (MQ=255)
             ggggCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTgg                                                                                                                                                                                                                              <  1:1121134/239‑1 (MQ=255)
               ggCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGa                                                                                                                                                                                                                             <  2:336391/238‑1 (MQ=255)
               ggCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGGGAAACTCCCCGGTATGAGTGCCTGGAt                                                                                                                                                                                                                            >  1:969426/1‑239 (MQ=255)
                   ttATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGGAGGTTTGTCTGGGGATGGGGATGGTCGAAATCTTTGGGTGCGGCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGGTGGGc                                                                                                                                                                                                                        >  2:592483/1‑239 (MQ=255)
                    tATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGGTTTGTCTGGTGATGGGGATGATCGGAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGCATGGGTGCCTGGATtggga                                                                                                                                                                                                                       >  1:1208084/1‑237 (MQ=255)
                    tATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCa                                                                                                                                                                                                                       >  1:223670/1‑239 (MQ=255)
                      tgtgGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATg                                                                                                                                                                                                                     <  2:565961/239‑1 (MQ=255)
                          gggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGGCTGGTGATGGGGATGATCGGAATCTTTGGTTGCGGCTGGGCGGGGAACCTCCCCGGGATGAGGGCCTGGTTGGGGCTGGCgg                                                                                                                                                                                                                 >  2:945131/1‑239 (MQ=255)
                           ggCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCgg                                                                                                                                                                                                                 >  1:185964/1‑238 (MQ=255)
                              aCCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTTACCGGTATGAGTGCCTGGATGGGCATGGCGGACTg                                                                                                                                                                                                             >  1:201839/1‑239 (MQ=255)
                                cATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCCCTTCGGCGGCTGGCGTGCGGGGCTGGAGCTTTGTCTGGGGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCTTGGCGGACTggc                                                                                                                                                                                                           >  2:78722/1‑238 (MQ=255)
                                     gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGGGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCCCCGGTCAGGGTGCCTCGATGGGCATGGCGGGCTGGTCGaa                                                                                                                                                                                                       >  1:1128998/1‑237 (MQ=255)
                                     gTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGAt                                                                                                                                                                                                       >  1:84044/1‑238 (MQ=255)
                                      tCAGCCGTTTGCGGTTGGCGTCCTGGTCGCCTGCGGCCCCTCGGTGCTGATTATTCCTGGCTTTCGCCCGATGTTCTTCTCTAACGCCCCCATCGGTGTGTTCGCTAACCTCTGCGGCGGCTGGCGTGCGGCGCTGAAGAGTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGGGTCTGGGCTGTGAAACTCGCCGGTTTGAGGGCCTGGATGGGCCTGGCGGACTGGTCGAtt                                                                                                                                                                                                      <  2:991541/238‑1 (MQ=255)
                                          cTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTTGc                                                                                                                                                                                                 >  2:1064861/1‑239 (MQ=255)
                                                gCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGGTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGGCTGGTCGATTCTGGCAccgcc                                                                                                                                                                                           >  1:1224924/1‑239 (MQ=255)
                                                    ttGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatg                                                                                                                                                                                       >  1:1045588/1‑239 (MQ=255)
                                                       gCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCgatgat                                                                                                                                                                                     >  1:103292/1‑238 (MQ=255)
                                                              ggTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGATGGGGGTGATCGAAATCTTTGGCTGCGTTTGGGCGGTGAAACTCACCGGTATGAGAGCCTGGATGGGCATGGCGGGCTGGTCGATTCTGGCACCGCCGATGATGCAAGGct                                                                                                                                                                             >  1:660985/1‑239 (MQ=255)
                                                                ccGCCTCGGGCTCCTCTTCCTTGATTTTCTCTGGTTTTACCCCGGTTTTTTTCTCTAGCGCCCCCTCGGGCGTGTTCGCTAACCATGTCGGGGGGGGGGGGGCGGCGCTGAAGTTTTGTCGGGTGAGGGGGATGATCGAATTTTTTGGTTGCGTCGGGGCGGTAAACCCCCCGGGTATGAGTGCCGGGATGGGCTTGGCGGTCTGGGCGTTTCTGGCACCGCCGATGATGCAAGGcttc                                                                                                                                                                           <  2:1099605/238‑1 (MQ=255)
                                                                    cTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGCAACTCCCCGGTATGAGTGCCTGGATGGGCCTGGCGGACCGGGCGGTTCCGGCAACGCCGGTGATGCAAGGctgtctt                                                                                                                                                                        >  2:69546/1‑233 (MQ=255)
                                                                     tGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGGTGATGCAAGGCTtcttct                                                                                                                                                                       >  2:382131/1‑238 (MQ=255)
                                                                        gcctcctcGTCCCTGATTGTTCCTGTTTTTACCCCGATGTTCTCCCCTACCGCCACCACGGGTGTGTGCGCTACCCATTGCGGCGGGTGGGGTGGGGGGCTGAAGTTTTGTCTGGTGATGGGGATGATGAAAATCTTTGGTGGCGTTGGGGGGTGGAAACTCCCCGTTATGAGTGCCTGGATGGGCTTGGGGGCCTGGCCGATTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCa                                                                                                                                                                    <  2:1208428/236‑1 (MQ=255)
                                                                              tcGACCTTGATTATCCCTGTCTTTACCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTACCCACTGCGGCGGCGGGCGTGCGGCGCTGAAGATTTGTCGGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTAt                                                                                                                                                             <  1:964938/239‑1 (MQ=255)
                                                                                  tCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCg                                                                                                                                                           >  1:701608/1‑237 (MQ=255)
                                                                                    cTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACTCCACCATCGGCGTGTTCGCTCACCTCTTCGGCGGCTGGCGTGCCGCGCTGAATATTTTTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCATCGGTATGAGGGCCTTGATGTGCATGGCGGACTTGTCTGTTCTGGCGCCGCCGATGATGCAATGCATCTTTGCCTACGGTCTTGCCtt                                                                                                                                                       >  2:21678/1‑239 (MQ=255)
                                                                                       attattCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTAt                                                                                                                                                    >  1:926517/1‑239 (MQ=255)
                                                                                            tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGGCTGGTCGATTCTGGCACCGCCGATGATGCCAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCg                                                                                                                                               >  1:882086/1‑239 (MQ=255)
                                                                                            tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCg                                                                                                                                               >  1:190042/1‑239 (MQ=255)
                                                                                            tCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCg                                                                                                                                               >  1:926662/1‑239 (MQ=255)
                                                                                                   tttATCCCGATTTTTTTCTCTAACGCCACCATGGGCGTTTTCGCTACCCCCTTGGGCGGCTGGCGTGGGGCGCTGAAGATTTGTCGGGTGATGGGGATGATCAAAATCTTTGGTTGCGTCGGGGCGGTAAACCTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTcatcat                                                                                                                                        <  1:376798/239‑1 (MQ=255)
                                                                                                     tATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTg                                                                                                                                      >  2:857928/1‑239 (MQ=255)
                                                                                                      aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGGCTGGTCGAATCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGt                                                                                                                                     >  1:821381/1‑239 (MQ=255)
                                                                                                      aTCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATTATTGt                                                                                                                                     >  1:913474/1‑239 (MQ=255)
                                                                                                        cccGTTGTTCTTGCCTACCGCCCCCATCGGGTTGTCCGCAAACCCCTGGGGGGGCGGGGGGGGGGCGCGAAAGTTGTGTCGGGTGTGGGGGATGATGGAATTTTTTGGTTGCGTCTGGGCGGTGAAACTCCCCGTTATGAGTGCGTGGATGGGCATGGCGGCCTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTaa                                                                                                                                   <  2:292646/239‑1 (MQ=255)
                                                                                                             tGTTCTTCTCTACCGCCGCCTGGGGCGGGTGCTCTAACCTCTTGGGGGGTTGGCGTGGGGCGCTGAAGATTTGGTGGGTGAGGGGGATGATCGAAATCTTGGGTGGCGTCGGGGCGGTGAAACTCCCGGGTTTGGGTGCCTGGATGGGCATGGCGGACTGGTCTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                              <  2:628241/239‑1 (MQ=255)
                                                                                                             tGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGc                                                                                                                               >  1:90252/1‑238 (MQ=255)
                                                                                                             tGTTCTTCTCACGCGCCCCCATCGGCTTTTTCGCAAACCCCTGCGGCGGCTGGCGGGCGGCGCGGAAGTTTTTGTTGGTGATGGGGAGGATAGAATTCTTTGGTTGCGTCGGGGCGGTGAAACTCCCCGGTATGAGTGCCTGGATGGGCTTGGCGGCCTGGTCTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                              <  2:903201/239‑1 (MQ=255)
                                                                                                              tttctttCCTAACGCCACGGGCGGCTGGTTCGCTACCCTCTTGGGCGGCGGGCGGGGGGCGCGGAAGATTGGTTGGGGGGTGGGGAGGATGAAATTCTTGGTTGGCGTTGGGGGGGTGAAACTCCCCGATATGCGTTCCTGATGGGGCATGGCGGACGGGTCGTTTCGGGCGCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCAc                                                                                                                             <  2:855138/238‑1 (MQ=255)
                                                                                                              tttcttcTCTACGCCCCCATGCGGCTGTTTCGCTAACCTGTGCGGCGGTGGGCGGGGGGCGGTGAAGTTTTGTCGGGTGATGGGGATGATGAAATTTTTGGTTTGCGTCGGGGGGGGGAACCTCCCGGGTAGGAGTGCGTGGAGGGGCTTGGGGGCTTGGTCTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCAc                                                                                                                             <  2:1121235/238‑1 (MQ=255)
                                                                                                              gTTCCTCTATAACTCCACCGGCGTGTTTTTCGGTAACCCCTTGGGGGGCGGGGGGGCGGCGCTGATGTTTTGGTGGGTGATGGGGTTGATGGAAATCTTGGGTTGGGTTTGGGGGGGAAACCGCACCGGTATGAGTCCCTGGAGGGGCGTGGCGGTCTGTTCGTTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCa                                                                                                                              <  2:572825/238‑1 (MQ=255)
                                                                                                                 ctcctagacccccaccaTCGGTGTTTGCGGTACCCATTTCGGCGGCGGGCGTGCGAGGATGAAGTTTGGTCGGGGGTTGGTGAGACTAGAATTTTTGGGGTGCGGGTGGGCGGTGAACCTCCCGGTTGTGGGGCCCGGGATGGGCGTGGCGGCCTGGTCGATTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTgg                                                                                                                          <  2:660985/228‑1 (MQ=255)
                                                                                                                  ttctCTACCGCCCCCATCGGCGTGTTCGCTAACCGCTGGGGGGGCTGGCGGGGGGCGCGGAAGATTTGTCTGGTGATGGGGATGATCAAATTTTTTGGTGGGGTTGGGGCGGTGAACCTCCCCGGTATGGGTGCCGGGATGGGCATGGGGGACTGGTCGTTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGc                                                                                                                         <  2:1056797/239‑1 (MQ=255)
                                                                                                                  ttctCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTg                                                                                                                           >  1:602350/1‑237 (MQ=255)
                                                                                                                   tctcTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGGAATCTTTGGGTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTTGCACCGCCGAGGATGCAAGGCTTCTTCTCCATCGGTATCGACTTTATGGCCGTCATCCTTGTAATTCAACTGGCt                                                                                                                        >  2:104968/1‑239 (MQ=255)
                                                                                                                      ctcgcGCCCCCAGCGGCGTGGCCTCTACCCGCTGCGGGGTGTGGGGTGGGGCGGTGAATATTTGTTTGGTGATGGGGATGATCGAACTCTTGTGTTGCGTCGGGGCTGTAAATCTCAGCGGTATGAGTGCCTGGAGGGGCATGGCGGACTGGTCGATTCTGGCTCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTtat                                                                                                                     <  2:667824/235‑1 (MQ=255)
                                                                                                                       tAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCCTCGGTATCGCCTTTATGGCCGTCATCCTTGTAATTGCACTGGCTtata                                                                                                                    >  1:1164380/1‑239 (MQ=255)
                                                                                                                       tAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGGAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTtata                                                                                                                    >  2:1078768/1‑239 (MQ=255)
                                                                                                                       tAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTtata                                                                                                                    >  1:896154/1‑239 (MQ=255)
                                                                                                                       tAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGGATCGCCTTTATTGCCGTTCTCCTTGTAATTGCCCTGGCTtata                                                                                                                    >  2:988470/1‑239 (MQ=255)
                                                                                                                         aCGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATg                                                                                                                  >  2:486209/1‑239 (MQ=255)
                                                                                                                                      gCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTgg                                                                                                       >  1:324546/1‑238 (MQ=255)
                                                                                                                                           ttCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGGATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCTCTGGCcgcgc                                                                                                 >  1:520588/1‑239 (MQ=255)
                                                                                                                                                tAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgc                                                                                            >  1:1162191/1‑239 (MQ=255)
                                                                                                                                                 acccccTTGGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGGGAGGGGGATGATCGAATTTTTTGGTTGCGTCTGGGGGGTAAACCTCCCCGGTATGAGTGCCTGGAGGGGCATGGCGGACTGGTCGTTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcg                                                                                           <  1:592483/234‑1 (MQ=255)
                                                                                                                                                  aCCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcgc                                                                                          >  2:844168/1‑239 (MQ=255)
                                                                                                                                                   ccgttgcggcggTTGGCGGGCGGCGCTGAAGGTTTGTCGGGGGGGGGAGATGAGCGAATTCTTGGGTTGCGTGGGGGCGGAGAAGCCCCCCGGTATGATTGCTGGGAGGGGCAGGGCGGATGGGTAGAGTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTgcgcg                                                                                         <  2:734313/233‑1 (MQ=255)
                                                                                                                                                                     tggggcgGTGAAGTTGTGTGTGGTGAGGGGGATGATCAAAATCTTGGGTGGCGTCTGGGCGGTGAAACTCCCCGGTATGAGTCCCGGGGTGGGCATGGCGGACGGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGaa                                                                        <  2:1224924/235‑1 (MQ=255)
                                                                                                                                                                            tggAAGATTGGTCGGGTGGTGGGGGTGAGGGAATTCTTTGGTGGGGTCGGGGCGGAGAATCTCCCCGGTATGGGTGCCGGGAGGGGCATGGCGGACTGGTCGATTCTGGCTCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                                 <  2:105648/236‑1 (MQ=255)
                                                                                                                                                                            aggAAGTTTTGTCGGGTGATGGGGGGGATGAAAATCTTGGGTGGCGTCGGGGGGGTGAACCTCACCGTTATGGGTGCGTGGATGGGCAGGGCGGACTGGTCGATTGTGGCCCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAaacaac                                                                 <  2:1206652/236‑1 (MQ=255)
                                                                                                                                                                                agttttgtttgGTGAGGGGTATGATCAAATTTTTTGGTGGCGTGTGGGGGGTAAAACTCCCGGTTTGGGGTGCCTGGTGGGGCATGGCGGCCTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCa                                                            <  2:90252/236‑1 (MQ=255)
                                                                                                                                                                                     tgtgtgGGGATGGGGATCATCGAATTCTTTGGTTGCGTCTGGGCGGTGAACCTCCCCGGTATGAGTGCCTGGAGGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACa                                                       <  1:975702/239‑1 (MQ=255)
                                                                                                                                                                                       tCTGGTGGGGGGGATGATCGAAATCTTTGGTTGCGTCGGGGCGGTGAACCCCCCCGGTATGAGTGCCTGGAGGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                     <  1:945131/239‑1 (MQ=255)
                                                                                                                                                                                       ccgggggATGGGGAGGATCGAACTTTTGGGCGTCTTCGGGGCGGGCAACCCCCCGGGTAGGGGCGCCTGGCGGGGCTGCGCGGATGGGTCCATTCTGGCTCCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGt                                                     <  2:1164380/233‑1 (MQ=255)
                                                                                                                                                                                        ctgggtggggggATGATAGAAGTCTTGGTGTGGTTGTGGGGGGTAAAACCCTCCGTTATGAGTGCTGGGATGGGCGTGGGGGATTGGTGTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTc                                                    <  2:863531/231‑1 (MQ=255)
                                                                                                                                                                                             gATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGCTTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCCTCGGTATCGCCTTTATGGCCGTCATCCTTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAACTATATGCATTAAAACAACTGGCAGAACAGCCTGCtt                                               >  2:332381/1‑239 (MQ=255)
                                                                                                                                                                                                 gggATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTtaata                                           >  2:327984/1‑239 (MQ=255)
                                                                                                                                                                                                  gtatgatCGAATTTTTTGTTTGCGTCTGGGGGGTGAACCTCCCCGGTATGGGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTtaata                                           <  1:824170/236‑1 (MQ=255)
                                                                                                                                                                                                  ggttgagCGAATTCTTGGGTGGGGGCGGGGGGGTGAAACTCCCCGTTAGGAGTGCCTGGAGGGGCGTGGCGGACGGGTTGATTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTtaata                                           <  2:841579/235‑1 (MQ=255)
                                                                                                                                                                                                    atgatCAAATTCTTTGGTTGCGTGTGGGGGGTGAAACTCACCGGTATGAGGGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAg                                         <  1:460193/238‑1 (MQ=255)
                                                                                                                                                                                                    ataaccagaaTTTTGGGTTGGTTCGGGGGGGTAAAACTCGCCGGTAGGAGTGCCTGGTGGGGCTGGGCGGAGTGGTGTATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAg                                         <  2:821381/230‑1 (MQ=255)
                                                                                                                                                                                                      gatCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGCCTGCTTAATAAGGAg                                      >  1:1193059/1‑239 (MQ=255)
                                                                                                                                                                                                        tggatgttttgggTTGCGGCGGGGGGGTAAAACTCCCCGGTATGAGTGCGTGGATGGGCATGGCGGACTGGGCTATTCTGCCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                    <  2:520588/228‑1 (MQ=255)
                                                                                                                                                                                                        tCGAAGTCTTGGGTGGCGTGGGGGGGGTAAAACTCCCCGGTATGGGTTTCTGGATGGGCGTGGCGGACTGGGCTTTTCGGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCCTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                    <  2:871672/239‑1 (MQ=255)
                                                                                                                                                                                                          gAAATTTTTGGTTGCGTCGGGGGGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGtt                                    <  1:387356/237‑1 (MQ=255)
                                                                                                                                                                                                                  tGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGTGCTGCGCGCAGAAGCAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATg                           >  2:1122929/1‑238 (MQ=255)
                                                                                                                                                                                                                   ggTGGCGTCTGGGCGGTGAAACTCACCGGTATGAGGGCCTGGAGGGGCAGGGGGGACGGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                         <  1:207431/239‑1 (MQ=255)
                                                                                                                                                                                                                     ggtgggcgggggggggTGAATGTCACCGTTAGGAGGGCCAGGAGGGGAATGGCGAACAGGTCAATTCTGGCTCCGCCGATGATGCAGGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGAc                         <  2:1178487/229‑1 (MQ=255)
                                                                                                                                                                                                                         gctgcggggCGGTAAAACCCCCCGGTAGGAGTGCCTGGATGGGCGTGGCGGATTGGTCGTTTCTGCCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgta                     <  1:382131/233‑1 (MQ=255)
                                                                                                                                                                                                                            tatGGGCGGTGAAACTCACAGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCCTCTTCTCCCTCGGTATCGCCTTTATGGCCGTCCTCCTTGTAATTGCACTGGCTTATATGTTCTTCGCTGGGCGCGCGCTCCGCGCAGAACAAGATGCAGCAAAACAACGGGCAGCACAGTCTGCTTAATAAGGAGTTTTGATTATGACcgtacgt                  >  2:689467/3‑239 (MQ=255)
                                                                                                                                                                                                                                    gTGAACCTCACGGGTAGGAGGGCCGGGAGGGGGAGGGCGGACGGGTGGATCCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTgg           <  1:21678/238‑1 (MQ=255)
                                                                                                                                                                                                                                     tGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTg            >  1:164749/1‑236 (MQ=255)
                                                                                                                                                                                                                                             cACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGGTTTGATTATGACCGTACGTATTCTGGCtgtgtgtg  >  1:621212/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
AACGCCGTGGTCTGGGGCTTTATGTGGGGCACCATCGGTCAGCTGATTGCGGTTGGCATCCTGGTCGCCTGCGGCTCCTCGATCCTGATTATTCCTGGCTTTATCCCGATGTTCTTCTCTAACGCCACCATCGGCGTGTTCGCTAACCACTTCGGCGGCTGGCGTGCGGCGCTGAAGATTTGTCTGGTGATGGGGATGATCGAAATCTTTGGTTGCGTCTGGGCGGTGAAACTCACCGGTATGAGTGCCTGGATGGGCATGGCGGACTGGTCGATTCTGGCACCGCCGATGATGCAAGGCTTCTTCTCCATCGGTATCGCCTTTATGGCCGTCATCATTGTAATTGCACTGGCTTATATGTTCTTCGCTGGCCGCGCGCTGCGCGCAGAAGAAGATGCAGAAAAACAACTGGCAGAACAGTCTGCTTAATAAGGAGTTTTGATTATGACCGTACGTATTCTGGCTGTGTGTG  >  NC_000913/4420949‑4421420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: