Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,564,385 C→T S102N (AGC→AAC)  yjiL ← putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,564,3850CT99.2% 428.8 / NA 124S102N (AGC→AAC) yjiLputative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (60/63);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 127 total aligned reads displayed.

TCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCAC  >  NC_000913/4564138‑4564623
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
tCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTtgttgt                                                                                                                                                                                                                                               >  2:330126/1‑251 (MQ=255)
cccTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTtgttgt                                                                                                                                                                                                                                               <  1:292036/250‑1 (MQ=255)
cccTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTtgttgt                                                                                                                                                                                                                                               <  1:163057/250‑1 (MQ=255)
     cGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGa                                                                                                                                                                                                                                          <  2:152496/251‑1 (MQ=255)
          aTCGCTTCTGGCGCGACGCCCGCGGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGaccacc                                                                                                                                                                                                                                     <  2:268232/251‑1 (MQ=255)
             gCTTCCGGCGCGACGCCCGCTGAGCGCCGGCTGATCGCTCCTGATTCAGCAACCACTGTGCACATCCTCGGGATGGCGTGCGGCGTGACATTTTCGGAACTGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGAGCGCGAGATCACCTCCAGGAAACGCCCGGCGCCCGCCGCGCATTTGGCATTCATCAGGAAATCGCACAGGTTACCGGCATCATCCCGCTGAATCCCTTTGCtaccccgcccatcgct                                                                                                                                                                                                                                  >  2:55047/1‑236 (MQ=255)
                 cTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCg                                                                                                                                                                                                                              >  1:177040/1‑251 (MQ=255)
                   ggCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTc                                                                                                                                                                                                                               <  1:285132/248‑1 (MQ=255)
                    gcgcGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGAtt                                                                                                                                                                                                                           <  2:351559/251‑1 (MQ=255)
                    gcgcGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGAtt                                                                                                                                                                                                                           <  1:57515/251‑1 (MQ=255)
                        gACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACcg                                                                                                                                                                                                                       >  1:170520/1‑251 (MQ=255)
                         aCGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACcgc                                                                                                                                                                                                                      <  1:169627/251‑1 (MQ=255)
                                 tGAGCGCAGGCTGATAGATTCTGCTTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGCGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGACCGGTGCCAGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCCTCAAGCTGAATCACTTTGTTGTCATGCCCACCGATGTCGATTACCGCGCGCGTTg                                                                                                                                                                                                             >  2:445207/1‑251 (MQ=255)
                                       cAGGCTGGTCGCTTCTGATTCAGCAAACACTTGGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGgcg                                                                                                                                                                                                       <  2:392128/251‑1 (MQ=255)
                                          gCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCa                                                                                                                                                                                                     >  1:506564/1‑250 (MQ=255)
                                              aTCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGaa                                                                                                                                                                                                 >  2:372256/1‑250 (MQ=255)
                                                  cTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGt                                                                                                                                                                                            >  1:279807/1‑251 (MQ=255)
                                                       gATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTgcgc                                                                                                                                                                                        <  2:231622/250‑1 (MQ=255)
                                                       gATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTgcg                                                                                                                                                                                         <  1:132759/249‑1 (MQ=255)
                                                           cAGCAAACACTGTGCACATCCTCGTGATGGCGGGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAg                                                                                                                                                                                    <  2:185468/250‑1 (MQ=255)
                                                               aaaCACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGt                                                                                                                                                                               >  2:241134/1‑251 (MQ=255)
                                                                aaCACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTg                                                                                                                                                                              <  1:82728/251‑1 (MQ=255)
                                                                 acacTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTg                                                                                                                                                                              >  1:404857/1‑250 (MQ=255)
                                                                  cacTGTGCAAATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGAc                                                                                                                                                                            >  1:221072/1‑251 (MQ=255)
                                                                      gtgCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAgg                                                                                                                                                                         >  2:583863/1‑250 (MQ=255)
                                                                      gtgCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGACCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGCTTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAgg                                                                                                                                                                         >  2:627494/1‑250 (MQ=255)
                                                                              aCTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCg                                                                                                                                                                >  2:454992/1‑251 (MQ=255)
                                                                              aCTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCg                                                                                                                                                                <  1:201257/251‑1 (MQ=255)
                                                                                      tGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGtt                                                                                                                                                                                                              <  2:461663/197‑1 (MQ=255)
                                                                                      tGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGtt                                                                                                                                                                                                              >  1:461663/1‑197 (MQ=255)
                                                                                       ggCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGtt                                                                                                                                                       <  1:365189/251‑1 (MQ=255)
                                                                                          gggcggcgTGACATTGTCGGGAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTAt                                                                                                                                                    <  2:235524/249‑1 (MQ=255)
                                                                                           tGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATc                                                                                                                                                   >  2:490320/1‑251 (MQ=255)
                                                                                           tGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGGGACAGGAGATTTCCGTTACCTGTTTAtt                                                                                                                                                   >  2:521567/1‑250 (MQ=255)
                                                                                                 gTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCaaa                                                                                                                                             <  1:191604/251‑1 (MQ=255)
                                                                                                   gACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGCGATTTCCGTTACCTGTTTATCGGCAAAAt                                                                                                                                           >  2:735009/1‑251 (MQ=255)
                                                                                                           cGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGt                                                                                                                                   <  1:521567/251‑1 (MQ=255)
                                                                                                                  ggtggcGAGTGGCTCGACGCGGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATGTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTgccgcc                                                                                                                            <  2:445833/246‑1 (MQ=255)
                                                                                                                     gacGAGTGGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGGGCCCGCCGCGCATTTGTCAGTCATCAGGAAATCGCACAGGTGACCGTCATCATCAAGCGGAATCACTTTGTTGTCCTGACCACCGATGTCCATTACCTCGCGCGGTGCTGGCGCAAGAAACCGGGCGCCCAGGCCGTGACAGGAGATTTCCGTTACCTGTTGAGCGGCAAAATCCACCAGTTGGCGccc                                                                                                                           <  2:48192/247‑1 (MQ=255)
                                                                                                                     gTCGAGTTGGTCGACGCTGGTGCCAAGCGGGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGGTTATCGGAaaaagaaaaaaggggacgacc                                                                                                                           <  1:713909/249‑18 (MQ=255)
                                                                                                                     gTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGGCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCg                                                                                                                          <  1:650628/250‑1 (MQ=255)
                                                                                                                       cgggTTGCCCGACGCTGGTGCCAAGCGTGCGCGAGAGCACCTCCAGGAAACGCCCGGGGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCAGCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAg                                                                                                                       <  2:15501/248‑1 (MQ=255)
                                                                                                                        gAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGACCCGAACGGTGCCCGCCGAGCATTTGTCATTCATCAGGAACTCGCCCAGGTTACCGTCATCATCAAGCTCAATCACTTTGTTATCCTGCCCACCGATCTCGATTACCGCGCGCGTTGCTGGCGCAAGAAAACGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGc                                                                                                                      >  1:41561/1‑251 (MQ=255)
                                                                                                                        gAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGc                                                                                                                      >  1:565453/1‑251 (MQ=255)
                                                                                                                         ggTTGCTCGACGCTGGAGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCGCGGTTCCCGCCGCGCATTTGTCATTCAGCAGGAAATCGCACAGGTTACCGGCATCATCAAGCTGAACCACTTTGTTGTCCGGACCTCCGATGTCGATTACCGCGCGCGTTGCGTGCGCAAGAAGCCGTGCGCGCAGCCCGTGACAGGAGATTCCCGTTACCTGTGTATCGGCAAAATCCACCAGTTGCCGCCCGTAGcc                                                                                                                     <  1:100632/250‑1 (MQ=255)
                                                                                                                          gTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCg                                                                                                                    <  1:241134/251‑1 (MQ=255)
                                                                                                                          gTTGCGCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCg                                                                                                                    <  2:215573/251‑1 (MQ=255)
                                                                                                                                 gACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCg                                                                                                            >  1:377854/1‑251 (MQ=255)
                                                                                                                                    gCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCGCCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCa                                                                                                         <  2:177040/251‑1 (MQ=255)
                                                                                                                                       ggTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGaa                                                                                                      <  1:213797/251‑1 (MQ=255)
                                                                                                                                        gTGCCAAGCGTGCGCGAAATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTAATCAGGAAATCGCCCAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCAGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGaa                                                                                                      >  2:326683/1‑250 (MQ=255)
                                                                                                                                            cAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGc                                                                                                 >  1:239540/1‑251 (MQ=255)
                                                                                                                                             aaGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGCTGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACCGGCGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCg                                                                                                >  1:387545/1‑251 (MQ=255)
                                                                                                                                                 gTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCATTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGt                                                                                            >  2:273641/1‑251 (MQ=255)
                                                                                                                                                  tGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGtc                                                                                           >  2:377844/1‑251 (MQ=255)
                                                                                                                                                      cgAGATCACCTCCAGGACACGCCAAGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCTTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCTCGTAGCCGGTGAGCGTCCGAAACGGCGTTGTctct                                                                                        >  2:486844/1‑250 (MQ=255)
                                                                                                                                                            cACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCtt                                                                                 >  1:640671/1‑251 (MQ=255)
                                                                                                                                                             aCCTCCAGGCAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTc                                                                                >  2:251310/1‑251 (MQ=255)
                                                                                                                                                              ccTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTGGCTGGCGCAAGAAACCGTGCGCCCAGCTCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTcg                                                                               <  1:445207/251‑1 (MQ=255)
                                                                                                                                                                tCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTcgcg                                                                             >  1:201505/1‑251 (MQ=255)
                                                                                                                                                                 ccAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAACCCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCgcgc                                                                            <  2:565453/251‑1 (MQ=255)
                                                                                                                                                                    ggAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATAACCGGCGGCGTGGCTGCCGCAGGAACCGGGGCGCCCAGCCGGGGCCAGGAGTTTTCGGTTACCGGTTATTCGGCAAAACCCACCAGTTCCCCCCCGAAGCGGTGGGGCGTCGGAAACGGGGTTGTCTCTACCCCtcgggccggg                                                                         >  2:165976/1‑242 (MQ=255)
                                                                                                                                                                     gAAACGCCCGGTGCCCGCCGCGCATTTGTCACTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCagag                                                                        <  1:273641/251‑1 (MQ=255)
                                                                                                                                                                      aaaCGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGt                                                                       >  1:336705/1‑251 (MQ=255)
                                                                                                                                                                          gCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTc                                                                    >  2:579939/1‑250 (MQ=255)
                                                                                                                                                                           cccGGTGCCCGCCGCGCATTTGCGATTCATCAGGCAATCGCGCAGGTTACCGTCATCATCAGGCTGAATCCCTTTGTTGTCCTGACCACCGCTGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGACCAGCCCGTGACAGGAGCTTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTCGCCGGTGAGCGTCAGAAACGCCGTTGTCTCTAACCCTTCGCGCAGAGTTTccc                                                                  >  2:652482/1‑251 (MQ=255)
                                                                                                                                                                            ccGGGGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCa                                                                 <  2:22735/251‑1 (MQ=255)
                                                                                                                                                                              ggTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAgg                                                               >  2:272478/1‑251 (MQ=255)
                                                                                                                                                                              ggTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAgg                                                               <  1:520556/251‑1 (MQ=255)
                                                                                                                                                                               gTGCCAGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGc                                                              >  2:652752/1‑251 (MQ=255)
                                                                                                                                                                                 gCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCtt                                                            <  1:503799/251‑1 (MQ=255)
                                                                                                                                                                                 gCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCtt                                                            <  2:263998/251‑1 (MQ=255)
                                                                                                                                                                                 gCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCt                                                             <  1:432984/250‑1 (MQ=255)
                                                                                                                                                                                   cagcctcgcATTTGGCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATGACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTGATCGGCAAAATCCACCAGTTGTCGCCCGTAGCCGTTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTGCCCAGGCTTCa                                                          <  2:40192/245‑1 (MQ=255)
                                                                                                                                                                                    cgccgcGCATTTGACATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCACAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCAACAGTTGCCGCCCGTAGCCGGTAAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTACCAGGCTTc                                                           >  1:711638/1‑249 (MQ=255)
                                                                                                                                                                                     gccgcGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGt                                                        <  1:583422/251‑1 (MQ=255)
                                                                                                                                                                                      ccgcGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGt                                                        >  1:537776/1‑250 (MQ=255)
                                                                                                                                                                                      ccgcGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCGTAATCACTTTGTTTTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTa                                                       <  2:576319/251‑1 (MQ=255)
                                                                                                                                                                                            aTTTGTCATTCATCAGGAAATCGCCCAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCt                                                 >  2:71627/1‑251 (MQ=255)
                                                                                                                                                                                              ttGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTg                                                <  2:286891/250‑1 (MQ=255)
                                                                                                                                                                                              gtGTCATTCATACGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTGGGTGGCGCAAGAAACCGGGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCAGGAGAAGGTGAGCGGCAGAAACGGCGTTGTCTCTAACTCTTCGCGCAGAGTTTCCAAGGCTTCAGTAATTGCTg                                                <  2:325141/249‑1 (MQ=255)
                                                                                                                                                                                                     tcatcaGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCgg                                         <  1:455785/250‑1 (MQ=255)
                                                                                                                                                                                                         cagaaaaTCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATTTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACCGGAGATTTCCGTTACCTGTTTCTCGACACAATCCGCCAGTTGCCGCCCGTATCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCCGGGTTTCCAAAGCTTCAGTAATTGCTTTTTCCGGGCgcc                                    >  2:389499/5‑249 (MQ=255)
                                                                                                                                                                                                             aaaTCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAACTCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGCGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAggg                                >  2:242/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaTCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAggg                                <  1:258945/251‑1 (MQ=255)
                                                                                                                                                                                                              aaTCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGTCAAAAGCCACCAGTTGCCGCCCGTAGCCTGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTGCAGTAATTGCTGTTGCCTGGCGAAAgggg                               <  2:721985/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGt                             <  1:678029/250‑1 (MQ=255)
                                                                                                                                                                                                                  gCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGACAGGGGGTTg                           >  1:657490/1‑251 (MQ=255)
                                                                                                                                                                                                                    acaGGTAACAGTCATCAACAAGCAGAATCCCTTTGTTGTCCTGACCACCGATGGCGATTACCGCGCGCGTGGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTGCCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTGCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGa                         <  2:380149/251‑1 (MQ=255)
                                                                                                                                                                                                                      aGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAAc                       <  1:13168/251‑1 (MQ=255)
                                                                                                                                                                                                                           aCCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAgg                   >  1:506977/1‑250 (MQ=255)
                                                                                                                                                                                                                            acGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCAGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGCGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGCACGGGGGTTGGAACGAGGa                  >  2:520194/2‑250 (MQ=255)
                                                                                                                                                                                                                             cGTCATCATCACGGTGAATCGCTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTGCCGTTACCTGTTTATCGGCAAAATCCCCCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGaaa                <  2:699966/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGaa                 >  1:5829/1‑250 (MQ=255)
                                                                                                                                                                                                                             cGTCAACATCAGTCTGAATCTCTTTGTTGGCCTGACCACCGATGGCGAGTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTAGGCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGTCTTGGAACGAGGaaa                <  2:475918/251‑1 (MQ=255)
                                                                                                                                                                                                                                      caAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGt        <  1:372256/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCACAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCCGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATc    >  2:236242/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGCGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATc    >  2:197294/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCAc  <  1:490320/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCAc  <  1:710278/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
TCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCAC  >  NC_000913/4564138‑4564623

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