Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,231 T→A 6.7% T84S (ACT→TCT)  dcyD ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,2310TA 6.7% 742.4 / 31.3 223T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (101/102);  new base A (8/7);  total (113/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

Only 100 of 227 total aligned reads displayed.

GCGCACTCTCCACATAACCTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGAT  >  NC_000913/1998961‑1999493
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
gcgcACTCTCCACATAACCTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAGGCCTTGTGCTTCGACTCGCGTCGCAAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGAAGTGCAGACCGAGTAGGGAAGCGACTGCGGGAGTCTg                                                                                                                                                                                                                                                                     >  2:138779/1‑274 (MQ=255)
    ccctccgcacggaCCTAGCGCGACCAGTGCATTACAACCGCCAACCGGAATGAAAGCCGGGCGAAAGCGTGGGTGTTGGAGTGGCGTCGCCGGGTCGGCGAGTTGGGCATTGGGATCGGTGAGTGCGTCGCACATTGCAAGTTGGGAAGTGAACGGAGCCAGCAACAAACGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGGCAATAGGATTTTCCAGCAGCGCCACGCAGAGCAGACCGAGTTTGGCGGCGACTGCGGCAGTCTGCCGca                                                                                                                                                                                                                                                                <  4:191895/263‑1 (MQ=255)
     ccctcCACATCAACTAGAGCGCCACAGGCATTCGACCCGCCAACCCGAATTACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                                                                                                                                                                                           <  1:234308/277‑1 (MQ=255)
                    aGCGCGCAAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGc                                                                                                                                                                                                                                                                    <  2:3554/255‑1 (MQ=255)
                      cgcgcCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGCAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCGATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAAGAGGATTTTCCAGCAGCCCCACGCCGTGCAGCCCGAGTTTGGCCGCCACTGCGGCAGTCTTCCGCACCATGTTAGACTggtacc                                                                                                                                                                                                                                              >  4:388194/1‑271 (MQ=255)
                                                               ggCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGca                                                                                                                                                                                                                                                                <  2:65241/216‑1 (MQ=255)
                                                                  gAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCCGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATACGCGTATCCGCCCCTTTACGCagggacacc                                                                                                                                                                                                  >  4:636952/1‑268 (MQ=255)
                                                                          tgtgCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCCCGCAGAGCATCTGCCGCgag                                                                                                                                                                                          >  2:674875/1‑275 (MQ=255)
                                                                                                           tGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTGGGTTAGATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGAGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGGCGGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTCACGGCCCCTTCACGCATAGCGGCGCACGCGAGGAATTCCCGCTTCCGTCATTTATTGACCCCCAt                                                                                                                                                         >  4:581244/1‑275 (MQ=255)
                                                                                                                          gTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGt                                                                                                                                                                                                                                                                        >  6:1452183/1‑149 (MQ=255)
                                                                                                                                   tCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcac                                                                                                                                                                                                                                                               >  7:1223829/1‑149 (MQ=255)
                                                                                                                                   tCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                                                                                                                                                                                           <  2:414699/169‑1 (MQ=255)
                                                                                                                                    cGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcaca                                                                                                                                                                                                                                                              <  8:2328/149‑1 (MQ=255)
                                                                                                                                    cGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTCCGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCg                                                                                                                                 >  2:329039/1‑274 (MQ=255)
                                                                                                                                       acaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGGCGCACATgg                                                                                                                                                                                                                                                           <  5:39475/149‑1 (MQ=255)
                                                                                                                                       acaTTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                                                                                                                                                                                           <  5:1682405/149‑1 (MQ=255)
                                                                                                                                           ttCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAg                                                                                                                                                                                                                                                       >  7:1109636/1‑149 (MQ=255)
                                                                                                                                            tCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGa                                                                                                                                                                                                                                                      <  7:1500250/149‑1 (MQ=255)
                                                                                                                                              aaTCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACt                                                                                                                                                                                                                                                    >  6:903984/1‑149 (MQ=255)
                                                                                                                                                tCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCTGCAGTCTGCCGCACATNGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCACAGCATCACCCGCGCGAATTTCCCNCTACAGTAATTTCTTACCGCCCATTCACACNNGGCTGACGTAATCCCGTTTCATCAACAttt                                                                                                                   >  3:110728/1‑276 (MQ=255)
                                                                                                                                                   gggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                                                                                                                                                                                               <  7:693890/149‑1 (MQ=255)
                                                                                                                                                   gggTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                                                                                                                                                                                               <  5:53131/149‑1 (MQ=255)
                                                                                                                                                     gTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGc                                                                                                                                                                                                                                             <  7:17480/149‑1 (MQ=255)
                                                                                                                                                        aagAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                                                                                                                                                                                          <  8:738628/147‑1 (MQ=255)
                                                                                                                                                           aaCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                                                                                                                                       >  6:1486668/1‑149 (MQ=255)
                                                                                                                                                           aaCAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                                                                                                                                                                                       >  5:600591/1‑149 (MQ=255)
                                                                                                                                                                aTCCAGCCACAAGCGATTGCCTTTGGTTAAAGAGTTTTCTGGGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGGGTTTGGCAGCGACTGAGGCAGTCTTCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAAt                                                                                                    <  4:348786/275‑1 (MQ=255)
                                                                                                                                                                aTCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTaa                                                                                                                                                                                                                                  >  5:340247/1‑149 (MQ=255)
                                                                                                                                                                aTCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTaa                                                                                                                                                                                                                                  <  5:734059/149‑1 (MQ=255)
                                                                                                                                                                    aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                                                                                                                                              >  7:338198/1‑149 (MQ=255)
                                                                                                                                                                    aGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAg                                                                                                                                                                                                                              >  6:1672957/1‑149 (MQ=255)
                                                                                                                                                                     gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                                                                                                                                                                                           >  1:216338/1‑135 (MQ=255)
                                                                                                                                                                     gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                                                                                                                                                                                           <  3:216338/135‑1 (MQ=255)
                                                                                                                                                                     gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                                                                                                                                                                                             <  6:574128/149‑1 (MQ=255)
                                                                                                                                                                     gCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt                                                                                                                                                                               >  2:32150/1‑195 (MQ=255)
                                                                                                                                                                          aagagtgtGCGGTTGGTTAAATGGTTTTCTGGGGGTGTGCCGAAAGGATTTTCCAGCAGGGACACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGAACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAgcg                                                                                          <  4:502863/268‑1 (MQ=255)
                                                                                                                                                                          aaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                                                                                                                                                                                        >  8:1680548/1‑149 (MQ=255)
                                                                                                                                                                          aaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                                                                                                                                                                                        <  8:918752/149‑1 (MQ=255)
                                                                                                                                                                          aaGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                                                                                                                                                                                        >  5:1238992/1‑149 (MQ=255)
                                                                                                                                                                                ttGTCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                                                                                                                                                                                                  >  6:211482/1‑149 (MQ=255)
                                                                                                                                                                                     gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTgg                                                                                                                                                                                                                                                  >  8:1113451/1‑112 (MQ=255)
                                                                                                                                                                                        ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCa                                                                                                                                                                                                          >  7:572593/1‑149 (MQ=255)
                                                                                                                                                                                          ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGTCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                                                                                                                                        >  8:1057549/1‑149 (MQ=255)
                                                                                                                                                                                          ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                                                                                                                                        >  8:635383/1‑149 (MQ=255)
                                                                                                                                                                                          ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                                                                                                                                                                                        >  8:487248/1‑149 (MQ=255)
                                                                                                                                                                                            aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGc                                                                                                                                                                                                      <  7:818159/149‑1 (MQ=255)
                                                                                                                                                                                               tAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                                                                                                                                                                                                   >  5:1014787/1‑149 (MQ=255)
                                                                                                                                                                                                 gTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                                                                                                                                                                                                <  5:1507389/118‑1 (MQ=255)
                                                                                                                                                                                                   tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                                                                                                                                                                                             <  6:950592/135‑1 (MQ=255)
                                                                                                                                                                                                    ttCTGCGGTTGTTCCAATAGGATTTTCCAGCAGGGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGgcg                                                                <  2:265735/275‑1 (MQ=255)
                                                                                                                                                                                                       tGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAAAGCATCTGCCGCga                                                                                                                                                                                           <  7:128605/149‑1 (MQ=255)
                                                                                                                                                                                                           gTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaaa                                                                                                                                                                                       >  7:1600153/1‑149 (MQ=255)
                                                                                                                                                                                                               tgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc                                                                                                                                                                                   >  7:620335/1‑149 (MQ=255)
                                                                                                                                                                                                                  cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCTCGCAGAGCATCTGCCGCGAGAAATTCCAGc                                                                                                                                                                                >  6:1038713/1‑149 (MQ=255)
                                                                                                                                                                                                                   aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTa                                                                                                                                                                                                                                                        <  5:908524/76‑1 (MQ=255)
                                                                                                                                                                                                                   aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                                                                                                                                                                                               >  8:1312432/1‑85 (MQ=255)
                                                                                                                                                                                                                     tAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTa                                                                                                                                                                             >  8:1499047/1‑149 (MQ=255)
                                                                                                                                                                                                                      aGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGACATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGCGCGGCGTTGGCGCGCCGATAAACTCCAGCAGTgg                                          >  1:153510/1‑279 (MQ=255)
                                                                                                                                                                                                                        gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                                                                                                                                  >  7:1685027/1‑141 (MQ=255)
                                                                                                                                                                                                                        gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                                                                                                                                  <  8:1685027/141‑1 (MQ=255)
                                                                                                                                                                                                                               cagcagCGCCACGCAATGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAt                                                                                                                                                                   >  8:1616344/1‑149 (MQ=255)
                                                                                                                                                                                                                                agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                                                                                                                                  >  5:1122271/1‑149 (MQ=255)
                                                                                                                                                                                                                                agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                                                                                                                                  >  7:209620/1‑149 (MQ=255)
                                                                                                                                                                                                                                agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTAtt                                                                                                                                                                  >  6:1084259/1‑149 (MQ=255)
                                                                                                                                                                                                                                 gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                                                                                                                                                                                             >  5:1529824/1‑57 (MQ=255)
                                                                                                                                                                                                                                 gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                                                                                                                                                                                             <  6:1529824/57‑1 (MQ=255)
                                                                                                                                                                                                                                  ccgcgcCACGCAGGGCAAACCGAGTTTGGCATCGACATCGGCATGAGTACGAAAATGGTTGAACTGGGTAGCCCCGGCAGTAATAAGCGTATCTGCACCTTAACGCAGAGCATCTGCCGCGAGAAATGCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACggg                                  <  4:329035/273‑1 (MQ=255)
                                                                                                                                                                                                                                  cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgc                                                                                                                                                                <  5:279069/149‑1 (MQ=255)
                                                                                                                                                                                                                                   tgagtccgggcgggCAGATCGGGTTTGGGGTCGCGTGGGGCAGTCTGCTGCCAATGGTTAGACTGGATCGGCCCGTAAGGAATAACCGGATATGCACATTCACGCAGGGGATCTGCCGCGAGAAATTCCAGCGTACGTAATTTATTGCAGCCCATGGCCATGGGGGGGACGTAATCCCTTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACggg                                  <  4:422200/261‑1 (MQ=255)
                                                                                                                                                                                                                                   agcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgcc                                                                                                                                                               <  5:1084259/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGTCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTgccgcc                                                                                                                                                            <  8:1143509/149‑1 (MQ=255)
                                                                                                                                                                                                                                       ccACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGccc                                                                                                                                                           <  5:1224991/149‑1 (MQ=255)
                                                                                                                                                                                                                                          cGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCAtt                                                                                                                                                        <  5:115469/149‑1 (MQ=255)
                                                                                                                                                                                                                                           acAGTGCAGACCGAGTTTGGGAGCGACTGCGGCAGTGTGCCGAAAATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTTCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTg                                                                                                                                                       >  8:108092/2‑149 (MQ=255)
                                                                                                                                                                                                                                               tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                                                                                                                                   <  8:572593/149‑1 (MQ=255)
                                                                                                                                                                                                                                               tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                                                                                                                                   <  6:639459/149‑1 (MQ=255)
                                                                                                                                                                                                                                               tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCAt                                                                                                                                                   >  6:479842/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 cAGACCGNGTTTGGCGNCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCa                                                                                                                                                    <  3:54991/146‑1 (MQ=255)
                                                                                                                                                                                                                                                   gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATgggg                                                                                                                                               <  5:685330/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    aCCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGATTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGTTAATCAGAAAAGCGCGGCAGATCTTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTCTGCAGTGGCa                >  2:455962/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     ccGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGt                                                                                                                                             >  5:117449/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        aGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGCGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGAc                                                                                                                                          <  6:365943/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         gtctGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGCTAGTCCGAAAAGCGCGGCAGATATTCCCGCGGAGTTGGCCCGCCGATCAACTCCAGCCGTAGAAATCGGGTTAAATTATGCTATGGCaaaaac           >  4:80861/4‑273 (MQ=255)
                                                                                                                                                                                                                                                         gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt                                                                                                                                                                               >  2:510270/1‑111 (MQ=255)
                                                                                                                                                                                                                                                            tctaaccatgtGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTCAATTATGCAGTGGCATAACAGCCTCCg    >  1:293187/11‑279 (MQ=255)
                                                                                                                                                                                                                                                             ggCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCAt                                                                                                                                     <  5:843059/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                                                                                                                                  <  7:1417709/103‑1 (MQ=255)
                                                                                                                                                                                                                                                              gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                                                                                                                                  >  8:1417709/1‑103 (MQ=255)
                                                                                                                                                                                                                                                               cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATcc                                                                                                                                   <  8:811675/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATccc                                                                                                                                  >  6:1598656/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aGCGACGGCGGCAGTCTGCCGAACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTcc     <  2:636957/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCg                                                                                                                                 >  7:928176/6‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt                                                                                                                               <  6:1154692/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt                                                                                                                               <  8:445888/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGtt                                                                                                                               <  5:1452183/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg                                                                                                                             >  7:50453/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa                                                                                                                                                                                  >  6:314283/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                      tGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtga                                                                                                                            <  7:1049919/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           cAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGaaa                                                                                                                       >  6:1317201/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              tCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAtt                                                                                                                    >  5:874145/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                      
GCGCACTCTCCACATAACCTAGCGCGCCAAGAGCATTAGAACCGCCAACCGGAATGACATACGGGCGAAAGCCTTGTGCTTCGACTCGCGTCGCCAGCTCTTCCAGTTGGGCATTGGGATCGGTCAGTGCGTCGCACATTTCAATCTGGGTATTGAACAGATCCAGCAACAAGCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGAT  >  NC_000913/1998961‑1999493

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