Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,193,996 T→G D70A (GAT→GCT)  thiE ← thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,193,9960TG98.7% 284.4 / NA 79D70A (GAT→GCT) thiEthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (40/38);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCC  >  NC_000913/4193751‑4194239
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                 <  1:174407/250‑1 (MQ=255)
 ccttGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGGCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGGCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCAGCCCCAAATGGACGCCATAGGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGa                                                                                                                                                                                                                                               <  1:31936/248‑1 (MQ=255)
           tGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcc                                                                                                                                                                                                                                     >  2:316591/1‑250 (MQ=255)
               ccAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATAGCGTCATGTGTCGAAACGCCCAACCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAAGCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATcgcgcg                                                                                                                                                                                                                                 <  2:404226/251‑1 (MQ=255)
                  gCCCCTGCGGTGAAGAAGGTATCTGTTTGGTTTGCGTCGGGAACACATTTCCCAGCCCGAAATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCGTGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAAGcgcgcggta                                                                                                                                                                                                                              <  1:156611/251‑4 (MQ=255)
                    cccTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTtat                                                                                                                                                                                                                             >  2:251767/1‑250 (MQ=255)
                          ggTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcggcg                                                                                                                                                                                                                      <  1:43072/251‑1 (MQ=255)
                              cAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCCGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                                                                                                                                                  >  2:554178/1‑251 (MQ=255)
                                gAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                                                                                                                                                  >  1:214116/1‑249 (MQ=255)
                                 aaGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAgc                                                                                                                                                                                                               >  2:103479/1‑251 (MQ=255)
                                     cATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCa                                                                                                                                                                                                            >  1:372514/1‑250 (MQ=255)
                                          gtttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGACCAGCGCAATTgcc                                                                                                                                                                                                      >  2:109087/1‑251 (MQ=255)
                                           tttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCTCGTTATAGCGGCGGCCCAGCGCCATTgccc                                                                                                                                                                                                     >  2:504194/1‑250 (MQ=255)
                                            ttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccgc                                                                                                                                                                                                    >  1:404499/1‑251 (MQ=255)
                                                tttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACg                                                                                                                                                                                                >  1:109088/1‑251 (MQ=255)
                                                     gTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATc                                                                                                                                                                                           >  2:448199/1‑251 (MQ=255)
                                                        gggAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCACTTGCCGCCCCGACATCGGc                                                                                                                                                                                        >  2:451118/1‑251 (MQ=255)
                                                         ggAACACAGGTCCCAGCGCGATATAAGAGGGGGGGGCTGCCAGCGCGACGTCTATTTCCATAGCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCt                                                                                                                                                                                       <  2:2857/251‑1 (MQ=255)
                                                              acaTGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAc                                                                                                                                                                                  >  1:445395/1‑251 (MQ=255)
                                                                  gTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACctct                                                                                                                                                                              <  1:209528/251‑1 (MQ=255)
                                                                            gATATCAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGcc                                                                                                                                                                    >  1:472828/1‑251 (MQ=255)
                                                                             atatAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCg                                                                                                                                                                   <  1:478474/251‑1 (MQ=255)
                                                                               ataAGGGGGGCGTGCTGCCAGCGCGCCTTCGATTTCCATATCGTCATGTGTCGAACCCCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGAt                                                                                                                                                                 <  2:157999/251‑1 (MQ=255)
                                                                               ataAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGAt                                                                                                                                                                 >  1:430849/1‑251 (MQ=255)
                                                                                  agagGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCtt                                                                                                                                                              >  1:343021/1‑251 (MQ=255)
                                                                                  agagGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCtt                                                                                                                                                              >  1:583985/1‑251 (MQ=255)
                                                                                    gggggCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGCCGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTg                                                                                                                                                            <  1:514614/250‑1 (MQ=255)
                                                                                       ggCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                 >  2:484409/1‑163 (MQ=255)
                                                                                       ggCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                 <  1:484409/163‑1 (MQ=255)
                                                                                          gTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCCCCTCTTCATCGCGCCGATCTTTGATGCGt                                                                                                                                                      >  1:32758/1‑251 (MQ=255)
                                                                                          gTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGAGTCAAAACGCCCAGCCGCAGGCATGCCAAGAGGATAGAATTGAGATCGGTGGCTTGCAAATCTTCCTGCCACAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTCAGCGTTGATAAACAATCGCGCGCTATAGCGGCGGCCCAGCGAAATTGCTGCCACGACATCGGCCTCCACCTCTTCATCGCGACGATCTTTGATGCgg                                                                                                                                                      >  1:32764/1‑250 (MQ=255)
                                                                                                ccAGCGCGACGTCGGTTTCCATATCGACATGTGTCGACACGCCCAGCCGCAGGGCTGACGCGCGGATGGCATTGAGATCGGAGGATTGCAAATCTTACTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATGGCGGCGGCACAGCGCAATTGCCGCCACGAAATCGGCTTCCACCACTTCATCGCGACGATCGTTGATGAGAAGCTgg                                                                                                                                                >  2:542477/1‑251 (MQ=255)
                                                                                                 cAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTgg                                                                                                                                                >  1:104535/1‑250 (MQ=255)
                                                                                                   gcgcGGCGTCGATTTCCATATCGTAAGGTGTCGAAACGCCCGGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgaga                                                                                                                                             <  2:410068/251‑1 (MQ=255)
                                                                                                   gcgcGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgaga                                                                                                                                             <  1:284772/251‑1 (MQ=255)
                                                                                                             gATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGTCCGc                                                                                                                                   >  2:212009/1‑251 (MQ=255)
                                                                                                                  ccATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTTCCGCCACGACATCGGCTTCCACCTCTTCCTCGCGCCGATCTTTGATGCGTACCTGGAGAGTACGTCCGCCTGcc                                                                                                                              >  2:433387/1‑250 (MQ=255)
                                                                                                                  ccATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCACCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGc                                                                                                                               >  1:86873/1‑250 (MQ=255)
                                                                                                                      atgggcatggggggAAACGCCCAGCATCAGGGCTGCCGCGCGGAGGGCATGGAGATCGGTGTCTTGAAAATCTTCCTGCCCAAAAGGGACGCCATCCGCCTGAGGC‑TGGATCGCCGGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCGTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCGGCATCCa                                                                                                                          <  2:516369/238‑1 (MQ=255)
                                                                                                                        cGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAAc                                                                                                                        <  1:176596/251‑1 (MQ=255)
                                                                                                                            ctgtgTTGAAAAGCCCATCCGCAAGCCTGCCGCGCCGATCGCATTGAGATTGGTGGCTTGCAAATCTTCCTGCCCCACACGGACGCCATACGCCTGATGC‑TTGATCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGAGCAATTGCCGACACGACATCGGCTTTCACCTCTTCATCGCGCCGATCTTTTCTGAGGAGCTGGAGAGTACGTATGCCTGCATTCAACAgcc                                                                                                                    >  2:375808/2‑249 (MQ=255)
                                                                                                                             ggggtggaaacgacaaggcGCAGGCCTGCGGCGAGGAGGGAATTGAGATGGGGGGCTGGCAAATCTCCCTGCCCAAAATGGACGCACTACGCATGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTGGAGAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCAGCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACg                                                                                                                   <  2:134723/233‑1 (MQ=255)
                                                                                                                              gtgccGAACCGCCTAGCCGCAGCCCGGCCGCGCGGATGGCTTTGAGATCGGTGGCGTGCAAATCGTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGt                                                                                                                  <  2:107686/247‑1 (MQ=255)
                                                                                                                                        ccccAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATCCGCCTGATGC‑TTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACt                                                                                                        <  2:234332/250‑1 (MQ=255)
                                                                                                                                              ccGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGc                                                                                                  >  1:301094/1‑251 (MQ=255)
                                                                                                                                                      cTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTCGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTccacc                                                                                          >  2:252437/1‑251 (MQ=255)
                                                                                                                                                           ccggggATGGCATTGAGAGCGGTGGCTTGCAAATCTTCCTGCCCCAAAGGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACAGCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACTTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCgg                                                                                     <  1:572272/247‑1 (MQ=255)
                                                                                                                                                              cgGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTa                                                                                  <  1:105855/251‑1 (MQ=255)
                                                                                                                                                                    gcttggagaTCGGTGCCTGGCAAATCTTCCGGCCCCAAATGGACGCCATCCGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTcc                                                                            <  2:40356/246‑1 (MQ=255)
                                                                                                                                                                    gCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACCGACGTTCGATCCACTGTACGCTTTCCACCCCCGGGTACAGTcc                                                                            >  2:47931/1‑251 (MQ=255)
                                                                                                                                                                     cATTGAGATCGGGGGCGTGCAAATCTCCCTGCCCCAAATGGACGCCAGACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCATCACCGGGTACAGTCCt                                                                           <  2:335654/251‑1 (MQ=255)
                                                                                                                                                                      aTTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTg                                                                          <  1:221555/251‑1 (MQ=255)
                                                                                                                                                                         gagaGCGGTGGCTTGCAAATCTTCCGGCCCCAAATGGACGCCATACGCCGGATGC‑TGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGaa                                                                        <  2:3470/250‑1 (MQ=255)
                                                                                                                                                                                gTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAgg                                                                >  2:181900/1‑251 (MQ=255)
                                                                                                                                                                                  ggCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTa                                                              <  1:244225/251‑1 (MQ=255)
                                                                                                                                                                                  ggCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTa                                                              <  1:310170/251‑1 (MQ=255)
                                                                                                                                                                                   gCTTGCAATTCTCCCTGCCCCAAATGGAAGCCGTACGCCTGATGC‑TTGATCGCCGGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTAc                                                             <  2:175357/251‑1 (MQ=255)
                                                                                                                                                                                        cAAATCTTCCGGCCCCAAATGGACGCCATACGCCTGAGGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAAGCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACaggag                                                        <  2:445395/251‑1 (MQ=255)
                                                                                                                                                                                            tCTTCCTGCACAAAAGGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGaa                                                     <  2:18788/250‑1 (MQ=255)
                                                                                                                                                                                               accTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGGAGCTGGGGAGTAAGTACGCCCGCATCCAACCGCCGTTCGAACCACTGTACGCTGTCCACCACCTGGGACAGTCCCGCACCAAAAGGGACCGGAGGAAAATc                                                 >  2:304826/2‑251 (MQ=255)
                                                                                                                                                                                                   gCCCCAAATGGACGCCATTCGCCTGATGC‑TTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCAGCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAgg                                              <  1:387989/250‑1 (MQ=255)
                                                                                                                                                                                                       cAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGAAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGAt                                         >  2:537319/1‑251 (MQ=255)
                                                                                                                                                                                                         aatgacccccATACGCCTGA‑GCGTGGAGCGCCAGCCGCCGGTAAGGGTTGAGAACCAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCATGCTGATAc                                       <  2:214504/243‑1 (MQ=255)
                                                                                                                                                                                                            ggACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATAGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTCCGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATc                                    >  2:270232/1‑251 (MQ=255)
                                                                                                                                                                                                                gccattccccTGATGC‑GTGGTCGCCAGCGGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGTCCCAGCGCAATTGCCGCCACGACATCGGCTTCCCCCTCTTCATCGCGCCGATCTTTGAGGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTGCGATCCACTGTTCGCTGTCCACCACCGGGTACAGTCCTGGACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGc                                <  1:237895/243‑1 (MQ=255)
                                                                                                                                                                                                                gCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGGACGCTGTCCACCACCGGGTACAGTCCTGAACGCAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGc                                >  1:128355/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCt                               >  2:197901/1‑251 (MQ=255)
                                                                                                                                                                                                                     aCGCCTGATGC‑TTGATCGCCAGCCGCCCGTACGCGCTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACCTCACGCTtcct                           >  2:355407/1‑251 (MQ=255)
                                                                                                                                                                                                                     aCGCCTGATGC‑TTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcc                            <  2:140066/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGCCTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcctc                          <  1:528789/251‑1 (MQ=255)
                                                                                                                                                                                                                         cTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCCACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGACAATCAGGCTGATACATCCCGCTTCCTCCtt                       >  2:168005/1‑251 (MQ=255)
                                                                                                                                                                                                                         cTGATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCtt                       >  1:9917/1‑251 (MQ=255)
                                                                                                                                                                                                                          ggATGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTa                      <  2:175781/250‑1 (MQ=255)
                                                                                                                                                                                                                            cgtactgtctCGCCAACCGCCATTTATCGTTGATAAACAATCGCGCGTTATAGCGGCCGCCCAGCGCAATTTCCGCCACGACATCGGCTTCCACCTCTTCATCGTGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACTCCTGCATCCAACAGACGTTCGATCCACTGTACGCGGTCCACCACCGGGTACATTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGTTACATCACGCTTCCTCCtt                       <  2:347143/240‑1 (MQ=255)
                                                                                                                                                                                                                            aTGC‑TTGATCGCCAGCCGCCAGTAAGCGTTGATAACCAATCGCGCGTTATACCGGCGGCCCAGCGCAATTGCCGCCACGACATCCGCTTCCACCTCTTCATCGCCACGACCTTTGATGCGTAGCTGGGGAGAACGTACGCCTGCATCCCAAAGACGCTCCAACCACTGTACGCTTTCCACCACCGGTTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAAGCTGGTACATCTCGCTTCCTCCttccg                    >  2:412498/1‑248 (MQ=255)
                                                                                                                                                                                                                                       gCCAGCCGCCAGTAAGCGGTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCTCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCAGCCAACAGGCGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAATGTACAGGAGTAAAATCAGGCTTATACATCACGCTTCCTCCTTACGCATGTAGAtt          <  1:574867/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ccgccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTcc     <  1:251767/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cgccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCg    <  2:9917/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGcc  >  1:444597/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGC‑TTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCC  >  NC_000913/4193751‑4194239

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