Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 83,182 G→A D116D (GAC→GAT leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,1820GA92.7% 237.6 / 7.8 83D116D (GAC→GATleuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (31/45);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

GCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTC  >  NC_000913/82940‑83424
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
gCTTTGTCGATGTTAGGTCCGCGTTCATACAGGCCGCTGATGATTCCGGCGCACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGAGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCCGATCGGCAATGGGTGTACGCCCGGCATCTTCGCACGAAACTTCAACATCATCGGTGTCATTACGGGCGCGTTTCACCATCTAGATAGCGCGTTCGATCACTTCATCCAGCGt                                                                                                                                                                                                                                            >  2:69574/1‑251 (MQ=255)
 cTTTGTCGATGTTAGGTCCGCGTCCATACAGGCCGCTGATGATTCCGGCCAACTCAACCGGAATGCGGTAGCCCACGGCGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCATCGCCCACTCGCGCCAGATCGGCAATGGGTACACGCCCGGCTTCTTCGCCAGAAAATTCAACATCATCGGTGTAGTCACGGGCGAGTTTCACCATATCGCAAGCACGTTCGATCACGTACTCCAGCGtg                                                                                                                                                                                                                                           >  2:601103/1‑251 (MQ=255)
 cTTTGTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCAATTCATCCAGCGtg                                                                                                                                                                                                                                           >  2:463404/1‑251 (MQ=255)
     gTCGATGTTAGGTACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGtgctgc                                                                                                                                                                                                                                       <  1:297145/251‑1 (MQ=255)
                gtaCGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGCGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGCGCTGCGCAGCTtggtg                                                                                                                                                                                                                            >  2:502005/3‑251 (MQ=255)
                gtaCGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGACCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTtggtg                                                                                                                                                                                                                            <  2:348951/249‑1 (MQ=255)
                 cACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTCATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTtggtgg                                                                                                                                                                                                                           >  2:125975/1‑251 (MQ=255)
                  aCGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGc                                                                                                                                                                                                                          <  1:692890/251‑1 (MQ=255)
                    gcgTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCgg                                                                                                                                                                                                                        >  1:574614/1‑250 (MQ=255)
                     cgTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGAt                                                                                                                                                                                                                       >  1:582559/1‑251 (MQ=255)
                         cATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGAtgtg                                                                                                                                                                                                                    <  1:554002/250‑1 (MQ=255)
                          aTACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGCTGTGc                                                                                                                                                                                                                   >  1:14704/1‑250 (MQ=255)
                               ggCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGc                                                                                                                                                                                                             <  1:463404/251‑1 (MQ=255)
                                gCCGCTGATGAGTCCGGCGAACTCAACCGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCAACGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCg                                                                                                                                                                                                            <  2:327290/251‑1 (MQ=255)
                                 ccGCGGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACGTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGa                                                                                                                                                                                                           <  2:592660/251‑1 (MQ=255)
                                       atgatTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGc                                                                                                                                                                                                     >  1:225277/1‑251 (MQ=255)
                                                  caaCTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACCCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCTCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTAt                                                                                                                                                                                          <  2:662928/250‑1 (MQ=255)
                                                   aaCTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATg                                                                                                                                                                                         <  1:224543/251‑1 (MQ=255)
                                                    aCTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAAGCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTTCGCCCGGCATCTGCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCGTGGTGGCGATGTGCATTGGCGAAGTTGCAATAAAGGTATGa                                                                                                                                                                                        <  2:37263/251‑1 (MQ=255)
                                                        aaaCGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGAGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATAc                                                                                                                                                                                    <  1:502005/251‑1 (MQ=255)
                                                          aCGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACgg                                                                                                                                                                                  <  1:327453/251‑1 (MQ=255)
                                                            ggCATGGTGTAGCCCACGGTGCCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGaa                                                                                                                                                                                <  2:239218/251‑1 (MQ=255)
                                                                       gCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGa                                                                                                                                                                     <  1:63653/251‑1 (MQ=255)
                                                                           ccggggtCCGGACTGTTGAGGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCTCAAGAAAAACCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCAGTGGCGAAGTGGCAATAAAGGTATGAATCCGGAAGGCTTCGGCGACttg                                                                                                                                                                 <  2:108774/246‑2 (MQ=255)
                                                                            cGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCGTCGGCGACTTTc                                                                                                                                                                <  1:212198/251‑1 (MQ=255)
                                                                                     gAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCg                                                                                                                                                       <  1:527731/251‑1 (MQ=255)
                                                                                     gAATGTTGAGGGTGGTGGCATCGGCATTAACCTCCGCTTCGACCACTCGCTCCAGATCGGCAATGGGTGTACGCTCGGCATCTTCGCAAGAAAATCCAACATCATCGGTGTAATTACGGGCTCGTTTCACCATATAGACAGCGCTTTCGATCACTTCATCCAGCGTGCTGCTCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCTATAAAGGTATGAATACGGAAGGCTTCGTCGACTTTCAGGGAtt                                                                                                                                                         <  2:702109/249‑1 (MQ=255)
                                                                                      aaTGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCgg                                                                                                                                                      <  1:259906/251‑1 (MQ=255)
                                                                                                         ccGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCttt                                                                                                                                   >  1:460180/1‑251 (MQ=255)
                                                                                                          cGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGGCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCtttc                                                                                                                                  <  1:688845/251‑2 (MQ=255)
                                                                                                            gCATTAATCGCCGCTTCGACCACTCGCGCCATATCGGCAATGGGTGTACTCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCTCCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGTTGCGCAGCTTGGTGGCGATTTGCATGGGTGAAGTGGTAATAAAGGTATGAATACGGAAGGCTTCTGCGACTGTCAGGGATTCGGCCGCAACGTCGATATCtctctc                                                                                                                                <  2:13109/251‑6 (MQ=255)
                                                                                                            gCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCtttctc                                                                                                                                >  1:507375/1‑248 (MQ=255)
                                                                                                             cATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTcc                                                                                                                               <  1:18316/251‑1 (MQ=255)
                                                                                                                tAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGGATGAATACGGAAGGCTACGGCGACTTTCAGTGATTCGGCCGCCACGTCGATATCtttatcaacg                                                                                                                            >  2:301160/1‑244 (MQ=255)
                                                                                                                  aTCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCa                                                                                                                          <  1:274896/251‑1 (MQ=255)
                                                                                                                       cgcTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAg                                                                                                                     <  1:174577/251‑1 (MQ=255)
                                                                                                                       cacTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACCTCATCGGTGTCATTACGGGCGCGTTTCCCCATATAGATAGCGCGTTCGCTCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGCATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTATACACGCAGCggg                                                                                                                     >  2:639096/3‑249 (MQ=255)
                                                                                                                           gcccccACTCGCGCCAGATTGGTAAATCGTTTACGCCCGGCATCTTCGCAAGAAAATTCAACAGCATCGGGTAAAGTCAGGGCGCGTTGCACCATAGAGATATCGCGTTCGATCACTCCATCCAGCGTGCTGCGCAGCTTGGGGGCGATGTGCAGTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGAATCGTCCGCCACGTCGATAGCTTTCTCCACGCATCGAGCTaa                                                                                                                 <  2:484214/247‑1 (MQ=255)
                                                                                                                             gtccACGCGCGCCAGAGCGGCAATGGGTGTACGCCCGGCATCGTCGCAAGAAAATTCAACATCATCGGTGGAATTACGGTCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACg                                                                                                               <  2:39286/249‑1 (MQ=255)
                                                                                                                             gACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACg                                                                                                               >  1:632180/1‑251 (MQ=255)
                                                                                                                              aCCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGc                                                                                                              >  1:93486/1‑251 (MQ=255)
                                                                                                                               cccagcgcgCCATATCGGCAATGGTGGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGTGCGTCCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAGGTGGCAATAAAGGTGTGAATCCGGAAGGCTTCGGCGCCTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGc                                                                                                              <  2:168196/245‑1 (MQ=255)
                                                                                                                                      cgcCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATAcgc                                                                                                      >  1:77596/1‑251 (MQ=255)
                                                                                                                                       gcaagaccGGCAATGTGTGTACGCCCGGCAATTTCGCAAGAAAACTCAATATCATCGGGGTATTGTCGGGCGAGTTTCACCATATAGATAGCGCGTTAGATCACTTCAACCAGCGTGTGGCGCAGCTTGGTGGCGAGGTGTATTGGCGAAGTGGCAATAAAGGTATGAATATGGAAGGCTTCGGCGCCCTTCAGGGAGTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACTCGCATAcgcg                                                                                                     <  2:420152/244‑1 (MQ=255)
                                                                                                                                        ccAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCACGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGACGTGGCACTAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATAcgc                                                                                                      >  1:307604/1‑249 (MQ=255)
                                                                                                                                        ccAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCgg                                                                                                    <  1:572364/251‑1 (MQ=255)
                                                                                                                                            aTCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGt                                                                                                <  1:125975/251‑1 (MQ=255)
                                                                                                                                              cGGCAATGGGTGTACGCCCTGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCCCCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGttt                                                                                              <  1:260883/251‑1 (MQ=255)
                                                                                                                                              cGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACGTCATCCAGCGGTCTGCGCAGCTTGGTGGCGATGTGCATTTGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCTACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGAATACGCGGCTGtt                                                                                               <  2:41399/250‑1 (MQ=255)
                                                                                                                                               ggCAAGGGGTGTTCGTCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTGCACCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGtttt                                                                                             <  2:96507/251‑1 (MQ=255)
                                                                                                                                                   aTGGGTTTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTTATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGTGAGCTAACGCGCATACGCGGCTGTTTTTAAc                                                                                         <  2:562203/251‑1 (MQ=255)
                                                                                                                                                      ggggtACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATCCGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGCTATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTg                                                                                      >  2:548360/4‑251 (MQ=255)
                                                                                                                                                         gtACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAACGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGAGCATTGGCGAAGTGGCAATAAAGGGATGAATACGGAAGGCTTCGGCAACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCTGCGAGCTACCGCGCATACACGGCTGTTTTTAACCTggcg                                                                                   >  2:153570/1‑251 (MQ=255)
                                                                                                                                                         gtACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGGAATTACGGGCGCGTTTCACCATATAGATAGCGCGTGCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTggcg                                                                                   <  2:477979/251‑1 (MQ=255)
                                                                                                                                                             gCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTggcgagcg                                                                               >  1:383829/1‑251 (MQ=255)
                                                                                                                                                               ccGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGAt                                                                             >  2:40821/1‑251 (MQ=255)
                                                                                                                                                                cGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATg                                                                            <  1:479724/251‑1 (MQ=255)
                                                                                                                                                                   cATCTTCGCAAGAAAATTCAACACACACGGTGTAATTACGGGCGCGTTTCACCATAGAGATAGCGCGTTCGAGCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATTTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGGCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGtt                                                                         <  2:632180/251‑1 (MQ=255)
                                                                                                                                                                      cTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAACGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCCTACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGc                                                                      >  2:343315/1‑251 (MQ=255)
                                                                                                                                                                        tCGCAAGAAAATTCAACATCATCGTTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTGTGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATTCGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTg                                                                       <  1:443677/248‑1 (MQ=255)
                                                                                                                                                                        tCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGttt                                                                        <  1:117590/247‑1 (MQ=255)
                                                                                                                                                                          gCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCg                                                                  >  1:149244/1‑251 (MQ=255)
                                                                                                                                                                            aaGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGAt                                                                >  1:221789/1‑251 (MQ=255)
                                                                                                                                                                            aaGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGAt                                                                >  1:61402/1‑251 (MQ=255)
                                                                                                                                                                                  aTTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCaaa                                                           >  2:145222/1‑250 (MQ=255)
                                                                                                                                                                                       aCATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGCTATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATcacc                                                     >  2:185919/1‑248 (MQ=255)
                                                                                                                                                                                                         tACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACggg                                      <  1:349086/250‑1 (MQ=255)
                                                                                                                                                                                                          aCGGGCGCTTTTCACCATATAGATAGCGCGTTCGATCACGTCACACAGCGTGCTGCGCAGCTTGGGGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGATTACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgggg                                     <  2:77596/250‑1 (MQ=255)
                                                                                                                                                                                                          aCGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgggg                                     <  1:81305/250‑1 (MQ=255)
                                                                                                                                                                                                          aCGGCCGCGTTTCACCATATAGATAGCGCTTTCGATAACTTCATCCAGCGGGCGGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATCCGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCTAGCTAACGCGCATACTCGGCTGTTTTTAACCTGGCGATCGATGGTTTGCACTGATTCAAAATCACCCGTCGAAGAAACGGGGa                                    <  2:647212/251‑1 (MQ=255)
                                                                                                                                                                                                                   tttCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGGGGCGATGTGCATTGGCGAAGTGGCAATAATGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGGTTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACtt                           >  1:55223/1‑251 (MQ=255)
                                                                                                                                                                                                                      cACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTcc                         <  1:185919/250‑1 (MQ=255)
                                                                                                                                                                                                                        ccATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCAc                                                                                                                             >  2:45644/1‑148 (MQ=255)
                                                                                                                                                                                                                        ccATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCAc                                                                                                                             <  1:45644/148‑1 (MQ=255)
                                                                                                                                                                                                                              agatagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTc                 >  1:447966/1‑250 (MQ=255)
                                                                                                                                                                                                                               gatagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCCAAATCACCCGGGGAAGCAACGGGGAAACCGCCTTCCATCACGTCa                >  2:21517/1‑250 (MQ=255)
                                                                                                                                                                                                                                atagCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAac              >  1:241202/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgcgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACAccc          >  1:185003/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgcgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGAGGCGTTGTGAATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTACCGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACAccc          >  1:557466/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCAt        >  1:337034/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATa       <  1:604572/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCAGTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATAc      <  1:452708/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTc  >  1:366018/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
GCTTTGTCGATGTTAGGCACGCGTTCATACAGGCCGCTGATGATTCCGGCGAACTCAAACGGCATGGTGTAGCCCACGGTGTCCGGAATGTTGATGGTGGTGGCACCGGCATTAATCGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTC  >  NC_000913/82940‑83424

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