Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,168 G→A 6.0% intergenic (+287/‑291) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,1680GA 6.0% 231.7 / 6.0 84intergenic (+287/‑291)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base A (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGA  >  NC_000913/374937‑375399
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
ggggcgTGACCGCCTTATCCGGCCTACGGAAGGAGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGcgcg                                                                                                                                                                                                                                   <  2:191032/237‑1 (MQ=255)
   gcgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGGCTGGTAAGACGCGTTAGCGTCGCCTCCGGGCTCTGCGCACGACTGCCGGATGCGGGGTGAACGCCCTGTTCGGCCTCCGGGCGGCGCggga                                                                                                                                                                                                                                >  2:258449/1‑239 (MQ=11)
    cgTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGGGAGGCGCGGGAATTTGTAGGCCTGGTAAGACGCGTTAGCGCTGCATCAGGCATCTGCGCACGACTGCCCGATGTGGCGTGGACTGCTTAGACCGCCTCCGGGGGGCGCgaggg                                                                                                                                                                                                                               >  2:936343/1‑237 (MQ=11)
                aTCCGGCCTACGGATGGCGTGCGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATCAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGGGGCCGCCTTATCCGGCCTCAGGATGGCGCGGGCATTTGTAGGGCTg                                                                                                                                                                                                                   >  2:564365/1‑239 (MQ=11)
                      ccTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGCGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAg                                                                                                                                                                                                              >  1:1090769/1‑238 (MQ=255)
                       cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                                                             >  1:933603/1‑238 (MQ=255)
                       agaCGGATGGCGGGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                                                                                                                             <  2:216093/236‑1 (MQ=255)
                           ggAGGGCGTGAGAATTTGTAGGTGTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                                                         <  2:769161/238‑1 (MQ=12)
                            gATGGCGTGAGAAATTGTAGGTCTGAGAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGt                                                                                                                                                                                                        <  1:223530/238‑1 (MQ=255)
                               ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAg                                                                                                                                                                                                                        <  2:723115/219‑1 (MQ=255)
                               ggCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAg                                                                                                                                                                                                                        >  1:723115/1‑219 (MQ=255)
                                 cGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                                                                                                                                                   <  1:258449/238‑1 (MQ=40)
                                  gTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGt                                                                                                                                                                                                  <  1:1016169/238‑1 (MQ=40)
                                    gagaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTCTTAAGCACCGTTAGCGTccg                                                                                                                                                                                               >  1:2587/1‑237 (MQ=255)
                                     agaATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa                                                                                                                                                                                              >  1:820452/1‑239 (MQ=39)
                                        aTTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                                                                                                                           <  1:522109/239‑1 (MQ=38)
                                           tGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGGATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCCGGCATCTGCGCACGACCGCCCGATGCGGCGTGAACGCCGTCTCCCGCCCACTGGGGCCGCGGGGAATTGTATACCTGATATGACGCGCTAGCGTCGCCTCCGGc                                                                                                                                                                                        >  2:962146/1‑239 (MQ=11)
                                              aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                                                                                      <  2:1133983/238‑1 (MQ=38)
                                              aGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                                                                                                                                                      <  2:943677/238‑1 (MQ=38)
                                                          aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGGGGCGCGAGGATTTGTAGGCCTGACAAGACGCCTTAGCGTCGCCTCAGGCATCTGCTCACGGC‑t                                                                                                                                                                         >  2:190559/1‑238 (MQ=255)
                                                               ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGCCGAATCAGGCATCTGCGCACGACTGCCGGATGCGGCGCGGACGCCTTATCCGGCCTTCGGGTGGCGGGGGAGCTTGTAGGCCAGATAAGCCGCGTTAGCGTCGCAATAGGCATCTGCGCACGAC‑AGCCGGa                                                                                                                                                                   >  2:1157195/1‑239 (MQ=11)
                                                                    gTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATgcg                                                                                                                                                               <  1:403020/238‑1 (MQ=38)
                                                                        cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTa                                                                                                                                                          >  1:387553/1‑239 (MQ=40)
                                                                         gtgaggcgttcgcgcACGACTGGCGGATGCGGCGGAAACGGCTTGTCCGGCCTACGGATGGCGCGGGAATTTGTCGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTaa                                                                                                                                                         <  2:99686/228‑1 (MQ=14)
                                                                             ggCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACg                                                                                                                                                      >  1:141420/1‑238 (MQ=255)
                                                                                 tcccagcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTa                                                                                                                                                 >  1:957138/6‑239 (MQ=25)
                                                                                   tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                                               <  2:94553/239‑1 (MQ=40)
                                                                                     cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGCGCGGGAATTTGTAGGCCTGATAAGACGCGGTAGCGGCGCCTCAGGCATCTGCGCCCGAC‑TGCCGGAGGCGGCGCAAACCCCTTATccc                                                                                                                                             >  2:593937/1‑238 (MQ=11)
                                                                                        aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACGGCTGGCGCGGGAATTTGTCAGCCCGTTTAGATCCGCTATCTTCCCATCACGCATCCGTGCACCAC‑TGCCGGCCGCGCCGTGACCGCCTTTTCCgcc                                                                                                                                           >  2:268250/1‑236 (MQ=11)
                                                                                                         cgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAAAGCCTTCTCCCGCCTACCGGGGGCGCGAGAAGTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCCGGCCTCTCCCGACGACATGCGGGAGGC‑GCGAAACCGGCTTATCCGGCCTACGGCGTGGGCGGGaa                                                                                                                         >  2:244938/1‑239 (MQ=11)
                                                                                                         agtgaccggagtATCCGGCTGACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                          <  2:1068134/227‑1 (MQ=255)
                                                                                                                  cTTATCCGGCGTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGc                                                                                                                <  2:957138/239‑1 (MQ=255)
                                                                                                                  cTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTCGCCTGTGCATTTGTAgc                                                                                                                 >  2:468245/1‑237 (MQ=255)
                                                                                                                   ttATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAgg                                                                                                                 >  1:121491/1‑237 (MQ=255)
                                                                                                                   ttATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGGCGCATCAGGCATCTGCGCACGAC‑TGCCGGGTGCGCCGTAAACGCCTTATCCGGCCTACGGCTGGCGGGGGGATTCGTTGGc                                                                                                                >  2:1056094/1‑238 (MQ=14)
                                                                                                                        cGGCCTGCGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGCCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                                                          <  1:1107475/239‑1 (MQ=255)
                                                                                                                          gCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAg                                                                                                                                                                                                                        <  1:411480/128‑1 (MQ=255)
                                                                                                                           ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGA‑AATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAc                                                                                                      >  2:35823/1‑239 (MQ=255)
                                                                                                                            cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                       <  1:118298/238‑1 (MQ=25)
                                                                                                                            cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGa                                                                                                       <  1:576383/238‑1 (MQ=25)
                                                                                                                                gggtGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGt                                                                                                  >  1:691181/4‑239 (MQ=255)
                                                                                                                                 ggtGGAGCGAGAATGTGTAAGCCTGCTGAGAAGCGTTAGCGTCGCATCAGGCGTCTGCGCACGGCTGCCGGATGCAGGGTGAACGCCTTATCCGGGCTTCGGGTGGCGCGAGAAGTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGt                                                                                                  <  2:714769/236‑1 (MQ=11)
                                                                                                                                  gtGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                                                                 <  2:691181/237‑1 (MQ=255)
                                                                                                                                  atggggggagAATTTGTAGACCTGCGTAGACGCGTTAGCGGCGGATCAGGCGGCTGCGCACGACTGCAGGATGCGGCGTGAACGGCTTATCCTGCCTACGGGGGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                                                                 <  2:359401/229‑1 (MQ=11)
                                                                                                                                    ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAg                                                                                                                                                                                                              >  2:411480/1‑128 (MQ=40)
                                                                                                                                    ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TTCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                                              >  2:1023389/1‑239 (MQ=255)
                                                                                                                                                          aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGCATTTGTAGGCCTGATCAGACGCGTTCGCCTCGCATCCGCCCTCTGccccc                                                                        >  2:474884/1‑235 (MQ=11)
                                                                                                                                                          aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGAAAACCCCTTATCCGGCCTACGGATGGCGGGGGAATTTGTAGGCTTGATAAGACTCGTTAGCGTCGCATCATGCATCTGcaaa                                                                         >  2:539354/1‑236 (MQ=255)
                                                                                                                                                            aaGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGa                                                                      <  1:936343/239‑1 (MQ=38)
                                                                                                                                                               aCGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTCGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAACAGCCTTATCCGGCCTACGATTCGCGCGGGAATTGGCAGGGCTGATAAGACGCCTCAGCGTCCGATCAAGCATCCGCACACGAcc                                                                    >  2:340857/1‑237 (MQ=11)
                                                                                                                                                                   gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGGCTGATAAAACGCGGTAGCGTCGCATCAGGCATCTGCGCCCGAC‑TGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGCTGGCGCGGTCATTCGTAGGCCTGATACGTCGCGTAACCGCGGCACCCGGCATCAGCACCCGCCTGCCgg                                                               >  2:772453/1‑239 (MQ=11)
                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCCCACGACTGCCGGa                                                              >  1:924419/1‑239 (MQ=37)
                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGa                                                              >  1:927406/1‑239 (MQ=38)
                                                                                                                                                                    ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTAACCGGCCTTCGGGTGGGGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGTGCAAGAC‑TGCCGGATGCGGCGTAAACGTCTTATCCGGACTCCGGAAAGCGCGGGAAATGTTACGCCTGATTAGTCGCGTTAACGTCACATCAGGCTACTGCACACGACTCACGGa                                                              >  2:921671/1‑239 (MQ=11)
                                                                                                                                                                           cGCATCAGGCATCTGCGCCCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                       <  1:921671/239‑1 (MQ=255)
                                                                                                                                                                            gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                      <  1:1062052/239‑1 (MQ=40)
                                                                                                                                                                            gCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                      <  1:474884/239‑1 (MQ=40)
                                                                                                                                                                             cATCGGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGTCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                      <  1:1157195/238‑1 (MQ=255)
                                                                                                                                                                             cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                      >  1:731114/1‑238 (MQ=40)
                                                                                                                                                                             cATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATaa                                                      >  1:1110573/1‑238 (MQ=40)
                                                                                                                                                                                     aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGt                                             <  1:910391/239‑1 (MQ=255)
                                                                                                                                                                                     aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGA‑AATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                            <  1:35823/239‑1 (MQ=255)
                                                                                                                                                                                      tCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTGCGGGTGGCGCGAGAAGTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                            <  2:673806/239‑1 (MQ=255)
                                                                                                                                                                                        tGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAACCACGCCTCAGGCCTCCGCCCACAACTGCCGGCTGCGCTACACGTCTTGccc                                           >  2:902308/1‑238 (MQ=14)
                                                                                                                                                                                          cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCCTCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATCAGACGCGTTGGCGTCGCATGAGGCATCTGCACCCCACTTCCCGATGCCATACAAGTATTGTCAgg                                         >  2:159003/1‑238 (MQ=14)
                                                                                                                                                                                          cgcACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCgg                                         >  1:738189/1‑238 (MQ=255)
                                                                                                                                                                                                 cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGTCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAtt                                 <  1:466282/239‑1 (MQ=255)
                                                                                                                                                                                                 cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAtt                                 <  1:772453/239‑1 (MQ=255)
                                                                                                                                                                                                    ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                          >  2:487926/1‑147 (MQ=39)
                                                                                                                                                                                                    ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                                                          <  1:487926/147‑1 (MQ=39)
                                                                                                                                                                                                         tGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATTAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACCTTTCACCCgc                         >  2:130693/1‑238 (MQ=255)
                                                                                                                                                                                                              cggaaaCGCCTTATCCGGACTACGGGTAGCGCGAGAATTTGTAGGCCTGATGAGAGGCGTTAGCGTCGCATCAGGCAGCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACGCGACTGCCGGATGCGATAAACGTCTTGTCCGGGCGACATTTCGCTCGTAGGc                     <  2:524731/234‑1 (MQ=17)
                                                                                                                                                                                                              cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc                                                                                                                                                                                                                      <  2:526115/46‑1 (MQ=255)
                                                                                                                                                                                                              cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc                                                                                                                                                                                                                      >  1:526115/1‑46 (MQ=255)
                                                                                                                                                                                                              cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                     <  1:962146/238‑1 (MQ=255)
                                                                                                                                                                                                              cgTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                     <  2:927406/238‑1 (MQ=255)
                                                                                                                                                                                                                ggatcGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                   <  2:1090769/234‑1 (MQ=255)
                                                                                                                                                                                                                  aaCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                >  2:268214/1‑239 (MQ=255)
                                                                                                                                                                                                                   gcGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                <  1:902308/237‑1 (MQ=255)
                                                                                                                                                                                                                   aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                <  1:121999/238‑1 (MQ=255)
                                                                                                                                                                                                                   aCGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                <  2:280934/238‑1 (MQ=255)
                                                                                                                                                                                                                        ttATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAAt          >  1:20250/1‑239 (MQ=255)
                                                                                                                                                                                                                           tCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGt       >  1:782376/1‑239 (MQ=255)
                                                                                                                                                                                                                                ccTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTACGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCGGGATGCGCTAAACGTCTTGTCCGGCCTCCATTTCGCCCGTAGGCAGTCCTTACATAGTTCTGa  >  2:311985/1‑239 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                        
GCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGAC‑TGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGA  >  NC_000913/374937‑375399

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