Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,177 C→T E134K (GAG→AAG)  thiF ← sulfur carrier protein ThiS adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,1770CT98.7% 251.0 / NA 76E134K (GAG→AAG) thiFsulfur carrier protein ThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/45);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTCCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCC  >  NC_000913/4192937‑4193419
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
tctcTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTc                                                                                                                                                                                                                                            >  1:190881/1‑251 (MQ=255)
     ataCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCTATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCa                                                                                                                                                                                                                                        <  1:164201/250‑1 (MQ=255)
         cGCTTAATACCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTg                                                                                                                                                                                                                                   <  2:497967/251‑1 (MQ=255)
         cGCTTAATACCTTAATGGCTTCCAGTCCCTGCAAAGTGCACATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTCTCTGGCCACACGCAGCGGTAACACCACTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCCGCGGGTTGGGTGCCACGCACGCGGCATTAAT‑CTTCCGGCGAGTCGCCCTAttt                                                                                                                                                                                                                                   >  2:242339/1‑250 (MQ=255)
           cTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCACGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAAATGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTTAGTGCCCCGCAGGCCGCATTAAT‑CTTCTGGCGAGTCGCCCTATTGTc                                                                                                                                                                                                                                 >  1:479850/1‑251 (MQ=255)
                     aaaTGGCTACCAGGGCCTGAAAAGTGCCCATACCCCCGCCCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGGCGGTA‑CAGTc                                                                                                                                                                                                                       <  2:466853/249‑1 (MQ=255)
                      aaTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGGTAAGCGGCGTGTTGGGTGCCAAGCAGGCGGCATTAAT‑CTTCTGGCGAGt                                                                                                                                                                                                                                             >  2:13366/1‑228 (MQ=255)
                      aaTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGt                                                                                                                                                                                                                                             <  1:13366/228‑1 (MQ=255)
                       aTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGCGATAAGAGGCGTGTTGAGTGCCACGCAGGCGGCATTCAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCg                                                                                                                                                                                                                      >  1:409491/1‑249 (MQ=255)
                               cAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACAt                                                                                                                                                                                                            >  1:404527/1‑251 (MQ=255)
                                     cgacaaggtccccATACCCCCGACCCCCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAAGTCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGccc                                                                                                                                                                                                      <  2:362031/241‑1 (MQ=255)
                                       gCAAAGTGCCCATACCCCCGACCACCGGGCAAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGt                                                                                                                                                                                                    <  2:173833/251‑1 (MQ=255)
                                       caaaaGTGCCCATACCCCCGACCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCg                                                                                                                                                                                                     <  2:18669/248‑1 (MQ=255)
                                             tGCCCATAACCCCGACCACCGGGCCAACCACGCCAGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGGGATAAGCGGCGTGTTGGGTGCCCCGCAGGACGCATTAATACTC‑TGGCGAGTAGCCATATTGTCGGTA‑CCATTGAGCACCAACTCGGccccgtccacc                                                                                                                                                                                              >  2:244087/1‑243 (MQ=255)
                                             ggCCCATACCCCCGACCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCTTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAAcc                                                                                                                                                                                              <  2:108351/250‑1 (MQ=255)
                                              gCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGGA‑CAGTCGGGCACCACATCGGCCCGGGCAAccc                                                                                                                                                                                             >  2:286687/1‑250 (MQ=255)
                                               cccATCACCACGACCACCGCGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCACATCCGACCGCGCTGGCGGAGATAAGAGGCGTGCTGCGGGCCACGCAGGCGGCATTCAT‑CTTCTGGCGAGTCGCCATATTTTCGGTA‑CAGTCGAGCACCCCATCCGCCCGTGccaacgc                                                                                                                                                                                            >  2:418565/1‑245 (MQ=255)
                                                 cctacccccGCCCCCCGGGCCACCCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCAt                                                                                                                                                                                          <  2:320065/246‑1 (MQ=255)
                                                   tAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGGGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGGCGACATATTGTCGGTA‑CAGTCGAGCACCACCACGGCCCGTGCCACCGCATCt                                                                                                                                                                                        >  1:310810/1‑251 (MQ=255)
                                                              accGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTc                                                                                                                                                                             >  1:444106/1‑251 (MQ=255)
                                                                 gggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCAcc                                                                                                                                                                          <  1:471861/251‑1 (MQ=255)
                                                                  ggCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCACCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCACCTGACCGCCAACTCGGACCGCGCTGGCGGGGATAAGCGGGGGGTTGGGGGCCACGCAGGCGGGCTTAAA‑CTTCTGGCGAGTCGCCATTTTTGCGGTAACGGT‑GAGAACCCCATCGGCCCGCGGAACCGCCGCTTTTACCGCGTCAcac                                                                                                                                                                         >  2:260585/1‑249 (MQ=255)
                                                                      aaCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCCGCCCGTGCCACCCCCTCTTTTAACGCCTCCCCCGttt                                                                                                                                                                     >  2:261183/1‑250 (MQ=255)
                                                                      aaCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑AAGTCGAGCACCACATCGGCCCGTGCAAACGCATCTTTTAACGCCTCACCCGtt                                                                                                                                                                      >  2:108103/1‑250 (MQ=255)
                                                                      aaCCACGACCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACCGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCACCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTTGCGAGTCGCCATATTGGCGGTA‑CAGTCGAGCACCCCATCGGCCCGTGCAACCGCATCTTTTACCGCCTCAcccctgt                                                                                                                                                                     >  2:215243/1‑247 (MQ=255)
                                                                            gcccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgtt                                                                                                                                                               >  1:29253/1‑251 (MQ=255)
                                                                               cgccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgttgtt                                                                                                                                                            <  2:190881/251‑1 (MQ=255)
                                                                                 ccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGACACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCACCCCTTAACCGTtgttgt                                                                                                                                                          >  2:405546/1‑251 (MQ=255)
                                                                                             tGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAg                                                                                                                                               <  1:505599/250‑1 (MQ=255)
                                                                                              gCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATCTTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCCCATCTTTTAACGCCTCACCCGTTTACCGTTGTTGGAAACCTGTCAGtt                                                                                                                                             >  1:199166/1‑251 (MQ=255)
                                                                                               cGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTg                                                                                                                                            <  2:201966/251‑1 (MQ=255)
                                                                                                     ggCTCCTTGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCACAGGGCAGCGTCAGTACCATCACCTGACCGCCAACTCAGACCGCGCTGGCGGTGATAAGAGGCGTGTTGAGTGCCACCCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCCGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTc                                                                                                                                      >  2:67965/1‑251 (MQ=255)
                                                                                                     ggCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGt                                                                                                                                       <  1:362897/250‑1 (MQ=255)
                                                                                                     ggCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCGGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTc                                                                                                                                      <  1:540554/251‑1 (MQ=255)
                                                                                                         ccTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGGTGTTGTAATGCTGTCAGTTGAATGTCggg                                                                                                                                   <  2:273574/250‑1 (MQ=255)
                                                                                                            ggTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGGCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTc                                                                                                                               <  1:166921/251‑1 (MQ=255)
                                                                                                            ggTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTc                                                                                                                               <  1:458713/251‑1 (MQ=255)
                                                                                                                aTCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAc                                                                                                                            <  2:538686/250‑1 (MQ=255)
                                                                                                                    ggCCACAAGCAGCGGTATCACCCCTGCTCCCAGGTCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTTGCGGTGATAAGCGGCGTGTTGATGGCCTCGTAGTCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                         <  2:486952/249‑1 (MQ=255)
                                                                                                                     gCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCGGGCGAGTCGCCATATTGTCGGGA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCGTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTc                                                                                                                      <  1:242339/251‑1 (MQ=255)
                                                                                                                        ccaGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACTATCAACTGACCGCTAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCGCGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGt                                                                                                                   <  2:580307/250‑1 (MQ=255)
                                                                                                                                gggTAACACCCCTGCTTCCAGGGCGGCGTCAATACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGCTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTg                                                                                                           <  2:18036/250‑1 (MQ=255)
                                                                                                                                  gTAATACCCCTGCTCCCGGGGCGGCGTCAGTACCATCAACGGCCCGCCAA‑TCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                        <  2:391699/251‑1 (MQ=255)
                                                                                                                                    aaCACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTg                                                                                                       <  1:165896/251‑1 (MQ=255)
                                                                                                                                        ccccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGTCCGCGCTGGCGGTGATAAGCGGCTTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                    <  1:58554/250‑1 (MQ=255)
                                                                                                                                        ccccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCt                                                                                                   <  1:28411/251‑1 (MQ=255)
                                                                                                                                             gCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCg                                                                                                <  1:470163/249‑1 (MQ=255)
                                                                                                                                              cTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATTCTGTCAGTTGAATGTCGGGATTCCACTGTGTCAGTCGCTTGTGGCTGACCTGCGAtt                                                                                             >  2:325313/1‑251 (MQ=255)
                                                                                                                                               ccccAGGGCGGCGTCAGTACCAGCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGAttt                                                                                            <  2:404527/250‑1 (MQ=255)
                                                                                                                                                  cAGGGCGCCGTCAGTCCCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCgg                                                                                         <  2:208940/251‑1 (MQ=255)
                                                                                                                                                   aGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCggg                                                                                        <  1:483655/251‑1 (MQ=255)
                                                                                                                                                   aGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCggg                                                                                        >  1:94501/1‑251 (MQ=255)
                                                                                                                                                    gggCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                       <  1:189353/251‑1 (MQ=255)
                                                                                                                                                       cggcgTCAGTACCATCAACTGACCGCAAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCCCGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGGAATGCTGTCAGTTGAATGTAGGGATTCCACTGTGTTCGTAGCTGTTGGCTGACCTGCGGTTTCGGGcgat                                                                                    >  2:68383/1‑251 (MQ=255)
                                                                                                                                                        ggcgTCAGTCCCATCATCTGACCTCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                   <  2:310810/251‑1 (MQ=255)
                                                                                                                                                        ggcgTCAGTACCAACAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                   <  2:428097/251‑1 (MQ=255)
                                                                                                                                                          cgTCAGTACCATCAACTGCCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTACT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAAGCCTCCCCAGTTAACCGTTGTTGTAATGCTGTCCGTTGAATGTCGGGATTTAACTGTGTCAATCGCTGTTGGCTGACCTGCGATTTCGGGcgctcgt                                                                                 >  2:371466/1‑250 (MQ=255)
                                                                                                                                                                  ccATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCa                                                                         <  1:157993/251‑1 (MQ=255)
                                                                                                                                                                       aaCTGACCGCCCAATCCGTCCGCGATGGCGGTGATAAGCGGCGTGTTGAGTGCCTCGCCGGCGGCATTAAT‑CTTCTGGCGAGTCGTTATATTGTCGGTA‑CAGGCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTTAATGTGGGGAGGCAACTGTGTCAGTCGCTGTTGGCTGACATGCGATTTCGGGCGGTCGATATCTGCAG‑Tggc                                                                   <  2:171964/251‑2 (MQ=255)
                                                                                                                                                                                ccAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑AAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCCGTTGAATGTCGGGATTCAACTGTTTACGTCGATGTTGGCTGAACTGCGGTTTCGGGCGATCGATATATTCAGGTGG‑AAAGAGGAt                                                          >  2:129343/1‑251 (MQ=255)
                                                                                                                                                                                  aaaTCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAAGGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGAttt                                                        <  2:476656/251‑1 (MQ=255)
                                                                                                                                                                                       cgcccGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGt                                                   <  1:363913/248‑1 (MQ=255)
                                                                                                                                                                                         aCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTg                                                 <  1:491585/251‑1 (MQ=255)
                                                                                                                                                                                            gcgcTGGCGGCGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCCTATTGTCGGGA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGACTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCCGTCGCTGTTGTCTGGCCTGCGATTTTGGGCGATCGATATCTTCAG‑TGGGAAAGAGGATTTGTCGTTGCa                                               >  2:264618/1‑250 (MQ=255)
                                                                                                                                                                                                 ggCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGc                                         >  1:437028/1‑251 (MQ=255)
                                                                                                                                                                                                  gCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGGA‑CAGTCGAGCCCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCt                                        <  1:108103/251‑1 (MQ=255)
                                                                                                                                                                                                    ggttaTAAGCGGCGTGTTGATTGCCACGCAGGCGGCATGAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                      <  2:187770/247‑1 (MQ=255)
                                                                                                                                                                                                    ggTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                      <  1:286687/251‑1 (MQ=255)
                                                                                                                                                                                                      tGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTaa                                     >  1:157336/1‑250 (MQ=255)
                                                                                                                                                                                                        aTAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATg                                  >  1:70916/1‑251 (MQ=255)
                                                                                                                                                                                                                        gTGCCACGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGGAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGAACATCGTCGTCATCt                  >  2:191950/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCCAAGCAGGCGGCATTAAT‑CGTCTGGCGAGGCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGACCGTGCAACCGCATCTTTTAACGCCTCACCCGTTACCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTg                 <  1:376780/250‑1 (MQ=255)
                                                                                                                                                                                                                             aCGCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAg             <  1:405546/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCAGGCGGCATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGTATTTGTCTTTGCAGATTTCTTAAATGCACATCGTCGTCATCTGCCAGTa           <  1:164143/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ggcATTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCg     >  1:90100/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTTAAT‑CTTCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTcc  >  1:528280/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
TCTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAT‑CTCCTGGCGAGTCGCCATATTGTCGGTA‑CAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAG‑TGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCC  >  NC_000913/4192937‑4193419

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