Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,320,344 G→A A107V (GCC→GTC)  phnI ← carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,320,3440GA95.6% 183.7 / NA 70A107V (GCC→GTC) phnIcarbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (41/24);  minor base C (1/2);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCG  >  NC_000913/4320104‑4320594
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                              
cAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGgagctc                                                                                                                                                                                                                                                       >  1:333733/1‑246 (MQ=255)
    ggTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGGAAAATCGAGCAGGCGATGGGGGTCGTCGTAGGTTGGGCCAAGCAGCTGGGCGCCGGGGATGTCTTGATAAACGACGgagaactcgggg                                                                                                                                                                                                                                                 >  2:392656/1‑243 (MQ=255)
            ggcggcGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCCGGGGAAAATCGAGCAGGCGGTGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTTTAAACGACGGAGATACGGCGttggggacg                                                                                                                                                                                                                                         >  1:29904/1‑244 (MQ=255)
             gggggtggcggtggaggTCATCCGGCTGGGGGCCGTTATCCTCTTCAAACGTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGGTGGCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGTAATGGCTTTATAAACGACGGAGATACGGCGTTCGAGACGc                                                                                                                                                                                                                                       <  2:256585/235‑1 (MQ=255)
                 cgTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGGAAAATCGGGCAGGCGATGGGTGTAGTCGTAGGGTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATcc                                                                                                                                                                                                                                   >  1:205002/1‑251 (MQ=255)
                  gggctggggATGGCATCCGGCTGGGCGCTGCTACCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATccc                                                                                                                                                                                                                                  <  2:57881/243‑1 (MQ=255)
                   ggcgggttaggccATCCGGCTGTGCGCCGCTATCCTCGTCAATCTTCGCCAGCCCCGGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCg                                                                                                                                                                                                                                 <  2:142519/239‑1 (MQ=255)
                          gATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGGGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGgtgt                                                                                                                                                                                                                         >  1:323333/1‑251 (MQ=255)
                          gATGTCACCCGGCTGGGCGCCCCTATCATCTTAAACCTTCCCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGgtgt                                                                                                                                                                                                                         <  2:346676/251‑1 (MQ=255)
                                aTCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGGAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCg                                                                                                                                                                                                                   >  1:429627/1‑251 (MQ=255)
                                   cGGCTGTGGGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACTTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCt                                                                                                                                                                                                                <  1:157139/251‑1 (MQ=255)
                                    ggCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGGAAAATCGAGCAGGCGATGGGGGGAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGGAAGTCTTTATAAACGACGGAGATACGGCGTTCGAGGCGCATCCCGGTGGTGTcaggcggccc                                                                                                                                                                                                               >  1:18133/1‑242 (MQ=255)
                                       tgtgCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGCCGGCTGCTGTTCGCTGTCGGCGGGGGTCAGCGGCGGCGCTTCGCCGTTTGCCCGCCGGGGAAAATCGAGCAGGCGGTGGGGGTAGTCGGTGGTTGGGCCCAGCAGCTGGCCGCCGGGGATGTCTTTATAAACGACGGGGATACCGCGGTCGAGACGCATCCCGGTGGTGGCGGGCGGCCCGCt                                                                                                                                                                                                            >  2:325778/1‑251 (MQ=255)
                                                 tATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGGTGGGGGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGc                                                                                                                                                                                                  >  1:56985/1‑251 (MQ=255)
                                                    ccTCTTCAAACTCCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGCCGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTc                                                                                                                                                                                              <  2:346290/251‑1 (MQ=255)
                                                     ctctTCACACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCCTTTGCCAGCAAGGTAAAATCGAGCAGGCGATGGGTGTAGTCGGAGGGTGGGGCAAGCCGCTGGCCGCCGGGACTGTCTTTTTACACGACGGAGATACGGCGTTCGAGCCGAATCCCCGTGGTGGCGAGCGGccccccttcctccgcgttgt                                                                                                                                                                                              >  2:496911/1‑232 (MQ=255)
                                                       ctTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGGTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGGCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGcc                                                                                                                                                                                           >  1:211134/1‑251 (MQ=255)
                                                           aaaCTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACg                                                                                                                                                                                       >  1:313394/1‑251 (MQ=255)
                                                            aaCTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCAGCAGGGGAAAATCGAGCAGGCGGTGGGGGTAGTCGTCGGTTGGGCCAAGCCGCTTGACGCCGGGAAAGTCTTTATAAACGACGGAGATACGGCGGTCGAGACGCAGCCCGGGGTGGTCGAGGGGGTCGCTTGCCGCCACCGTCGcccagcgg                                                                                                                                                                                     >  2:413574/1‑245 (MQ=255)
                                                                 ccccAGCCCCTGACGCGCCAGCAGGCTGAAAGCGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGAGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACCACGGAGATACGGCTTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGc                                                                                                                                                                                 <  1:334559/249‑1 (MQ=255)
                                                                       ccccTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACTGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGTTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCTTGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAgg                                                                                                                                                                           <  1:104918/251‑1 (MQ=255)
                                                                                   cagcagGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAACAgg                                                                                                                                                                >  1:335974/1‑250 (MQ=255)
                                                                                       agGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCCGCAGGGTAAAATCAAGCAGGCGATGGGGGTAGTCGTAGGTTGGGGCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGGTACGGCGTTCGAGTCGCATCCCGGTGGTGTCGGGCGGCCCGCTTACCGCCCGCTTCGCCAACGTGATGCGGTAGGGGCGCAGCAGGaaaa                                                                                                                                                            >  1:254792/1‑250 (MQ=255)
                                                                                            gAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTGGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGc                                                                                                                                                        <  1:478602/249‑1 (MQ=255)
                                                                                                aCGTGCGGCGACGGCTGCGGTTCGGTGGCGGCGGGGGTCAGGGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAGCGAGCAGGCGATGGGTGTAGGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAAc                                                                                                                                                  <  2:250112/251‑1 (MQ=255)
                                                                                                   tGCGGCGCCGGCTGCGGTTCGCTGCCGGCGGGGGTCGGGGTCGGCGCTTCGCCGTTTGCCAGCAGGGAAAAATCGGGCAGGCGATGGGGGTAGTCGTAGGTGGGGCCAGGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGt                                                                                                                                                <  2:445778/250‑1 (MQ=255)
                                                                                                        cgACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGGCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGCGGTACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCCCGCTTCCCCCCACCCTCCCCCAACTAGTGCGGTAGGCGCCCAGCAGGAAAAACGCTTCACCGTTATccc                                                                                                                                          >  2:277779/1‑249 (MQ=255)
                                                                                                               tgctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGc                                                                                                                                   >  1:301844/1‑251 (MQ=255)
                                                                                                                 ctgTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCt                                                                                                                                 >  1:374277/1‑251 (MQ=255)
                                                                                                                   gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGGAAATCGCTTCAACGTTATCGCCGCTGGCCTGt                                                                                                                               >  1:483489/1‑251 (MQ=255)
                                                                                                                   gTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGt                                                                                                                               >  1:209385/1‑251 (MQ=255)
                                                                                                                      cGCTGTCGGGGGTGGGCGGCGGCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAAACGAGCGGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTAGAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTc                                                                                                                            <  2:320673/251‑1 (MQ=255)
                                                                                                                              gcggGGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGGGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcagcgccag                                                                                                                    <  2:446268/251‑1 (MQ=255)
                                                                                                                                    gtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCTCCGCTGGCCTGTTTCAACGCCAGcgccgc                                                                                                              >  2:137614/1‑251 (MQ=255)
                                                                                                                                    gtCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCCCCAACGTGATGCTGTAGGCGCCCAGCAGGAAAATCGCTTCAACAGTATCGCCGCTGGCCGGTTTTCGCGCCAGcgccgc                                                                                                              >  2:191077/1‑251 (MQ=255)
                                                                                                                                                   gcGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAACCGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCgg                                                                                                <  1:132534/249‑1 (MQ=255)
                                                                                                                                                      ccGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGCGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGGGATGCGGTAGGCGCGCCGCAGGACAATCGCTTCACAGTTATCACCGCTGGCCTGTTTCAGACCCCGCGCCGCCAGTTCGCGGTCGGCact                                                                                            >  2:132589/1‑249 (MQ=255)
                                                                                                                                                       cGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCCACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCCATg                                                                                           >  2:119279/1‑251 (MQ=255)
                                                                                                                                                        gtgtgccggcggggTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCAAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGAGACGGCGTTCGAGCCGCATCCCGGTGGTGCCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATgc                                                                                          <  2:313394/240‑1 (MQ=255)
                                                                                                                                                        gTTTGCCGGCAGGGTAAACTCGAGCAGGCGCTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATACACGACGGAGCTGCGGCGTTCGAGCCGCATCCCGGCGGCGTCGAACGGCTCGCTTACCGCACGCCTCGCCAACTCGATGCGGTAGGCGCGCGGCAAGACACTCGCTTCAACGTTCTCGGCGCTGGCCTGGTTCAGCGCCCGCGCCGCCAGTTCCCTGTAGGCAATgc                                                                                          >  2:45508/1‑251 (MQ=255)
                                                                                                                                                             ccAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCt                                                                                     >  1:464886/1‑251 (MQ=255)
                                                                                                                                                              cagcagGGTAAAACCGAGCAGGCGAGGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTCCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCtt                                                                                    <  1:170550/251‑1 (MQ=255)
                                                                                                                                                                gcagGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTCCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCg                                                                                  <  1:55902/251‑1 (MQ=255)
                                                                                                                                                                  agGGTAAAATCGAGAAGGAGGTGGGAGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGAATCCCGGTGGTGTAGAGAGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGTGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCCGAGCCAGCGCCGCCAGTTCGCGGTCGGGCATGCCGCCTTCggg                                                                                >  1:111540/1‑250 (MQ=255)
                                                                                                                                                                  agGGGAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGATGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCggg                                                                                >  2:452897/1‑250 (MQ=255)
                                                                                                                                                                       aaaaGCGAGCGGGGGAGGGGTGGAGTCGTGGGTGGGGCGAAGCGGCGGGCCGCCGGGACGGTCTTTATAAACTACGGAGATACGGGTTCCGGGACGCATCCCGGGGGGGTCGCGCGGCTCGCTGACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTCCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGtcatca                                                                           <  2:85012/251‑1 (MQ=255)
                                                                                                                                                                           ccggggaggggatggggtaagtggTAGGTGGGGCCAAGCAGCGGGCGGTCGGGAAGTTCTTTAAAAACGACGGAGATCCGGCGTCCGAGCCGCACCCGGGGGGTGTCGGGCGGCCCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCAcgcg                                                                       <  2:162606/228‑1 (MQ=255)
                                                                                                                                                                             gAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTCCGGTCATCACGCGAt                                                                     >  1:442742/1‑251 (MQ=255)
                                                                                                                                                                             gAGCAGGCGAGGGGTGTAGTCGTAGGTGGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGAt                                                                     <  1:285431/251‑1 (MQ=255)
                                                                                                                                                                              ggggggggatggggtgggTCGGAGGGTGGGCCAAGGGGCGGGCCGCCGGGAATGTCTTTATAACCGCGGGAGTTACGGGGCTCGAGACGCATCCCGGTGGTTTCGAGCGGTTCGCTTACCGCCAGCTTCGCCAACGTGAGGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTAGCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATc                                                                    <  2:254792/234‑1 (MQ=255)
                                                                                                                                                                                cAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCt                                                                   >  1:373128/1‑250 (MQ=255)
                                                                                                                                                                                 aGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAACTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAACTAGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCCGTTCGCGGTTGGCAAGTCCGCCCTCGGTCATCCCGCGATCTAc                                                                 >  2:363935/1‑251 (MQ=255)
                                                                                                                                                                                    cGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACcgc                                                              >  1:237470/1‑251 (MQ=255)
                                                                                                                                                                                     gATGGGTGTAGTCGTAGGTTGGGCCAAGCAGATGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCACCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCCACCTACcgcg                                                             >  2:498710/1‑251 (MQ=255)
                                                                                                                                                                                       tGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACcgcg                                                             >  1:357914/1‑249 (MQ=255)
                                                                                                                                                                                        gggTGTAGTCGTAGGTTGGGCCAAGCAGCGGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGACGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGTTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAg                                                           <  1:52396/250‑1 (MQ=255)
                                                                                                                                                                                                    aGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAActgctg                                               <  1:113479/250‑1 (MQ=255)
                                                                                                                                                                                                        tGGGCCAAGCATCGGGCCGCCGGGAAGTTATTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGGCGAGCGGCTCGCTTACCGCTAGCTTCTCCAACGTGATGCGGTAGGCGTGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCTTCACGCGCGCTACCGCGAGGTTAAACTTCTGTTCa                                           <  1:366598/250‑1 (MQ=255)
                                                                                                                                                                                                         gggCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAAt                                         >  1:35867/1‑251 (MQ=255)
                                                                                                                                                                                                              aaGCAGCTGGCCGCCGGGAATGTCTTTATAAACGCCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTggg                                    >  1:392160/1‑251 (MQ=255)
                                                                                                                                                                                                               agcagcTGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGGCATCTACCGCGGGGTTAAACAGCTGTTCAATCTgggg                                   >  2:386036/1‑250 (MQ=255)
                                                                                                                                                                                                                     tGGCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCt                             >  1:343039/1‑251 (MQ=255)
                                                                                                                                                                                                                      ggCCGCCGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTc                            >  1:302226/1‑251 (MQ=255)
                                                                                                                                                                                                                       gccgccGGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCCGGCCTGTTTCAGCGACAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCCAGGTTAATCTGCTGTTCAACCTGGGCGACGccc                            >  2:42310/1‑248 (MQ=255)
                                                                                                                                                                                                                         cgccgGGAATGTCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCCGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCCCGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGt                         >  2:194211/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccgGGAATGTCTTTATAAACGCGGTGGGTACGGCGTCCGTGCCGCATCCCGGTGGTTGCGAGCGGCTCGCTTCCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGAAGCAGGAAAATCGCTTCAACGTTATCGCCGCGGGCCTGTTTCAGCGCAAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTc                       <  2:501682/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggAATGCCTTTATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCACGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCggg                    <  1:194920/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tcttttATAACCGACGGACATACGGCTTCCGAGACGCATCCCGGTGGTGCCGAGCGGCGCGCCTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCaaa                <  2:50629/247‑1 (MQ=255)
                                                                                                                                                                                                                                      tttATAAACGACGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCaaa                >  1:50629/1‑247 (MQ=255)
                                                                                                                                                                                                                                       ttATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCGCGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGt           >  1:190061/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gcggGGATACGGCGTGCGGGACGCACCCCGGGGGTGTCGAGCGCCCCGCTACCCGCCAGCTTCGCCAACGTGATGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCGGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAACTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCg  <  2:429627/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                              
CAGGGGTAAACCGGCGGCGTGCGGGTGATGTCATCCGGCTGTGCGCCGCTATCCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCG  >  NC_000913/4320104‑4320594

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