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,1770CT97.6% 262.2 / ‑6.2 85E134K (GAG→AAG) thiFsulfur carrier protein ThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/37);  minor base A (0/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGT  >  NC_000913/4192938‑4193417
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
ccctATACCGCTTAATAACTTAATGGCTTCCAGGGCCTGCAAAGGGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTc                                                                                                                                                                                                                                         <  1:313971/248‑1 (MQ=255)
  ctATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTc                                                                                                                                                                                                                                         <  1:329025/248‑1 (MQ=255)
    ctaCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAACGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATCTTCTGGCGAGTCGCCat                                                                                                                                                                                                                                    >  2:83776/2‑251 (MQ=255)
    ataCCTCTTAATAACTTAAGGGCTTCCAGTGCCTGGAAAGGGCCCATAACCCCGACCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCa                                                                                                                                                                                                                                     <  1:475423/250‑1 (MQ=255)
    ataCCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAACGAGCCCATAACCCCGACCACAGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCCGCGGGAACACCCCTGCTCCCAGGGCGGCGTCCGGACCAGCAACTGACCGCCAAATCCGCCCGCGCTGGCGGTGAAAAACGGCGCGGTCGGCGCCACGCAGGCGGAATTCATCCTCTGGCgcggacccct                                                                                                                                                                                                                                    >  2:426293/1‑242 (MQ=255)
      cccGCTTAATAACTTAATGCCTCCCAGTTCCTGCAAAGGGCCCATAACCCCGCCCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCatat                                                                                                                                                                                                                                  <  2:491600/250‑1 (MQ=255)
        cGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTg                                                                                                                                                                                                                                >  1:57329/1‑251 (MQ=255)
        cGCTTAATAACTTAATGGCTTCCAGTCCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGATGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTg                                                                                                                                                                                                                                <  1:71446/251‑1 (MQ=255)
             cataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGGTGAGGGCCACGCAGGCGGCATTAATCTTCTGGCGCGGCGCCATATTGTCggg                                                                                                                                                                                                                           >  2:265652/2‑250 (MQ=255)
             aaaaccTTAATGGCTTCCGGTGCCTGCAAAGTGCCCATACCCCCGCCCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGt                                                                                                                                                                                                                           <  2:207052/246‑1 (MQ=255)
               taaCTTAATGGCATCCAGTGCATGCAAAGTGCCCATAACCCAGCCCACCGGGCCAACCACGCCCGCCGAGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCCGTACCATCAACTGACCGCCAAATCCGACCGCGATGGCGGTGATAAGCGGCGTGTTGAGGGCCACGCAGGCGGCATTAATCTTCTGGCGATTCGCCATATTTTCGGTAc                                                                                                                                                                                                                         >  2:44637/1‑251 (MQ=255)
               aaaCTTAATGCTGTCCAGTGCTGGCAAGGTGCCCACACCCCCGCCGTCCGGGCCAACCCCGCCCGCCGGGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGGAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGGTAAGCGTCTTGTTGAGTGCCACGCAGGCGGCAGTAATCTTCTGGCGAGTCGCCATATTGTCGGTAc                                                                                                                                                                                                                         <  2:540456/250‑1 (MQ=255)
                aacgtAATGGCTTCCAGGGCCTGCAAAGTGCCCATACCCCCGACCACCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGGCGGTACa                                                                                                                                                                                                                        <  2:39011/247‑1 (MQ=255)
                     aaTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCCGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGGCGGTACAGTCgg                                                                                                                                                                                                                   >  1:331849/1‑250 (MQ=255)
                      aTGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGCCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCCGGCGGCATTAATCTTCTGGCGAGGCGCCATATTGTCGGGACatgcg                                                                                                                                                                                                                    >  2:584411/1‑245 (MQ=255)
                          cTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCAGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGGGATAAGCGGCGTGTTGGGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGGCGGTAACGACGAGcacc                                                                                                                                                                                                              >  1:450452/1‑251 (MQ=255)
                                 gttccagCAAAGTGCCCATAACCCGGCCCGCCGGGCCCACCAGGCCCGCCGGGCAGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAATACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGGTGAGTTCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATGGCCGGTACAGTCTAGCTCCACATCg                                                                                                                                                                                                       <  2:307628/245‑1 (MQ=255)
                                  tGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGCGCACCCCATCgg                                                                                                                                                                                                      >  1:355256/1‑251 (MQ=255)
                                     cTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGccc                                                                                                                                                                                                   >  1:81257/1‑251 (MQ=255)
                                      tGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACAATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATCAGCGGCGTGTTGAGTGACCCGCAGGCGGCCATAATCTTCTGGCGAGTCGACATATTGTCGGTACAGTCGAGAACCACATCGGcacg                                                                                                                                                                                                  >  1:45095/1‑248 (MQ=255)
                                        cAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCCCATCGGCCCGTg                                                                                                                                                                                                >  1:158149/1‑251 (MQ=255)
                                           aGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCCAATCCGACCGCGCTGGCGGCCATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGTGTCGCCATATTGTCGGTACAGTCGAGCTCCACCACCGCCCGTGcc                                                                                                                                                                                              >  1:303086/1‑249 (MQ=255)
                                             tGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAcc                                                                                                                                                                                           >  1:168195/1‑251 (MQ=255)
                                             tGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTAAAATCAACTGACCGCCAAATCCGACAGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCCTTAATATTCTGGCGAGTCGCCATATTGTCGGTAAAGTCGAGCACCACATCGGCCCGTGCAAcc                                                                                                                                                                                           >  1:93827/1‑251 (MQ=255)
                                                 cATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACCGTCGAGCACCACATCGGCCCGTGCAACCGCAt                                                                                                                                                                                       >  1:447068/1‑251 (MQ=255)
                                                 cATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGAGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCAt                                                                                                                                                                                       >  1:531555/1‑251 (MQ=255)
                                                     aCCCCGACCACCGGGCCGACCCCGCCCGCCGTGCGGCTGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATAAACTGACCGCCAAATCCGAACGCGCTGGCGGTGGTAAGCGGCGGGTTGAGTGCCACGCGGGCGGCATTAAACTTCTGGCGGGGCGCCATATTGTCGGGCCAGGCGTGCACCTCATCGGcgcgagcgaccgagccttt                                                                                                                                                                                   >  2:543470/1‑233 (MQ=255)
                                                       cccGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGGGATAAGCGGCGTGTTGGGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTAAAGTCGAGCACCAACACGGCCCGTGCAACCGCATCttttt                                                                                                                                                                                 >  2:209326/1‑250 (MQ=255)
                                                        ccGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCCCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCACGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCACGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTa                                                                                                                                                                                 >  1:268714/1‑250 (MQ=255)
                                                                         cACGCCCGCCGTGCGGCAGTTGCGTTCTGGCGCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGGCAGTACCATCAACTGACCGCCAAATCCGAACGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAcc                                                                                                                                                               <  1:74310/251‑1 (MQ=255)
                                                                            gcccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgtt                                                                                                                                                            <  1:485705/251‑1 (MQ=255)
                                                                            gcccgccGTGCGGCAGTTGCGTTCTGCCTCCTGTTTACCTGGCCACAGGCAGCGGTAACACTCCTGGTACCAGGGCGGCTGCAGTACCATCAACTGACCGCCCAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGGCGGTACAGGCGAGTACCATATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCgtt                                                                                                                                                            <  1:357071/251‑1 (MQ=255)
                                                                                 ccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGGGTCGCCATATTGTCGGTACAGGCGAGCCCCCCATCGGCCCGTGCAACCGCCTCTTTTAACCCCTCACCCGTTTACCGAtgttgt                                                                                                                                                       >  2:155925/1‑251 (MQ=255)
                                                                                          aGTTGCGTTCTGGCTCCTGGTTACCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCGGGGCGGCGTTAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGTGGCATTAATCTTCTGGCGAGTCTCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCATCCGCATCTTTTAACGCCTCACCCGTTAGCCGTGGTTGTAATGCTGTc                                                                                                                                              <  1:291267/251‑1 (MQ=255)
                                                                                             tGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAg                                                                                                                                            <  1:482044/250‑1 (MQ=255)
                                                                                               cGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGGTGTAATGCTGTCAGTTg                                                                                                                                         <  1:478932/251‑1 (MQ=255)
                                                                                                gTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGa                                                                                                                                        <  1:584411/251‑1 (MQ=255)
                                                                                                    tGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAAAGCCTCACCCGTTAAACGTTGTTGTAATGCTGTCCGTTGAATgg                                                                                                                                    >  2:313896/1‑250 (MQ=255)
                                                                                                     ggCTCCTGTGTATCTGGCCACAGGCAGCGGTAACGCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTTACCGCCAAATCCGACCTCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCCCGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCTCCTCATCGGCTCTTGCAACTGCATCTTTTAACTCCACACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAAt                                                                                                                                      <  2:270579/248‑1 (MQ=255)
                                                                                                       cTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGCTGAATGTCgg                                                                                                                                 >  2:397842/1‑251 (MQ=255)
                                                                                                         ccTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGCGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATTCTGTCAGTTGAATGTCggg                                                                                                                                >  1:581882/1‑250 (MQ=255)
                                                                                                                  cTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCAACACATCGGCCCGTGCAACCGCATATTTTAAAGCCTCACCCGTTAACCGTTGGTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtgt                                                                                                                      >  1:541795/1‑251 (MQ=255)
                                                                                                                     gCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTCATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTCACCGTTGTTGTACTGCTGTCCGTTGAATGTCGGGATTCCACTGTGTc                                                                                                                   >  1:367143/1‑251 (MQ=255)
                                                                                                                      ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCCGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTTGCGAGTCGCCATATTGTCGGGACAGTCGAGCACCCCATCGGCCCGTGCAACCGCCACTTTTAAAGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGCTTCCACTGTGTc                                                                                                                   >  2:118586/1‑250 (MQ=255)
                                                                                                                                cGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCCTCAACTGACCGCCAAATCCGACCGCGCTGGCGGAGATAAGCGGCGTGTTGCGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTCCAGTCGGGCCCCACATCGGCCCGTGCAACCGCATCTTTTTACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTTTCCTTCGCTGTTg                                                                                                        >  2:477577/1‑251 (MQ=255)
                                                                                                                                   tAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGCTAAGCGGAGCGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTAGAGCACCCCATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                     >  1:387917/1‑251 (MQ=255)
                                                                                                                                   tAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                     >  1:184255/1‑251 (MQ=255)
                                                                                                                                   tAACACCCCCGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTAAAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCACCCGGTAACCGTTGTTGTAATGCTTTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCt                                                                                                     >  2:565951/1‑251 (MQ=255)
                                                                                                                                    aaCACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTg                                                                                                    <  1:367374/251‑1 (MQ=255)
                                                                                                                                       acCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGAcc                                                                                                 >  1:443153/1‑251 (MQ=255)
                                                                                                                                         cccTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTg                                                                                               <  1:139808/251‑1 (MQ=255)
                                                                                                                                            tGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGGTGGCTGACCTGCgg                                                                                            >  1:560843/1‑250 (MQ=255)
                                                                                                                                             gCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCg                                                                                             <  2:330303/249‑1 (MQ=255)
                                                                                                                                             gCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGAt                                                                                           <  1:10564/251‑1 (MQ=255)
                                                                                                                                                        ggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatc                                                                                <  2:42238/251‑1 (MQ=255)
                                                                                                                                                        ggcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAATGTGTCAGTTGCTGTTGGCTGACCTGCGATTTCGGGcggtc                                                                                >  2:41023/1‑248 (MQ=255)
                                                                                                                                                         gcgTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGGTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGGAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcg                                                                               >  1:435016/1‑251 (MQ=255)
                                                                                                                                                           gTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTGAACGCCTCACCCGTTAACCGTGGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcgat                                                                             <  1:250382/251‑1 (MQ=255)
                                                                                                                                                               gTACAATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCt                                                                         <  1:23158/251‑1 (MQ=255)
                                                                                                                                                                      cAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATCTTGTCGGGACAGTCGAGCACCACCTCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGGAATGCTGTCAGTTGAATGGCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGcgatcggtacattcagggg                                                                  >  2:389577/1‑238 (MQ=255)
                                                                                                                                                                      cAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAAGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCATATGAATGTCCGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGAGATCTTCCGTgg                                                                  >  2:160309/1‑251 (MQ=255)
                                                                                                                                                                       caCTGACCGCCAAATCCGACCGCGCTGGCGGTGTTCAGCGGCGTGTTGTGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCTCCATATTGTCGGTACAGTCGAGCACCACATCGTCCCGTGCATCCGCATCTTTTAACGCCTCTCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGTTATCTTCAGTGGt                                                                 <  2:481181/250‑1 (MQ=255)
                                                                                                                                                                       aaCTGACCGCCAAATCCGACCGCGCTGGCGGGGATAAGCGGCGTGTGGAGTGCGACGCAGGCGTCAGTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGGTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCATTCACTGTTGGCTGACCTGCGTTATCGCGCGACGGATATCTTCAGTGGt                                                                 <  2:586147/251‑1 (MQ=255)
                                                                                                                                                                        aatgcccGCCAAATCCGCCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGGCGAGCCCCACAGCGGCCCGTGCAACCGCATCTTTTAACGCCTCCCCCGTTAACCGTTGTTGTAATGCTGTCAGGTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTa                                                                <  2:450452/246‑1 (MQ=255)
                                                                                                                                                                            aCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAAACTTCTGGCGAGTCGCCATATTGTCGGTAAAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGAGTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATGTCGGGCGATCGATATCTTCAGTGGTAAaga                                                            <  2:142325/251‑1 (MQ=255)
                                                                                                                                                                                  aaatctgcccGCGCTGGCGGCGATACGCGGCGTGTTGAGAGCCACGCCCGCGGCATTAATCTTCTGGCGCGTCGCCATATTCTCGGTACAGTCGCGCACCCCAACGGCCCCGGCCACCGCCTCTTTTAAAGCATCCCCCGTTAAACGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTTAACTGTGTCAGTCGCTGTTGGCTGTCCTGCGTATTCGTGCGATCGATCTCTTCCGTGGTTAAGAGGAttt                                                      >  2:34850/9‑251 (MQ=255)
                                                                                                                                                                                  aaaTCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAttt                                                      <  1:558642/251‑1 (MQ=255)
                                                                                                                                                                                       cGACCGCGCTGGCGTTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGt                                                 <  1:62475/251‑1 (MQ=255)
                                                                                                                                                                                         cccGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                               <  1:482592/250‑1 (MQ=255)
                                                                                                                                                                                         aCCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGAGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTTAGTTGAATGTCGGGATTCCAATGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                               >  1:542230/1‑251 (MQ=255)
                                                                                                                                                                                         aCCGCGCTGGCGGTGATAAACGGCGTGATGCGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGGACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGGCCTGCGATTTCGGGCGATAGATATCTTCAGTGGTAAAGAGGATTTGTCGTTg                                               >  2:529443/1‑251 (MQ=255)
                                                                                                                                                                                                    ggTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTa                                    <  2:81257/251‑1 (MQ=255)
                                                                                                                                                                                                       gATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCCGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCGTAAAt                                 <  2:39958/251‑1 (MQ=255)
                                                                                                                                                                                                           aGCGGCGTTGTGAGTGCCACGCAGGCGGCATCAATCTACCGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAAACGCATCTTTTAACGCCTCTCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGTTCGATAACTCCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                             <  2:187600/251‑1 (MQ=255)
                                                                                                                                                                                                           aGCGGCGTGTTGAGTGCCCCGCAGGCGGCACTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGcac                             <  2:142375/251‑1 (MQ=255)
                                                                                                                                                                                                            gcggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTACTGGCGAGGCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGAGTGCTTAAATGcaca                            <  1:433741/251‑1 (MQ=255)
                                                                                                                                                                                                                cgTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTTGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgt                        >  2:284122/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtgtTGAGTGCCACGCAGGGGGCATTAATCCTCTGGCGAGGCGCCATATTGTCGGTACAGTCGAGCGCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTTGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtc                       <  1:507632/251‑1 (MQ=255)
                                                                                                                                                                                                                   gttgacggCCACGCAGGCGGCAGTAATCTGCTGGCGAGGCGCCATATTGTCGGTACAGTCGAGCTCCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCTCCCGTGAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATcgtcgt                     <  1:260498/244‑1 (MQ=255)
                                                                                                                                                                                                                   gatGAGTGCCACGCAGACGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCAAGCACCACATCGGCCCGTGCAAACGCATCTTTTAAAGCCACACCCGTTAACCGTAGTTGTAATGCTGTCAGTTGAATGTCGGGATACAACTGCGTCAGTCGCTGTTTGCTGAGTTGCGATTTCGGGCGATCGGAATCTTCAGTGGTAAAGAGGATGTGTCGTTGCAGACTGCTTAAATGCACATcggcgt                     >  2:443722/3‑251 (MQ=255)
                                                                                                                                                                                                                       aGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATc                 >  1:4910/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAg           <  1:44637/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTAcc       >  1:327614/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cggcATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTACAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGACATCTGCCAGTACCAGc    >  1:451066/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gcATTAATCTTCTGGCGAGTCGCCATATTGTCGTTACATTCGAGCACCACATCGGTCCGTGCAACCGCATCTTTTAACGCCTCACCCTTTAACCGTTGTGGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCTATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGAACATCGTCGTCATCTGCCAGTACCAGCGt  <  2:462430/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
CTCTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGT  >  NC_000913/4192938‑4193417

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