Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,426,929 C→T C165Y (TGC→TAC)  yjfZ ← DUF2686 domain‑containing protein YjfZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,9290CT97.8% 305.8 / ‑2.0 89C165Y (TGC→TAC) yjfZDUF2686 domain‑containing protein YjfZ
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (39/48);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGAATTA  >  NC_000913/4426689‑4427174
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
ggCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCAt                                                                                                                                                                                                                                             >  1:384500/1‑251 (MQ=255)
      aGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAgg                                                                                                                                                                                                                                                >  1:929623/1‑242 (MQ=255)
      aGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAgg                                                                                                                                                                                                                                                <  2:929623/242‑1 (MQ=255)
                 aGCGCAACTGATGTCTCGCTCCTTCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGATGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTATTTCAGATGACCAGTAATTAGTCAGTTGGGCAGGAAGGGTTa                                                                                                                                                                                                                            <  2:183461/251‑1 (MQ=255)
                                 cGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGt                                                                                                                                                                                                            >  1:614248/1‑251 (MQ=255)
                                 cGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCAGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCg                                                                                                                                                                                                             <  2:404469/250‑1 (MQ=255)
                                    tCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTg                                                                                                                                                                                                         >  2:633393/1‑251 (MQ=255)
                                        gCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCg                                                                                                                                                                                                     <  2:289705/251‑1 (MQ=255)
                                            cACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGtat                                                                                                                                                                                                 <  2:924835/251‑1 (MQ=255)
                                             aCTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGtat                                                                                                                                                                                                 >  1:751592/1‑250 (MQ=255)
                                                    cgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCCGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCCGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTa                                                                                                                                                                                                                                                  >  2:1146244/1‑194 (MQ=255)
                                                    cgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTa                                                                                                                                                                                                                                                  <  1:1146244/194‑1 (MQ=255)
                                                    cgAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGa                                                                                                                                                                                         <  2:755022/251‑1 (MQ=255)
                                                     gAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTgg                                                                                                                                                                                          >  1:183397/1‑249 (MQ=255)
                                                      aaacccaaacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGCGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGt                                                                                                                                                                                       >  1:118927/1‑251 (MQ=255)
                                                      aaacccaaacccGAGTAGATTGCGGCTGAACGTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGt                                                                                                                                                                                       <  2:751592/251‑1 (MQ=255)
                                                        acccaaacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGtgc                                                                                                                                                                                     >  2:466667/1‑249 (MQ=255)
                                                             aacccGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATa                                                                                                                                                                                 >  2:319080/1‑250 (MQ=255)
                                                                 cGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCa                                                                                                                                                                                                                                      >  1:438254/1‑193 (MQ=255)
                                                                 cGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCa                                                                                                                                                                                                                                      <  2:438254/193‑1 (MQ=255)
                                                                   aGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGt                                                                                                                                                                          <  1:10781/251‑1 (MQ=255)
                                                                    gTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTc                                                                                                                                                                         <  2:328078/251‑1 (MQ=255)
                                                                     tAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATc                                                                                                                                                                                                                   <  1:793830/208‑1 (MQ=255)
                                                                     tAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATc                                                                                                                                                                                                                   >  2:793830/1‑208 (MQ=255)
                                                                        aTTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCaa                                                                                                                                                                     <  1:707610/251‑1 (MQ=255)
                                                                             gggtGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAAcc                                                                                                                                                                 >  2:770010/4‑250 (MQ=255)
                                                                                tGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGCATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACtt                                                                                                                                                             >  2:479358/1‑251 (MQ=255)
                                                                                   aCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                                          >  1:126648/1‑203 (MQ=255)
                                                                                   aCTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGGTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCATATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGttt                                                                                                                                                                                                          <  2:126648/203‑1 (MQ=255)
                                                                                         tAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtatt                                                                                                                                                     <  1:428921/250‑1 (MQ=255)
                                                                                         tAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtatt                                                                                                                                                     >  2:428921/1‑250 (MQ=255)
                                                                                           aGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtattgtg                                                                                                                                                  <  2:614248/251‑1 (MQ=255)
                                                                                             gAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCgtgtattgtgta                                                                                                                                                <  1:500138/251‑1 (MQ=255)
                                                                                                cGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTAttt                                                                                                                                             >  2:252209/1‑251 (MQ=255)
                                                                                                 gCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTAttt                                                                                                                                             >  2:579554/1‑250 (MQ=255)
                                                                                                  cAGACCAATTTTCAAATCATGGGGTAAATCAGAATCCGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGACTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGAAGGAAGGGTTACAACGTATACTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTAc                                                                                                                                           >  1:279194/1‑251 (MQ=255)
                                                                                                   aGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACt                                                                                                                                          >  1:977901/1‑251 (MQ=255)
                                                                                                      ccAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGtt                                                                                                                                       <  2:1242830/251‑1 (MQ=255)
                                                                                                        aaTTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTgg                                                                                                                                     <  1:770010/251‑1 (MQ=255)
                                                                                                           tttCAAATCATGGGGTAAATCAGAATTCGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGt                                                                                                                                   <  1:185730/250‑1 (MQ=255)
                                                                                                             tCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTc                                                                                                                                 <  1:1082266/250‑1 (MQ=255)
                                                                                                              cAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCg                                                                                                                                <  2:583580/250‑1 (MQ=255)
                                                                                                              cAAATCATGGGGTAAATCAGAATAAGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGGTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACCACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCg                                                                                                                                >  1:525339/1‑250 (MQ=255)
                                                                                                               aaaTCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGaa                                                                                                                              <  2:977901/251‑1 (MQ=255)
                                                                                                                  tCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTACTGTTGGGTTCGAAGGAt                                                                                                                          >  2:163956/1‑251 (MQ=255)
                                                                                                                      ggggTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTc                                                                                                                        <  1:319080/250‑1 (MQ=255)
                                                                                                                         gTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTATAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTTCCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCgcca                                                                                                                    <  2:279194/251‑4 (MQ=255)
                                                                                                                           aaaTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATc                                                                                                                  >  1:41809/1‑251 (MQ=255)
                                                                                                                           aaaTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGGCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATc                                                                                                                  <  1:53773/251‑1 (MQ=255)
                                                                                                                            aaTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCt                                                                                                                 <  2:41809/251‑1 (MQ=255)
                                                                                                                             aTCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCtt                                                                                                                >  2:232962/1‑251 (MQ=255)
                                                                                                                                    tACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaata                                                                                                          >  1:951073/1‑250 (MQ=255)
                                                                                                                                      cgcgTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaataaat                                                                                                       >  1:1236643/1‑251 (MQ=255)
                                                                                                                                         gTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTCTCCTATCGTTTGGCCTGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATAtt                                                                                                    <  2:525339/251‑1 (MQ=255)
                                                                                                                                               cAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaata                                                                                                          >  1:1220819/1‑239 (MQ=255)
                                                                                                                                               cAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCaata                                                                                                          <  2:1220819/239‑1 (MQ=255)
                                                                                                                                                 aTGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCt                                                                                            <  1:685244/251‑1 (MQ=255)
                                                                                                                                                                     aGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCcgc                                                                        <  1:633393/251‑1 (MQ=255)
                                                                                                                                                                      gggCATTAGAACCGACTGTAACAGTCCGAAGTTGGTCGTTATAAAGATCGCTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAACACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCATTAAATATTTCAGTTCTGCGGTTAACCTCGCACCcgcg                                                                       >  2:523810/1‑251 (MQ=255)
                                                                                                                                                                      gggCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCTCACCcgcg                                                                       <  1:252209/251‑1 (MQ=255)
                                                                                                                                                                       ggCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACGTCGTGTATTTTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGtttgcttttcacccgcgt                                                                      <  2:118927/251‑17 (MQ=255)
                                                                                                                                                                         cATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTaa                                                                    <  1:852243/251‑1 (MQ=255)
                                                                                                                                                                               aaaCGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTa                                                              <  1:232962/251‑1 (MQ=255)
                                                                                                                                                                                 aCGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGTGTAAGCTTTAtt                                                            <  1:579554/251‑1 (MQ=255)
                                                                                                                                                                                       gTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTg                                                                                                                                         <  2:882901/168‑1 (MQ=255)
                                                                                                                                                                                       gTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTg                                                                                                                                         >  1:882901/1‑168 (MQ=255)
                                                                                                                                                                                          aCAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTTGGCAGGAAGGGTTACAACGCATCCTATCGTTTGGCAGGTATAGTCTGGAGTACCATAACCAGGTCCCCAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTg                                                   >  2:410520/1‑251 (MQ=255)
                                                                                                                                                                                               gCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTTCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGt                                              <  1:163956/250‑1 (MQ=255)
                                                                                                                                                                                                       ggTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTa                                                         >  1:626810/1‑232 (MQ=255)
                                                                                                                                                                                                       ggTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTa                                                         <  2:626810/232‑1 (MQ=255)
                                                                                                                                                                                                       ggTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCAGTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTTCCAAAACCACTGCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTGCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAAtatt                                       <  1:759625/250‑1 (MQ=255)
                                                                                                                                                                                                         tCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCa                                                                                                 >  1:143998/1‑190 (MQ=255)
                                                                                                                                                                                                         tCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCa                                                                                                 <  2:143998/190‑1 (MQ=255)
                                                                                                                                                                                                            atatAAAGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTTGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTaa                                 <  1:410520/251‑1 (MQ=255)
                                                                                                                                                                                                                aaaGATCGATTGTTGGTAGTTCAGATGACCATTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTTGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTGCTGGGCGTGAATATTATCTAATGt                              <  2:688064/250‑1 (MQ=255)
                                                                                                                                                                                                                 aaGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTTGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCt                                                                 <  2:367525/214‑1 (MQ=255)
                                                                                                                                                                                                                 aaGATCGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCt                                                                 >  1:367525/1‑214 (MQ=255)
                                                                                                                                                                                                                      cGATTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTc                        >  1:607706/1‑250 (MQ=255)
                                                                                                                                                                                                                        aTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAtt                                                                                                                         >  2:666651/1‑151 (MQ=255)
                                                                                                                                                                                                                        aTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATg                     >  2:1048189/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTTGTTGGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTGCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGAtt                                                                                                                         <  1:666651/151‑1 (MQ=255)
                                                                                                                                                                                                                           gttgGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTa                                                              >  2:198018/1‑207 (MQ=255)
                                                                                                                                                                                                                           gttgGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTa                                                              <  1:198018/207‑1 (MQ=255)
                                                                                                                                                                                                                             tgGTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTtcc                >  1:345266/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGCAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGCAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTAtt                                                       >  2:1171969/1‑211 (MQ=255)
                                                                                                                                                                                                                              ggTAGTTCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTGCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTAtt                                                       <  1:1171969/211‑1 (MQ=255)
                                                                                                                                                                                                                                    tCAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTGCCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCaa          <  1:466667/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGATGACCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCaaa         >  1:261461/1‑250 (MQ=255)
                                                                                                                                                                                                                                           aCCAGTAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGAATTa  <  2:261461/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
GGCGTTAGATACATACGAGCGCAACTGATGTCTCGCTCCTGCGGCACTCGCTCGAAACCCAAACCCGAGTAGATTGCGGCTGAACTTGATAAGGAACGCAGACCAATTTTCAAATCATGGGGTAAATCAGAATACGCGTTTAGCAATGCAGCAATGATGAATCCGAGGGCATTAGAAACGCCTGTAACAGTGCGAAGTTGGTCGATATAAAGATCGATTGTTGGTAGTTCAGATGACCAGCAATTAGGCATTTGGGCAGGAAGGGTTACAACGTATCCTATCGTTTGGCCGGTATAGTCTGGAGTACCATAACCAGGTCCCAAAACCACTTCGTGTATTGTGTATTTACTGTTGGGTTCGAAGGATTCGACATCTTGCAATAAATATTTCAGTTCTGCGGTTAACCTCGCACCCGCGTAAGCTTTATTGTATTTCTGGGCGTGAATATTATCTAATGTAGCTTCATGGTTCCATCCAAAAGAATTA  >  NC_000913/4426689‑4427174

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