Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,256,157 C→T Q165Q (CAG→CAA cyuA ← putative L‑cysteine desulfidase CyuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,256,1570CT98.7% 246.4 / ‑3.9 76Q165Q (CAG→CAAcyuAputative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/35);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCT  >  NC_000913/3255911‑3256397
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
ggATGAAGAGAGAATTGTCGCCAGCAAACCGCGTGCGCACTTTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCTCTTACCGCTCAAACCTACCTGCGATAACGCACAATTGAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCGCAGCCAGCGTCGTTCTGGAAAGCGCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACGCACGCCTGGTggg                                                                                                                                                                                                                                              <  2:502803/251‑1 (MQ=255)
  aTGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGt                                                                                                                                                                                                                                             >  1:101/1‑250 (MQ=255)
  aTGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGt                                                                                                                                                                                                                                             >  2:15635/1‑250 (MQ=255)
     aaGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAac                                                                                                                                                                                                                                         <  2:311211/251‑1 (MQ=255)
             tCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCTCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCCCAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCCTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGACAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAacaccacacc                                                                                                                                                                                                                                 >  2:361299/1‑251 (MQ=255)
                   gCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTg                                                                                                                                                                                                                           >  2:399402/1‑251 (MQ=255)
                   gCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTg                                                                                                                                                                                                                           <  2:471793/251‑1 (MQ=255)
                     cAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGGAAACACCACACCATCGTg                                                                                                                                                                                                                           >  1:289578/1‑249 (MQ=255)
                            ccGCGCTCGCCCTGATTTTACAGCGTCGCGCCACTATGCAGCCCACACTTACCGCTCAAACATTCCTGCGATATAGCACAATTTAGCTTCGCGGCCTCGAGCATAAAGCGGCTCGCCGCAAACGGGACTTCATTGACGAACTTCAGGCTCTCAGCCAGCGTCGTTCTGGAAGGCCCCGTCAGCGGAGACTCCTGCTCGCCCTCAGCCACACACGCCCGTTGGATAAACACCACACCATCGTGCgacacga                                                                                                                                                                                                                   >  2:532421/1‑244 (MQ=255)
                              gcgcTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACTGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGatat                                                                                                                                                                                                                <  2:267979/251‑1 (MQ=255)
                               cgcTCGCACTGTTTTTCCAGCGTCGCGCCATTATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGAAAACTCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCATCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCGCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATg                                                                                                                                                                                                               <  2:450143/251‑1 (MQ=255)
                                  ccGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCGTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACGTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACTCCATCGTGCGTCTCGATATGcac                                                                                                                                                                                                            <  2:69262/250‑1 (MQ=255)
                                    gCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATCTGcaca                                                                                                                                                                                                           >  2:443768/1‑250 (MQ=255)
                                         gTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGt                                                                                                                                                                                                       <  1:472041/249‑1 (MQ=255)
                                            gttcccgagttgcgcCAATATGCAGCCCCCACTTACCGCGCAAACCTTCCTGCGATAACGCACAAGTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTAATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCGGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGt                                                                                                                                                                                                   <  2:478897/239‑1 (MQ=255)
                                              tCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGCTAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGGTCTCAGCCAGCGTCGTTCTGGCTAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAGACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATg                                                                                                                                                                                                >  2:206348/1‑251 (MQ=255)
                                                   cGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAAGTTAGCGTCGCGGAATCGAGAATAAAGCGGATCGCCGCAATCGGGACTTCAGTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTGGGTATGCccgc                                                                                                                                                                                           <  2:101/251‑1 (MQ=255)
                                                        cgcCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgacg                                                                                                                                                                                      <  1:328678/251‑1 (MQ=255)
                                                                      ccccACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGc                                                                                                                                                                        >  2:557238/1‑251 (MQ=255)
                                                                           cTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACt                                                                                                                                                                   >  2:277508/1‑251 (MQ=255)
                                                                             tACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTtc                                                                                                                                                                 >  1:224436/1‑251 (MQ=255)
                                                                             tACCGCTCAAACCTTCATGAGAAAACGCACACTTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTATGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTtc                                                                                                                                                                 >  2:511006/1‑251 (MQ=255)
                                                                               cctctccaaCCTTCCTGCGATAACGAACCATTTAGCTTCGCGGAATCGAGAATAAAGCGCATCGCCGCAAACGGGACTTCATTGACGAACGTCAGGACCTCAGCCAGCGTCGTTCTGGAAAGCTCCGTCAGCGGAGACTCTTGCTCGCCCCCTGCCACTCACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTtctc                                                                                                                                                               <  2:12999/244‑1 (MQ=255)
                                                                                   tCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCCCTTCTCACCg                                                                                                                                                           >  1:309278/1‑251 (MQ=255)
                                                                                     aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGACTAAAGCGGATCGCCGCAAACGGGACTACATTGACGAACTGCAGGATCTCAGCCAGAGGCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGACCTCTGCCACACACTCCTGTTGGGTAAACACCACACAATCGTGCGTCTCGATATGCACAATGTTGGAAAGCCCGCCGACGATGGTGACACACGCCCACTactccacgtt                                                                                                                                                         >  1:494813/1‑245 (MQ=255)
                                                                                     aaaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGtt                                                                                                                                                         >  1:76059/1‑251 (MQ=255)
                                                                                      aaCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTc                                                                                                                                                        >  1:385455/1‑251 (MQ=255)
                                                                                       aCCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGAGCTCAGCCAGCGTCGTTCTGGAAAGCACCTTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTcc                                                                                                                                                       <  1:411819/251‑1 (MQ=255)
                                                                                        ccTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCa                                                                                                                                                      <  1:216454/251‑1 (MQ=255)
                                                                                           tccggCGATAACGCACAATTGATCCTCGCGGAATCGGGAATAAAGCGTATCTCCGCAAACGGGACTTCATTAACGAACTTCAGTATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTTATCGGATGCTCGTGGTCGCCCTCCGCCACACACTCGTGTGGGGTAAACACCACTCCATCGTGCGTCGCGATATGCACAATGGGGGTATGTCCGCCGACGATGGTGACACACGCCCGCTtctcccccctccag                                                                                                                                                     <  2:554831/245‑11 (MQ=255)
                                                                                           tCCTGCGATAACGCACAATGTAGCTTCGCGGAATCGAGAATAACGCGGACCGCCCAAAACGGGTCTTCATGGACGAACGTCAGGAGCTCAGCCGGCGTCGTTCTGGAAAGCTCCGTCAGCGGAGCTTCTTGCTCGCCCTCTGCCTCACACGCCGTTTGGGTAATCACCATCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCGCACCGTTCCAGa                                                                                                                                                    <  2:324240/250‑1 (MQ=255)
                                                                                               gCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTg                                                                                                                                               <  1:34589/251‑1 (MQ=255)
                                                                                                 gATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGCCTCTTGCTCGCCCCCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGgc                                                                                                                                             <  2:76059/251‑1 (MQ=255)
                                                                                                   tAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGCGACTCTTGCTCGCCCTCTGCCACACACGCCTATTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAAAGTTGGTATGCCAGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGgcgc                                                                                                                                           >  1:562904/1‑251 (MQ=255)
                                                                                                          caATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTgaga                                                                                                                                     >  1:108589/1‑250 (MQ=255)
                                                                                                              ttAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCTCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAagag                                                                                                                                 <  2:543224/250‑1 (MQ=255)
                                                                                                              ttAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCAGTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAagag                                                                                                                                 <  1:81559/250‑1 (MQ=255)
                                                                                                                           aTCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCa                                                                                                                     >  1:531244/1‑249 (MQ=255)
                                                                                                                            tCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGCCGATGGTGCCACACGCCCACTTCCCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTCATCGCAAgg                                                                                                                  >  2:142287/1‑251 (MQ=255)
                                                                                                                              gagaATAAAGGGGATCGCCGCAAACGGGACTGCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGtt                                                                                                                <  2:496817/251‑1 (MQ=255)
                                                                                                                              gagaATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATAACAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATACACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGCTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGtt                                                                                                                >  1:288558/1‑251 (MQ=255)
                                                                                                                                  aTAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTCATCGCGAGGTTCCTg                                                                                                            >  2:196447/1‑251 (MQ=255)
                                                                                                                                         ggATCGCCGAAATCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGAGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCCTCGTGCGTCTCGATATGCACACTGTTGGTATGCCCGCCGAAGATGGTGACACACGCCCACTTCTCACCGTACCATCCCTTGGCGCGTGAGAAGAGGATTTCCTCGCAAGGTTCCTGGATCtt                                                                                                      >  2:300773/1‑250 (MQ=255)
                                                                                                                                          gATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACTATGGTTACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTaa                                                                                                    <  1:360630/251‑1 (MQ=255)
                                                                                                                                             cgccgcAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACgg                                                                                                 >  1:205685/1‑251 (MQ=255)
                                                                                                                                                    aaCGGGACTTCCTTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                            >  2:220803/1‑217 (MQ=255)
                                                                                                                                                    aaCGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACTCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGAttt                                                                                                                            <  1:220803/217‑1 (MQ=255)
                                                                                                                                                      cGGGACTTCATTGAAGAACTTCAGGATCTCAGCCAGCGTCGTTATGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTACACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGATCTTAACGGAGACTTTcc                                                                                        >  2:246288/1‑251 (MQ=255)
                                                                                                                                                       gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGCAAGCACCGTCAGCGGAGACTCTTGCTCGCCATCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCCCAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTccc                                                                                       >  1:261503/1‑251 (MQ=255)
                                                                                                                                                       gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCACCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACCTCTCACCGTTCCAGACTTTGTCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTTGATCTTAACGGAGACTTTccc                                                                                       <  2:466421/251‑1 (MQ=255)
                                                                                                                                                       gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTccc                                                                                       <  1:191402/251‑1 (MQ=255)
                                                                                                                                                       gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGAACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTccc                                                                                       >  1:252490/1‑251 (MQ=255)
                                                                                                                                                          aCTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgc                                                                                    >  1:315716/1‑251 (MQ=255)
                                                                                                                                                            ttCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCccgccg                                                                                  >  1:110662/1‑251 (MQ=255)
                                                                                                                                                                 tGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcagc                                                                             >  2:429212/1‑251 (MQ=255)
                                                                                                                                                                   aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCcagc                                                                             >  2:471067/1‑249 (MQ=255)
                                                                                                                                                                    tgAACTTCAGGATGTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGAATCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCCTCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGCAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCACCACcagccgg                                                                          >  2:26403/2‑248 (MQ=255)
                                                                                                                                                                          ccAGGACCTCAGCCAGCGTCGTTCGGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACCCCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCAGTGCttt                                                                     <  2:108589/249‑1 (MQ=255)
                                                                                                                                                                           cAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTgg                                                                   >  1:325427/1‑251 (MQ=255)
                                                                                                                                                                               atatcacccagcGTCGTTCTCGCAAGGGCCGTCAGCGGAGACCATTGATCGCCCTATGCCACACACACCTGTTGGGTAAACACCACACCATCGTGCCTCTCGTTATGCCCAATGTTGGAATGCCCGCCGACGATGGAGACACACGCCCCCTTCTCACCGTTCCACACTTTGGCGCGTGAGCAGAGGATTTCATCGCAAGGTTCCTGCAGCTTCACGGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCagc                                                               >  2:63472/8‑249 (MQ=255)
                                                                                                                                                                                  tcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCCGCAGTGCTTTGGCATCGGc                                                            >  2:212508/1‑251 (MQ=255)
                                                                                                                                                                                   cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATc                                                               <  1:451076/247‑1 (MQ=255)
                                                                                                                                                                                   cagcaagcGTCGTTCTGGAAAGCTCCGTCAGCGGAGCCTCTGGCTCGCCCTTGGCCCCACACGCCTGTTGGGTAAACACCGCCCCATCGTGCGTCGCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCACCGTTCCAGTCTTGGGCGCGTGAGAATAGGATTTCATCGCATGGTGCCTGGATCTGAACGGAGACGTGCCCCGCCGCCAGCAGTGCTTTTGCATCGGCa                                                           <  2:197786/251‑1 (MQ=255)
                                                                                                                                                                                        accgtcgtTCGGGAAAGCCCCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACTCGCCTGTTGGGTAAACACCACCCCATCGTGCGTCGCGATATGCACAATGTTGTTATGCCCGCCGACGATGGTGACACACGCCCTCTTCCCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGTAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGc                                                      <  1:511006/249‑1 (MQ=255)
                                                                                                                                                                                          cgtcgtTCTGGAAAGCCCCGTCAGCGGAGCCTCTTGCTCGCCCTCTGCCACACACCCCTGTTGGGTAAACCCCACCCCATCGTGCGTCCCGATATGCACAATGTTGGTATGCCCGCCGACGAGGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCt                                                    <  2:315716/251‑1 (MQ=255)
                                                                                                                                                                                              gtTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCcgac                                                                                                                                                                                       >  2:402382/1‑116 (MQ=255)
                                                                                                                                                                                                tCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACTACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCt                                               <  1:440240/250‑1 (MQ=255)
                                                                                                                                                                                                   ggAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCGGTTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTGCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGt                                             <  1:272324/249‑1 (MQ=255)
                                                                                                                                                                                                     aaaGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCg                                         >  1:245142/1‑251 (MQ=255)
                                                                                                                                                                                                        gCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGATGGGTAAACACCACACCATCGTGCGTCTCGATATGCAAACTGTTGGAATGCCCGCCGACGATGGTGCCACACGCCCACTTATCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCCGTTGCGTCt                                      >  2:563153/1‑251 (MQ=255)
                                                                                                                                                                                                             gTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCCCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGCTTTCATCGCAAGGTTCCTGGATCTTAACGGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTccg                                 >  2:412546/1‑249 (MQ=255)
                                                                                                                                                                                                                             gCTCGCCCTCTGCCACACCCGCCTGTTGGGTAAACCCCCCCCCATCGGGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg                 <  1:557238/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cccctGCCACACCCGCCGGTTGGGTAAACCCCACACCAGCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTGCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCa           <  1:206348/248‑1 (MQ=255)
                                                                                                                                                                                                                                        ccaacacacGCCTGTTGGTTAAACACCCCCCCATCGTGCTACCCGATACTCAAAATGTGGTTATGCCCGCCGTCGAGGGTGACACACGCCCACTTCTCACAGTACCAGACTTTGGCGCGTGAGAAGAGGATTCCACCGCAAGGTTCCGGGATCTTACCGGAGCCTTTCCCCGCCGCCAGCAGTGCTCTGGCATCGGCAATTGCTTGCGCTGTTGCGTCTTTCAGCACTTCCAGCACGGCGTCGGCAttt        <  2:392127/246‑1 (MQ=255)
                                                                                                                                                                                                                                           ccccacGCCTGTTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCAcc   <  1:399402/248‑1 (MQ=255)
                                                                                                                                                                                                                                            cacacGCCTGTTGGGTAAACACCACCCCACCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGAGGGTGAAACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCt  <  1:121227/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GGATGAAGAGAGATCTTTCGCCAGCAAACCGCGCTCGCACTGTTTTTCCAGCGTCGCGCCAATATGCAGCCCCCACTTACCGCTCAAACCTTCCTGCGATAACGCACAATTTAGCTTCGCGGAATCGAGAATAAAGCGGATCGCCGCAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCT  >  NC_000913/3255911‑3256397

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