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.8% 245.3 / NA 80E134K (GAG→AAG) thiFsulfur carrier protein ThiS adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/48);  minor base G (0/1);  total (31/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 = 9.92e-01

CTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCG  >  NC_000913/4192940‑4193421
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
ctATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGGGGTAAGCGGCGGGTTTAGTGCCCCGCAGGCGGCAATAATCTTCTGGGGAggcccc                                                                                                                                                                                                                                         >  2:15507/1‑246 (MQ=255)
  ataCCGCTTAATAACTTAATGGCTTCCAGGGCCTGCAAAGTGCCCATAACCCCGCCCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCa                                                                                                                                                                                                                                        <  1:517061/250‑1 (MQ=255)
            ataaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGt                                                                                                                                                                                                                              >  1:458464/1‑250 (MQ=255)
              aaCTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAAGCGGCATTAATCTTCTGGGGAGCCGCCCTATTGTCgggca                                                                                                                                                                                                                            >  2:10639/1‑247 (MQ=255)
                  tAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCCGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGGCGCCCTATTGTCGGTACAGTcc                                                                                                                                                                                                                       >  1:227621/1‑250 (MQ=255)
                  aaaTGACTTCCAGTGCATGCAAATGCCCCATCACCCCGACCCCGTGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGGGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTc                                                                                                                                                                                                                        <  2:229667/249‑1 (MQ=255)
                    attgCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTAATAAACGGCGGGGTGTGTGGCACGCAGGCGGGCATAAATTTATGGCGGGTCGCCATATTTTCGGTActcgcg                                                                                                                                                                                                                       >  2:475575/4‑244 (MQ=255)
                     tGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATGAGCGGCGTGTTGAGTGCCACGCAGGCGGGATTAATCTTCTGGCGAGTCGCCATATTGTCAGGACCGTCGAGc                                                                                                                                                                                                                    >  1:207777/1‑251 (MQ=255)
                      gctccccgggtctggataggctcccATCACCCCGATCCCCGCGCCACCCCGGCCCGCCGGGCGTCCGTTGTGTTCTGGTTCCGGTTTCTCCGGCCATAGGCAGCGATAACCCCCCTGCCTCCAGGGCGGCGTCAGTACCATCAACTTACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTTTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGGACAGTCGAGc                                                                                                                                                                                                                    <  2:423082/228‑1 (MQ=255)
                           ccatggCCTGCAAAGTGCCCATAACCCCGACCCCGGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGTCGCCATGACCATCACCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACcaca                                                                                                                                                                                                              <  2:480278/246‑1 (MQ=255)
                            cAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTTGCGAGTCGCCATATTGTCGGTACCGGCGAGcaccccct                                                                                                                                                                                                             >  1:17786/1‑247 (MQ=255)
                              gTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGAGGTGATAAGCGGAGGGTTTAGTGCCCAGCAGGAGGCATTAATCTTCTGGCCAGTCGCCCTAttttacggaaagtcgggcccccaatcg                                                                                                                                                                                                           >  2:395469/1‑226 (MQ=255)
                               tGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCgg                                                                                                                                                                                                          >  1:381772/1‑251 (MQ=255)
                               cccacgccgtcctcccgtcctccccGCCCCCCGGCCCACCCACGCCCGCCGTGCGGCAGGTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCAGCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCg                                                                                                                                                                                                           <  2:214041/230‑1 (MQ=255)
                                    gCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGccc                                                                                                                                                                                                       <  1:552921/249‑1 (MQ=255)
                                     cAAAGTGCCCATACCCCCGCCCCCCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTg                                                                                                                                                                                                    <  1:481538/251‑1 (MQ=255)
                                       aaagcacctaaatcccggccccccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCTTCCTGGCCCCTGGTTATCTGGCCACAGGCAGCGGTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGTCCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGc                                                                                                                                                                                                   <  1:214609/229‑1 (MQ=255)
                                           gCCCATAACCCCGCCCACCGGGCCACCCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCTTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCg                                                                                                                                                                                              <  1:41528/251‑1 (MQ=255)
                                                  ccccccGCCCCCCGGCCCACCCCCGCCGGCCGTGCCGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCCCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGATTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCTAGCACCACATCGGCCCGTGCAACCGCATCtt                                                                                                                                                                                        <  2:240885/249‑1 (MQ=255)
                                                  aaCCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAAGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAACGGCCCGTGCAACCGCAtatt                                                                                                                                                                                        >  1:497584/1‑248 (MQ=255)
                                                    ccccGACCACCGGGCCAACCCCGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCtttt                                                                                                                                                                                      <  1:456488/251‑1 (MQ=255)
                                                    ccccGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAACGGCGTGGTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATCTTGTCGGGACAGTCGAGCACCACATCg                                                                                                                                                                                                           >  1:214041/1‑230 (MQ=255)
                                                        gACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAcc                                                                                                                                                                                  >  1:48067/1‑250 (MQ=255)
                                                        gACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACAACCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGAGGAGATAAGCGGCGTGTTGAGTGCCCAGCAGGCGGCATTAATCTTCTGGCGCGTCGCCATATTGTCGGTACAGACGAGCACCACATCAGCCCGTGCAACCGCATCTTTTAACg                                                                                                                                                                                  >  1:531523/1‑251 (MQ=255)
                                                           caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCCCATCGGccc                                                                                                                                                                                                       >  1:29284/1‑227 (MQ=255)
                                                           caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAAGccc                                                                                                                                                                               >  1:279254/1‑250 (MQ=255)
                                                           caccGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGATGAGTGCCACGCAGGCGGCATTAATCCTCTGGCGAGTCGCCATCTTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCt                                                                                                                                                                               >  1:61628/1‑251 (MQ=255)
                                                                       ccccgcccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTCCCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAc                                                                                                                                                                    <  1:238912/248‑1 (MQ=255)
                                                                        cACGCCCGCCGTGCGGCAGTTGAGTTCCGGCTCCGGGTTATCTGGCCACAGGCAGCGTTAACCCCCCTGCTCCCAGGGCGGCGTCAGTACCACCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGACCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAAcc                                                                                                                                                                   <  2:381772/251‑1 (MQ=255)
                                                                            cccgccGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCggtt                                                                                                                                                               >  1:227692/1‑248 (MQ=255)
                                                                             ccgccgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGTTGAGTGCCACGCAGGCGGCCTTAATCCTCTGGCGAGTCGCCATATTGTTGGGACAGTCGCGCCCCACATAGGCCCGGGACACCCCATCTTTTAAAGCCCCACCCCTTaaacacggt                                                                                                                                                              >  2:142642/1‑244 (MQ=255)
                                                                                ccgTGCGGCAGTTGCGTTCTGGCTCATGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTtgttgt                                                                                                                                                           <  1:314791/251‑1 (MQ=255)
                                                                                 cgTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTTACCGTtgttgtt                                                                                                                                                          >  1:506065/1‑250 (MQ=255)
                                                                                    gCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATg                                                                                                                                                       <  1:20588/251‑1 (MQ=255)
                                                                                     cGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTTTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATgg                                                                                                                                                      >  1:541315/1‑250 (MQ=255)
                                                                                            tGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAg                                                                                                                                                <  1:244012/250‑1 (MQ=255)
                                                                                                    ggCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGt                                                                                                                                        <  1:383471/250‑1 (MQ=255)
                                                                                                        actggtgaccTGGCCCCCGGCAACGGTCACCCCCCTGTTCGCAGGGCGGCTTCAGGACTATCATCGGAACGCCTAATCCGACCGCGCTGGCGGGTATAAGTGGCGTGTTGAGTTCCACGCTGGCGGCATTAAGCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCATCCGCATCTTTTAACGCCGCACCCGTTAACCGTTCTTGTAATGCTGTCAGTTGAATGGCgg                                                                                                                                     <  2:470259/240‑1 (MQ=255)
                                                                                                         cTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGCGCACCCCATCGGCCCGGGCCACCGCATCTTTTTACGCCTCCCCCGTTTACCGTTGTTGTTATGCTGTCCGTTTGATGTCGGGAt                                                                                                                                  >  2:53277/1‑251 (MQ=255)
                                                                                                          tctttaccgggctcctggcgctggTAACCCCCATGCTCCCTGGGAGTCGTCAGTACCTACGCCTTACCTCCTAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGTGTGCTACGAAGGCGGCATTAATCTTCTGGCTGGTCGCCATATTGTAGGTGTAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGGTGAATGTCGGGAt                                                                                                                                  <  2:500769/228‑1 (MQ=255)
                                                                                                           ggTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTc                                                                                                                                <  1:537314/251‑1 (MQ=255)
                                                                                                               agcgggCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGACAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAAc                                                                                                                             <  1:148213/246‑1 (MQ=255)
                                                                                                               aTCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGGTAAGCGGCGTGTTGAGTGCCCCGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTTTCGGGACAGTCGGGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCCCACCCCTTAACCCTTGTTGGAATGCTGTCAGTTGAATTTCGGGGTTccaat                                                                                                                            >  2:213881/1‑247 (MQ=255)
                                                                                                                     ccacctgccatggTAACACCCCGGCGCGCAGGGCGGCGGCAGTACCATCAACGGACCGCCAAACCCGACCGCGCTGGCGGTGAGTAGTTGCGTGTTGAGTGCCACGCAGGCTGCATTAATCTGCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACAGCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTTACCGTTGTTCTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCa                                                                                                                      <  2:227621/240‑1 (MQ=255)
                                                                                                                     ccACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCCGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGGGTCGCCATATTGTCGGTACAGTCGAGCACCAACTCGGCCCGAGCAACCGCCTCTTTTAAACACTCCCCCGTTCACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGGTTCACCTGTGTc                                                                                                                       >  2:344980/1‑250 (MQ=255)
                                                                                                                                   aaCACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAACGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCCTCTTTTAACGCCTCCCCCGGTAACCGTTGTTGTAATGCTTTCAGTTGTAAGTCGGGATTCAACTTTGTCAGTCGCTTTTTGCTg                                                                                                        >  1:552827/1‑251 (MQ=255)
                                                                                                                                             cTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAATCGGACACCACATTGGACCGTGCAAACGCCTCTTTTAAAGCCTCACCCGTTAACCGTTTTTGTAATGCTGTCAGTTGTATGGCGGGATTCCACTGTGTccggctcgtgttgattgactgctgtt                                                                                              >  2:413334/1‑226 (MQ=255)
                                                                                                                                              tCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCAGCGTGTTGGGTGCCACGCCGGCGGCATTAATCTTCTGGCGAGTCGCCATAATGTCGGTTCAGTCGGGCCCCCCATCCGCCCCTGCCACCGCATCTTTTTACGCCTCACCCGTTAACCGTTTGTGTAATGCTGTCAGGTGAATGTAGGGATTCAACTTTGGTAGTCGCAGTTGGccgtccctggaatt                                                                                             >  2:497747/1‑238 (MQ=255)
                                                                                                                                               cccAGGGCGGTGTAAGTACCATCAACTGACCGCCAAAGCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTCCAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTc                                                                                            <  2:92640/251‑1 (MQ=255)
                                                                                                                                                ccAGGGCCGCGTCAGTACCATCAACTGACCCCCAAATCCGACCGCGCTGGTGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGGCGCTGTTGGCTGACCTGCGATTTCg                                                                                           <  2:313151/251‑1 (MQ=255)
                                                                                                                                                   gggCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                        <  1:53277/251‑1 (MQ=255)
                                                                                                                                                   gggCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGc                                                                                        <  1:213881/251‑1 (MQ=255)
                                                                                                                                                          ggcAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGcc                                                                                                                                                                                <  2:214974/154‑1 (MQ=255)
                                                                                                                                                          gTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGcc                                                                                                                                                                                >  1:214974/1‑156 (MQ=255)
                                                                                                                                                              ggtcaaTCAACTGACTGCCACATCAGACCGCGTTGGCGGTGATAAGCGGCGTGTTGAGGGCCCCGCAGGCGGAATTAATCTTCTGTCGAGTCGCCATATTGTCGGCACAGTCGAGCACCACATCGGCCCGTGCACCCGCATCTTTTAACGCCTCACACGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCt                                                                             <  2:143210/246‑1 (MQ=255)
                                                                                                                                                               tACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACtg                                                                                                                           <  1:51509/204‑1 (MQ=255)
                                                                                                                                                                tccattcattgcccGCCAAAACCTACCGCTCTGGCGGTGCTAAGCGGTGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGCGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGGTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTCCGGGCGATCGATATCtt                                                                            <  2:207777/238‑1 (MQ=255)
                                                                                                                                                                 ccATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGGGTTGAGTGCCACGCAGGCGGCCATAAACTTCTGGCGAGTCGCCAAATTGTCGGTACACGCCAGCACCACATCGGCCCGAGCAACCGCATCTTTTTACGCCTCACCCGATAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGTTTCAACTGTGTCAGTCGCTATTTGCGGACCTGCGATTTCGGGGGATCGTTATCttt                                                                           >  2:435688/1‑249 (MQ=255)
                                                                                                                                                                  cATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAg                                                                         <  1:151038/251‑1 (MQ=255)
                                                                                                                                                                  aaTCAACGGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGTGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCTCATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTATCCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAg                                                                         <  1:50895/250‑1 (MQ=255)
                                                                                                                                                                      aaCTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGt                                                                     >  1:440516/1‑251 (MQ=255)
                                                                                                                                                                           actcccAAATTCGACCGCGTGGTCGGGGATAAGCGGCGTGTTGAGTGCCACGCAGGCTGCATTAAGCTTCTGGCGAGTCGCCAGAGTGTCGGGTCAGTCGAGCACCACTTCTGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAGCCGGTGTTGTAATGCTGTCAGTTGAATGTCGGGAGTCAATTGTGTTAGTCGCTGTTGGCTGAGCTGCGATTTCTGGCGATCTAGATCTGCAGTGGGAAaga                                                                <  2:240506/247‑1 (MQ=255)
                                                                                                                                                                             cGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAgg                                                              <  1:503380/251‑1 (MQ=255)
                                                                                                                                                                              cgcAAATCCGACCGCGCGGGCGTGGATAAGCGGGGTTTGTAGTGCCAGGCAGGCGGCATTAATCTTCTGGCGTGTCGCCAGATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTTAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGa                                                             <  2:398302/249‑1 (MQ=255)
                                                                                                                                                                                 aaatcctcccGCGCTGGCGGTGATAAGGGGCGTGTGGAGTGCCACGCAGGCGGCAATAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCGCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGAttt                                                          <  2:458464/243‑1 (MQ=255)
                                                                                                                                                                                  aaTCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTg                                                         <  1:10639/251‑1 (MQ=255)
                                                                                                                                                                                     ccgtccGCGCTGGCGGTGATAAGCGGCGTGTTGATTTCCACGCAGTCCGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGGCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCGCACCCGTTAACCGTGGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTc                                                       <  2:552827/246‑1 (MQ=255)
                                                                                                                                                                                       gACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTTTAATGCTGTCAGTTGAATTTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGGTATCTTCCGTGGTAAAGAGGATTTTTCGtt                                                    >  1:17021/1‑251 (MQ=255)
                                                                                                                                                                                       gACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGtt                                                    <  1:437501/251‑1 (MQ=255)
                                                                                                                                                                                         ccGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGc                                                  <  2:48067/251‑1 (MQ=255)
                                                                                                                                                                                          cgcgCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCa                                                 <  1:395469/251‑1 (MQ=255)
                                                                                                                                                                                             gcTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGAt                                              <  1:12740/251‑1 (MQ=255)
                                                                                                                                                                                                ggCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGAtt                                             <  1:103500/249‑1 (MQ=255)
                                                                                                                                                                                                              ggcgTGTTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAAAGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGACAGTCGCTGTTTGCTGGCCTGCGATTTCGGGCGATCGATATCTTTCGTGGTCAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtc                              >  2:305212/1‑251 (MQ=255)
                                                                                                                                                                                                                 gtgtTGAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACAtcgtc                           <  1:405444/251‑1 (MQ=255)
                                                                                                                                                                                                                      gAGTGCCACGCAGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCa                       <  1:497747/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTGCCACGCAGGCGGCATTAAGCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCt                    <  1:571589/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGGCGGCATTAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTAcc           <  1:423516/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tAATCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTccc   <  1:15507/250‑1 (MQ=255)
                                                                                                                                                                                                                                             aaTCTTCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCg  <  1:381563/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
CTATACCGCTTAATAACTTAATGGCTTCCAGTGCCTGCAAAGTGCCCATAACCCCGACCACCGGGCCAACCACGCCCGCCGTGCGGCAGTTGCGTTCTGGCTCCTGGTTATCTGGCCACAGGCAGCGGTAACACCCCTGCTCCCAGGGCGGCGTCAGTACCATCAACTGACCGCCAAATCCGACCGCGCTGGCGGTGATAAGCGGCGTGTTGAGTGCCACGCAGGCGGCATTAATCTCCTGGCGAGTCGCCATATTGTCGGTACAGTCGAGCACCACATCGGCCCGTGCAACCGCATCTTTTAACGCCTCACCCGTTAACCGTTGTTGTAATGCTGTCAGTTGAATGTCGGGATTCAACTGTGTCAGTCGCTGTTGGCTGACCTGCGATTTCGGGCGATCGATATCTTCAGTGGTAAAGAGGATTTGTCGTTGCAGATTGCTTAAATGCACATCGTCGTCATCTGCCAGTACCAGCGTCCCG  >  NC_000913/4192940‑4193421

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