Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,216,473 A→G Y7H (TAC→CAC)  nfeF ← NADPH‑dependent ferric‑chelate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,216,4730AG98.9% 261.5 / NA 91Y7H (TAC→CAC) nfeFNADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (31/59);  minor base T (0/1);  total (31/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTT  >  NC_000913/3216226‑3216712
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
gTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCg                                                                                                                                                                                                                                              <  1:549350/251‑1 (MQ=255)
gTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCg                                                                                                                                                                                                                                              >  1:310191/1‑251 (MQ=255)
 tataGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCACTCTTCCGATACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGCAGAGTTTGCTGTGATCGGCAAAGCCACGCGATGGAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGGTGCGCTCAACGCGTAACACAGTCAGTTCCCGGAAGAGCAGATCATTGCGAACGCGCTGCGGGTGGCgg                                                                                                                                                                                                                                             >  2:555007/1‑251 (MQ=255)
  ataGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCgg                                                                                                                                                                                                                                             <  2:98986/250‑1 (MQ=255)
     gTCACGCGACGGTGGGCTTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCgggggg                                                                                                                                                                                                                                         <  1:578542/251‑1 (MQ=255)
       cactcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCAGTGCGAACGCGCTGTGGGTGGCGGGGGgt                                                                                                                                                                                                                                        <  2:214341/246‑1 (MQ=255)
         cgcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCGGATGCGCTCAACGCGTAACACAGTCAGTTCACGTAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGtt                                                                                                                                                                                                                                     <  2:573481/251‑1 (MQ=255)
          gcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGtt                                                                                                                                                                                                                                     >  1:488775/1‑250 (MQ=255)
           cgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTAt                                                                                                                                                                                                                                   <  2:188553/251‑1 (MQ=255)
             gcGGTGGGCGTGGTCCTTCCGGCCAGACGTTGCCCGCTTCCGTTACCGTTGGCTGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTTAACGCTTAACTCAGTCAGTTCACGGGAGCGCAGTTCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTc                                                                                                                                                                                                                                 <  1:261043/250‑1 (MQ=255)
             aCGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTc                                                                                                                                                                                                                                 <  2:133797/251‑1 (MQ=255)
               ggTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCAt                                                                                                                                                                                                                               >  1:461536/1‑251 (MQ=255)
                 ggggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGTGGGGTGTTATTCATTg                                                                                                                                                                                                                             <  1:140156/250‑1 (MQ=255)
                  gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGt                                                                                                                                                                                                                            <  1:7847/251‑1 (MQ=255)
                    gCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTa                                                                                                                                                                                                                          <  2:460494/251‑1 (MQ=255)
                      gTGGGCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGGGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTGCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGTTGTTATTCATTGTTATc                                                                                                                                                                                                                        <  2:577016/251‑1 (MQ=255)
                       tGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCg                                                                                                                                                                                                                       >  2:343938/1‑251 (MQ=255)
                        ggTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGc                                                                                                                                                                                                                      <  1:458519/251‑1 (MQ=255)
                            cTTCCGGCACGACGATGCCCTCTTCCGCTACCGATGGCGGCACCAAGAAAGCGTCAGGTTGAGGAAAGAAGCGTTTGCTGCAATCGTCAAACCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCCATGTTATTGCCTTc                                                                                                                                                                                                                  >  2:591411/1‑251 (MQ=255)
                              tCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTc                                                                                                                                                                                                                  <  2:449658/249‑1 (MQ=255)
                               ccGGCCCGACGATGCCCTCTTCCGTTACCGGTTGATGCACAAAGTGAGCGTCAGGTTGAGGAAAGACGAGTTTGCTTTGATCGTTAAAGCCACGCGATGTAAAACCGCCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAAGACAGTCAGTTCGCGGAAGCGCAGATCATTGCGGACGCGCTGCGGGTTGCGGGTGGTGTTATTCAtctgtatcccctgcgtg                                                                                                                                                                                                               <  2:19509/251‑17 (MQ=255)
                               ccGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTg                                                                                                                                                                                                               <  1:26421/251‑1 (MQ=255)
                                       aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                                                                                                                                                                       <  1:228752/251‑1 (MQ=255)
                                       aCGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAAt                                                                                                                                                                                                       <  2:585989/251‑1 (MQ=255)
                                        cGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATc                                                                                                                                                                                                      <  1:385839/251‑1 (MQ=255)
                                         gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGGAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCATTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGTGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCa                                                                                                                                                                                                     <  1:293608/251‑1 (MQ=255)
                                          aTGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAg                                                                                                                                                                                                    <  2:445045/251‑1 (MQ=255)
                                          aTGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAg                                                                                                                                                                                                    <  1:555017/251‑1 (MQ=255)
                                          aTGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAg                                                                                                                                                                                                    <  2:173750/251‑1 (MQ=255)
                                              ccTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatat                                                                                                                                                                                                <  1:20892/251‑1 (MQ=255)
                                              ccTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatat                                                                                                                                                                                                <  2:44864/251‑1 (MQ=255)
                                                   tCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATc                                                                                                                                                                                             <  1:146702/249‑1 (MQ=255)
                                                    ccGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCtaa                                                                                                                                                                                          <  1:606576/251‑1 (MQ=255)
                                                        tACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGGGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCGCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGTTTATTCATGGTTATCGCCTTCGGGATGGTAATCAGATATATCtaaata                                                                                                                                                                                       <  2:320053/250‑1 (MQ=255)
                                                         aCCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATAGCtaaataa                                                                                                                                                                                      <  1:145076/250‑1 (MQ=255)
                                                          ccGTTGGCGGCACAAAGGGAGGGGCAGGTTGAGGAAAGAAGAGTTTGCTGTGAGCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAGACGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGTGTGATATTCATTGTTAGCGCCTTGGTGATGTTAATCAGATATATCTAAATaaaa                                                                                                                                                                                    <  2:247129/251‑1 (MQ=255)
                                                             ttGGCGGCACAACGTGAGCGTCAGCTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTc                                                                                                                                                                                 >  2:378443/1‑251 (MQ=255)
                                                                             aGCGTCAGTTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGACGCTGAGGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCACAGATCAATGCTAACGCGCTGCTGGTGGCGATGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGAtt                                                                                                                                                                 <  1:405652/251‑1 (MQ=255)
                                                                                gTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGtt                                                                                                                                                              >  2:229457/1‑251 (MQ=255)
                                                                                gTCAGGTGGAGGAAAGAAGAGTTTGCTGGGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCGGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGGTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTGATTCATTGTTAGCTCCGGCGTGATGGTGATCGGATATATCTAAATAAAACTCGCAAACGAGAATGATTGtt                                                                                                                                                              <  2:535030/251‑1 (MQ=255)
                                                                                     gTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTACAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                                                                                                                                          >  1:367172/1‑250 (MQ=255)
                                                                                       tGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATACAACTCGCAAATGATAATGATTGTTAATCATg                                                                                                                                                       >  2:233804/1‑251 (MQ=255)
                                                                                        gAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGa                                                                                                                                                      <  1:431057/251‑1 (MQ=255)
                                                                                          ggAAAGAAGAGGTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAAGAAAACGCGCAAATGATAATGATTGTTAATCATGATa                                                                                                                                                    <  2:81574/251‑1 (MQ=255)
                                                                                           gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATa                                                                                                                                                    >  1:497660/1‑250 (MQ=255)
                                                                                               gaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGGCCAGCGCCTCGCCGCCGAGGACAATGCGCGGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCGGCGGGTGGTGGGGGGTGGTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAAGGATAATTATTGTTAGTCATGATAAATGc                                                                                                                                               <  2:585011/251‑1 (MQ=255)
                                                                                                 agagTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCaa                                                                                                                                             <  1:343938/251‑1 (MQ=255)
                                                                                                 agacgttGCTTTGAGCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCGCGCAGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCATATCGTTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCTTGGTTATCGCCTTCGGGATGGTTATCAGAGATAACAAAATAAAACTCGCTAAAGATAAGGATTGtgagttatgattaatgaaa                                                                                                                                             <  2:367237/246‑19 (MQ=255)
                                                                                                  gagTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAg                                                                                                                                            <  1:425557/251‑1 (MQ=255)
                                                                                                              aTCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGaa                                                                                                                                <  1:219665/251‑1 (MQ=255)
                                                                                                                cGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAAGACACTACGCTGAAAACCGGCGCTGATGAGATCACCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGCTATTCATTGTTATCGCCCTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACt                                                                                                                              >  1:283454/1‑251 (MQ=255)
                                                                                                                 gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                                             >  1:341494/1‑251 (MQ=255)
                                                                                                                  tCAAAGCCACGCGATGAAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGAGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCACATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTg                                                                                                                             >  1:610657/1‑250 (MQ=255)
                                                                                                                      aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGCGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACCGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATg                                                                                                                        >  1:253715/1‑251 (MQ=255)
                                                                                                                      aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATg                                                                                                                        <  1:427447/251‑1 (MQ=255)
                                                                                                                          aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCt                                                                                                                     >  1:512536/1‑250 (MQ=255)
                                                                                                                          aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTca                                                                                                                     >  2:82352/1‑249 (MQ=255)
                                                                                                                            gcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCtata                                                                                                                  <  1:142113/251‑1 (MQ=255)
                                                                                                                              gATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGt                                                                                                                >  2:373904/1‑251 (MQ=255)
                                                                                                                              gATGTAAAACCGTCCAGCGCCGCGCCGACGAGGACAAGGCGCTGAAAATCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCTCGGAAGAGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATACATTGTTATCGCCTTCGTTATGGTAATTAGATATAGCTAAATAAAACTCGCAAATGATAATGGTTGTTAATCATGATAAATGCAAGCGATGTGTAGAACGGATATGTCTATAGt                                                                                                                <  1:372798/251‑1 (MQ=255)
                                                                                                                                    aaaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                          <  1:98798/251‑1 (MQ=255)
                                                                                                                                     aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATaa                                                                                                         <  2:196921/251‑1 (MQ=255)
                                                                                                                                        cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGAc                                                                                                      <  1:330671/251‑1 (MQ=255)
                                                                                                                                         gTCCAGCGCCTCGACGGCGTGGACGATGCGCTGCAAACCGGCGCTGATGCGCCCAACACGTAACACAGTCACATCCCGGAGGCGAAGATCATTGCGCACGCGCTGCGGGTGGCGGGGGGGGTTATTCAATGGTATCGACTTCGAGATGGTAATACGATATCTCTCAATAAACATCGCGAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAAAACTGATATGTCTATAGTATGATAAGACg                                                                                                     >  2:307452/1‑251 (MQ=255)
                                                                                                                                             aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGCCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGAGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAAcc                                                                                                 >  1:335146/1‑251 (MQ=255)
                                                                                                                                               cgcCTCGCCGCCGAGGACACTGCGCTGAAAAACGGCGCTGCTGCGCTCCACGAATAACACAGTCCGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGATATTCATTGTTATCGCCTTCGCGATGGTAATCAGATATATCTAACTAAACCTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATACTCTGATAAGACGAACCGc                                                                                               >  2:599526/1‑251 (MQ=255)
                                                                                                                                                gcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGCTATATCTAAATAAAACTCGCAAATGATAATGCTTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGcc                                                                                              >  2:87179/1‑251 (MQ=255)
                                                                                                                                                    cgccgcCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCctc                                                                                            <  1:338804/249‑1 (MQ=255)
                                                                                                                                                    cgccgcCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCtctt                                                                                          <  2:35039/251‑1 (MQ=255)
                                                                                                                                                    cgccgcCAAGGACAATGCGCTGAAAACCGGCGCTGCTGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGAGCAGCTCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGATATTCCTTGTTATCGCCTTCGTGATGGTCATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGCTTTGTAGAACTGATCTGTTTATAGTCTGATAAGACGAACCGCCtctt                                                                                          >  2:516282/1‑251 (MQ=255)
                                                                                                                                                         ccgAGGACAATGCGCGGAAACCCGGGGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCTGTGCGAACGCGCTGGGGGTGGCGGGGGGTTTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTTTTAATAATGATAAATGCAGGCGATTTGTAGAACGGATAGGTTTATAGTCTGATAAAACGATCCGCCTCTTCTCAg                                                                                     <  2:73145/251‑1 (MQ=255)
                                                                                                                                                                       cTGAAAACCGGCGGTGATGCGCTCAACGCGTAACCCAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCAGCGGGTGGCGGGGGGGGATATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCCAAATAAAACTCGCAAACGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACAGATATGTCTATAGTCTGATAAGACGACCCTACTCTTCTCAAGCATCGTTAccca                                                                        >  2:472464/1‑247 (MQ=255)
                                                                                                                                                                         gAAACCCGGCGCTGAGGCGCTCAACGCGTAACACAGTCATTTCAGGGAAGCGCAGATCATTGCGAACGCGCTGCGGGGGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTTTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAAc                                                                      <  1:412525/250‑1 (MQ=255)
                                                                                                                                                                           aaaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCAt                                                                                                                                                        <  2:43774/166‑1 (MQ=255)
                                                                                                                                                                           aaaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCAt                                                                                                                                                        >  1:43774/1‑166 (MQ=255)
                                                                                                                                                                            aaCCGGCGCTGATGCGCTCAACGCGGAACACAGTCAGTTCACGAAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGGGTTATTCATTGTTATCGCCTTCGCGCTGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCg                                                                  >  1:310753/1‑251 (MQ=255)
                                                                                                                                                                              ccGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCgg                                                                 >  1:148065/1‑250 (MQ=255)
                                                                                                                                                                                      gATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGt                                                        >  1:14894/1‑251 (MQ=255)
                                                                                                                                                                                                  gcgTAACACAGTCAGTTCACGGAAGCGCAGTGCATTGCGTACGCGCTGCGGGTGGCGGGGGGTGTTATTCAGTGTTAGCGCCTTCCTGATGGTAATCAGat                                                                                                                                                                                                  <  2:236942/101‑1 (MQ=255)
                                                                                                                                                                                                  gcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTAATTGTTATCGCCTTCGTGATGGTAATCAGat                                                                                                                                                                                                  >  1:236942/1‑101 (MQ=255)
                                                                                                                                                                                                   ggTAACACGGGCAGTTAAGGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGGGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAAGGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATcc                                           <  2:488775/250‑1 (MQ=255)
                                                                                                                                                                                                         acaGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCCGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCGGATATATCTAAATAAAACTCGAAAATGATAAGGATTGTTAATCATGATAATTGCAAGCGATTTGTAGAACTGATAGGTCTATGGTCTGATAAGATGTACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                     <  1:317570/251‑1 (MQ=255)
                                                                                                                                                                                                          cagtcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGGGGCGGGGGGTGTTATTCATTTTTATCGCCTTCGGGATGGAAATCAGATATATCAAAAAAAACCTCGCAAATGATAATGATTGTTAATCAATATAAATGCAAGCGATGTGTAGACCGGTTATGTCTAGAGTCGGATAAGACGAACCGCCTCTTCTCAGTCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                    <  1:378443/251‑1 (MQ=255)
                                                                                                                                                                                                              catttCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTGGGGGGGGGTGTTATTCATTGTTAGCGCTGTCGTGATGGTAATCAGATATATCTAAAAAAAACTGGCAAATGATAAGGATTGTTAATCGTGATAAATGCAAGCGATTGGTAGAACTGATGGTTTTATAGTCTGATAAGAGGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                <  1:273421/248‑1 (MQ=255)
                                                                                                                                                                                                                ggtcccggaagcgcgaatatttGCGAACGCGTTGGGGGTGGCGGGGGGTTTTATTCATTGTTCTCGCCTTCGTGATGGTAATCAGATATATCTAAAAAAAACTCGAAAATGAAAAGGATTGTAAATCATGAAAAATGCAAGCTATTTGTAGAACGGATATGTTTATAGTCGGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTga                              <  2:148065/231‑1 (MQ=255)
                                                                                                                                                                                                                            gcAGATCGTTGCGACCGCGCTCGGGTTGGCGGGGGTTGTTTTTCATTGTTACCGCTTTCGTGTTGTTAATCAGATATATCTAAAAAAACCTCGCAAATTATAATTATTGTTAACCATGATAAAGGCAAGAGTTTTGTAGAACGGATATGTCATTAGTCTGATAAGTCGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCTGCCTGCGTGTGAGATGAGTCCACg                  <  2:610657/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAGATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACCCGCAAATGATAATGATTGTTAATCATGATAAATCCAAGCGTTTTGTAGAACTGATATTTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGt                 <  1:237592/251‑1 (MQ=255)
                                                                                                                                                                                                                               gATCATTGCGAACGCGCTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTAAAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGttt               <  1:233804/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTCATTGCGACCGCGCTGCGGGGGGGGGGGGGTGTTATTCATTGTTATCGCTTTCGTGATGGAAATAAGATATATATAAATAAAACTCGAAAATGAAAATGATTGTTAATCATGAAAAAGGAAAGCGATTTGTAGACCTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTa              <  2:158424/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgcTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCCAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGAATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGAGATGAGTCCCCGTTTAAAAATCAttttt  >  2:429303/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gcTGCGGGTGGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                               >  1:11576/1‑220 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
GTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTT  >  NC_000913/3216226‑3216712

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