Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,230 T→A 5.9% intergenic (+36/+193) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2300TA 5.9% 568.7 / 17.4 171intergenic (+36/+193)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (85/72);  new base A (5/5);  total (91/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 174 total aligned reads displayed.

AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGC  >  NC_000913/3422985‑3423502
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                 
                                                 tgttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCACCCCCCGGGCCGTTTCtt                                                                                                                                                                                                    >  4:520817/1‑275 (MQ=255)
                                                  gttgtGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAACCTCCGGGCCGTTTCtt                                                                                                                                                                                                    >  2:337228/1‑274 (MQ=255)
                                                      tgtgcAGCACATGAGATGGGGTTGGCACGAACCGTCGCTGACCGGGTAATTTTGATTGATCGTGTGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTGTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                <  1:1432/274‑1 (MQ=255)
                                                             cacaTGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCaa                                                                                                                                                                                        >  2:558981/1‑275 (MQ=255)
                                                                       ggggTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                                                              >  1:87900/1‑275 (MQ=255)
                                                                            ttGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCCCCCCcggcgg                                                                                                                                                                         >  2:594224/1‑275 (MQ=255)
                                                                             tGCACGAACCGTCGCTGACTGGGTAATTTTTAGGGATCGCGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTATTGCGCATCCTAAATCAGCGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATAGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTCGATGTCTGGCAGTCAATGGCGGGGGTCCTGCCCGCCACCCCACGGGCCGGTGCTTGGCTACGGTGAAACCCGCCACcggacgg                                                                                                                                                                        >  4:593825/1‑271 (MQ=255)
                                                                                       gtccgtGACAGGTTAATTTTTAGGGTTGAGGGGGAAAATGTTGAGCAAACGGACCCTGATGAATTTTTTTCGCATCCTAAATCATAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGCTTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                                                               <  4:294022/269‑1 (MQ=25)
                                                                                         tgCTGACCGGGGAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                                                                                                                                             <  2:418854/273‑1 (MQ=255)
                                                                                             gcacGGGTAATTTTTATGGTTCGTGGGGAAATAGTGGGGCAACCTGCAACTGATGAATTTTTTGCGCATCCTAAATCAGAGGGTACGGGGGCATTTTTAGCGCAGGTAATCCATTAATTAAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg                                                                                                                                                         <  4:558981/271‑1 (MQ=255)
                                                                                              aCCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCataaactgccag                                                                                                                                                                                                           >  2:486372/1‑212 (MQ=255)
                                                                                                      aTTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGa                                                                                                                                                                                                                                                                              <  6:536943/149‑1 (MQ=255)
                                                                                                          ttATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          <  8:835006/149‑1 (MQ=255)
                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCTTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCCGTTTATGGCGGGCGTCCTTCCCTCCACCCTCCGGGCCGTTGCTTTGCAACGTTCAAAGCCGCTCCCCGCGGGTTTGTCCTCCTCGCTCAAGTGTTCACCGaaca                                                                                                                                        >  4:593970/1‑272 (MQ=255)
                                                                                                             tGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCtt                                                                                                                                                                                                                                                                      <  5:305452/149‑1 (MQ=255)
                                                                                                               gATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTc                                                                                                                                                                                                                                                                    <  8:418133/149‑1 (MQ=255)
                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGc                                                                                                                                                                                                                                                                  <  6:763247/149‑1 (MQ=255)
                                                                                                                 tCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGc                                                                                                                                                                                                                                                                  >  8:126180/1‑149 (MQ=255)
                                                                                                                  cGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATAGCCTTTCGCt                                                                                                                                                                                                                                                                 <  8:311236/149‑1 (MQ=255)
                                                                                                                   agagacggaTAGTGGAGCAAGCTGAACCTTAGTGATTTTTTTGGAACCCTAAATAAGAGGGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAGAAGGCCATCATTTCGGATGGCCTTTCGCTTGATTTGAGGTCTGGCAGTTTATGGCGGGCGTCCTGTCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTAACCGACAaacaac                                                                                                                                   <  4:283140/266‑1 (MQ=21)
                                                                                                                    tGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttg                                                                                                                                                                                                                                                               >  7:591964/1‑149 (MQ=255)
                                                                                                                          aaTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttga                                                                                                                                                                                                                                                         <  5:572356/149‑1 (MQ=255)
                                                                                                                           aTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                                                                                                                                                                                        <  7:633928/149‑1 (MQ=255)
                                                                                                                           aTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttgat                                                                                                                                                                                                                                                        <  8:362829/149‑1 (MQ=255)
                                                                                                                              gTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGGGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCttt                                                                                                                                                                                                                                                                     <  7:755456/133‑1 (MQ=255)
                                                                                                                              gTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCttt                                                                                                                                                                                                                                                                     >  8:755456/1‑133 (MQ=255)
                                                                                                                                ggAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTg                                                                                                                                                                                                                                                   <  5:282269/149‑1 (MQ=255)
                                                                                                                                  aGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGc                                                                                                                                                                                                                                                 >  5:1003534/1‑149 (MQ=255)
                                                                                                                                   gCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGGGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                                                                <  4:305626/197‑1 (MQ=255)
                                                                                                                                   gCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          >  7:106151/1‑124 (MQ=255)
                                                                                                                                    cAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTACTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTTAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACCGATAACACAAAAGAcca                                                                                                                 >  2:145361/1‑274 (MQ=255)
                                                                                                                                          gCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATg                                                                                                                                                                                                                                         <  7:908095/149‑1 (MQ=255)
                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                     <  6:1003534/149‑1 (MQ=255)
                                                                                                                                              cTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgg                                                                                                                                                                                                                                     >  6:418727/1‑149 (MQ=255)
                                                                                                                                               tgatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                                                                                       <  2:520817/274‑1 (MQ=255)
                                                                                                                                                gatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                   >  7:999221/1‑149 (MQ=255)
                                                                                                                                                gatgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggc                                                                                                                                                                                                                                   >  8:591272/1‑149 (MQ=255)
                                                                                                                                                    gaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttga                                                                                                                                                                                                                                                         >  6:981818/1‑124 (MQ=255)
                                                                                                                                                    gaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttga                                                                                                                                                                                                                                                         <  5:981818/124‑1 (MQ=255)
                                                                                                                                                     aaTTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGc                                                                                                                                                                                                                                                                          <  5:132601/107‑1 (MQ=255)
                                                                                                                                                                 aTCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCAccc                                                                                                                                                                                                                   >  6:716192/1‑149 (MQ=255)
                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTGGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                                                                                                                                                >  5:1062505/1‑149 (MQ=255)
                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGt                                                                                                                                                                                                                                                      >  7:203006/1‑111 (MQ=255)
                                                                                                                                                                    cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGt                                                                                                                                                                                                                                                      <  8:203006/111‑1 (MQ=255)
                                                                                                                                                                     tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                                                                                                                                               >  5:832907/1‑149 (MQ=255)
                                                                                                                                                                        aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCGGGcc                                                                                                                                                                                                           <  5:418727/149‑1 (MQ=255)
                                                                                                                                                                            gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                        <  7:162410/149‑1 (MQ=255)
                                                                                                                                                                            gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                                                                        <  8:440588/149‑1 (MQ=255)
                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaaa                                                                                                                            >  4:194432/1‑224 (MQ=255)
                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                          >  2:598184/1‑274 (MQ=255)
                                                                                                                                                                             agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCAACTGAGCCTTTTGTTTTATTTGATTTCTgg                                                                          >  2:273027/1‑274 (MQ=255)
                                                                                                                                                                              gCGTACGGGGGTTTTTTTTTCGCAGGTAATCCATTAATTTAATTGTAGTTCGAAAAGCAAAAAGGTCATCCGTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                        <  4:145361/275‑1 (MQ=255)
                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATGCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  5:716192/149‑1 (MQ=255)
                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGACAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  8:562208/149‑1 (MQ=255)
                                                                                                                                                                              gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                                                                      <  6:135587/149‑1 (MQ=255)
                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCtt                                                                                                                                                                                                                                                                      >  1:5901/1‑82 (MQ=255)
                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                                                                                                                                                      <  1:126340/194‑1 (MQ=255)
                                                                                                                                                                                 tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaa                                                                                                                        <  2:194432/224‑1 (MQ=255)
                                                                                                                                                                                    ggtggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTCGTTCGAAAAGCAAAAAGGCCCTCCTGTCGTGTGGCCTTTCCCTTGATTTGATGTCTGTCAGTTTATGGCGGGCGTCCTGCCCTCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                  <  4:598184/272‑1 (MQ=255)
                                                                                                                                                                                     aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGc                                                                                                                                                                                                                                                                  <  3:5901/82‑1 (MQ=255)
                                                                                                                                                                                      gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                                                                              <  8:839344/149‑1 (MQ=255)
                                                                                                                                                                                        gCATTTTTATTGCAGGGAATACATTAATTGAATGTTAGGTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTAct                                                              <  1:22400/275‑1 (MQ=255)
                                                                                                                                                                                         cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                                                                           <  5:173899/149‑1 (MQ=255)
                                                                                                                                                                                          aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                          >  7:772629/1‑149 (MQ=255)
                                                                                                                                                                                          aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                                                                          >  7:377376/1‑149 (MQ=255)
                                                                                                                                                                                           tttttATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGCCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCAGTCCTGCCCGCCACCCTCCTGGCCGGCTCTTCCCAACGTTTAAATTCCCTCCCGTCCGTGTCTTCCTTCCCGGGGGAGTGTTCACCGACCAACAACCGCTAACACAAAAGCCCCCGTCTTCCGCTTTCACCTTCTGTTTTATTCTATTTCTGGCCGTTCCctcctcc                                                            >  4:300299/1‑273 (MQ=255)
                                                                                                                                                                                                    cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                                                                >  5:5755/1‑149 (MQ=255)
                                                                                                                                                                                                      ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgag                                                 >  2:432835/1‑274 (MQ=255)
                                                                                                                                                                                                       gTAATACATTAATTAAATGTTAGTGCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                <  1:165294/274‑1 (MQ=255)
                                                                                                                                                                                                           tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                                                                         >  8:272142/1‑149 (MQ=255)
                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  7:699050/1‑149 (MQ=255)
                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  6:627153/1‑149 (MQ=255)
                                                                                                                                                                                                              aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                                                                      >  8:662272/1‑149 (MQ=255)
                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                                                                     <  5:1102537/149‑1 (MQ=255)
                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                                                                                  >  2:288687/1‑184 (MQ=255)
                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                    >  2:467755/1‑262 (MQ=255)
                                                                                                                                                                                                               ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacacc                                       >  1:69112/1‑274 (MQ=255)
                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAATGCCATCCTTACGGATGGCCTTTCGCTTGATTTGATATCTGGCAGTTTATGGCGGACGGCCTTCCCGCCACCTTTCGGGCCGTTGCTTCTCAACGTTCAAATCCGCTCCCCCCCGCTTTGTcc                                                                                                                                                                 >  8:929686/1‑149 (MQ=18)
                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                                                 >  6:933647/1‑149 (MQ=255)
                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGgaga                                                <  4:467755/262‑1 (MQ=255)
                                                                                                                                                                                                                   ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGACAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                                                                 >  7:215526/1‑149 (MQ=255)
                                                                                                                                                                                                                          ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                                                                                          >  6:824060/1‑149 (MQ=255)
                                                                                                                                                                                                                            aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGga                                                                                                                                                        <  6:5755/149‑1 (MQ=255)
                                                                                                                                                                                                                               tagAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGagag                                                                                                                                                     <  6:651389/147‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaGCAAAAAGGCCATCCTGTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGATATTGtt                                                                                                                                                 <  8:105140/149‑1 (MQ=255)
                                                                                                                                                                                                                                    aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTc                                                                                                                                                >  7:832040/1‑149 (MQ=255)
                                                                                                                                                                                                                                    aaGCAAAAAGGCCATCCTTTCGGATGGCCCTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                  <  4:420096/275‑1 (MQ=255)
                                                                                                                                                                                                                                     aagcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                   <  6:488484/141‑1 (MQ=11)
                                                                                                                                                                                                                                     aagcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                                                                                   >  5:488484/6‑146 (MQ=11)
                                                                                                                                                                                                                                      gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg                                                                                                                     <  2:671337/174‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGa                                                                                                                                           >  7:992165/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaa                                                                                                                                        >  8:515901/1‑149 (MQ=255)
                                                                                                                                                                                                                                              gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                                                                                      <  8:992165/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                       <  8:543978/147‑1 (MQ=255)
                                                                                                                                                                                                                                               ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                                                                                       >  7:543978/1‑147 (MQ=255)
                                                                                                                                                                                                                                                cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa                                                                                                                                    >  7:179256/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAaacaac                                                                                                                                   >  6:443767/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTTGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCCCTACGGCGTTTACCTTCTGAGTTc     >  2:371001/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa                                                                                                                                <  6:285677/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tttGGGATGGCCTTTCGCTTGATTTGGGGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc  <  4:377037/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                 
AAAGAGGTGCTGGATACGATGATTGGGCTGGCGCAGTCGGGTATGACAATGTTGTGTGTAACACATGAGATGGGGTTTGCACGAACCGTCGCTGACCGGGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGC  >  NC_000913/3422985‑3423502

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: