Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,684,019 C→T 5.5% R380W (CGG→TGG)  rstB → sensory histidine kinase of RstAB two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,684,0190CT 5.5% 480.2 / 11.9 146R380W (CGG→TGG) rstBsensory histidine kinase of RstAB two‑component system
Reads supporting (aligned to +/- strand):  ref base C (68/70);  new base T (4/4);  total (72/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 146 total aligned reads displayed.

CATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTATGGCGCGTTACGATCTCGTAGAC  >  NC_000913/1683746‑1684282
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                       
cATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCAAGAAAACCGCGAACATATCTTTGAACCTTTTGTTACCCCCGATCCCACCCg                                                                                                                                                                                                                                                                        >  4:438768/1‑275 (MQ=255)
                            ccGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                                                                                                                                                            <  2:190555/275‑1 (MQ=255)
                                        tACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCACAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCTATCCCAACCGGGACCACTCAACCGGAGGCTGCGGTCTGGGGCTTGCCAtt                                                                                                                                                                                                                                >  4:131190/1‑275 (MQ=255)
                                           ggATAAAAACGATCGGGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTAATGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                                                                                                                                                                                                             <  2:457011/275‑1 (MQ=255)
                                           ggAAAAAAACGTACGTGCAAGGCCTTTATGCGGCGTTGGATATGCGCTTAATGGAGGCCGTGCTGGATAATTTGCTCAATACCGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTc                                                                                                                                                                                                                             <  2:429286/275‑1 (MQ=255)
                                              gagaaacgtcggTGCAAGGCAAGTAGGCGGTGGTGGAGTTGAGGTTAATGGGGCGCGTGCTGGATAATTTGCTCAATAACGGCCTGCGGTACGGCCAGTAAACAGTTGAACCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCGAGGGAGTGCCCCAGAAAACCGCGAACATATCGTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCa                                                                                                                                                                                                                           <  4:464006/263‑1 (MQ=255)
                                                       tCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGTTGCGGGCGGGGGCTGGCAATTGTCCACTCTATAGcc                                                                                                                                                                                                                 >  1:133773/1‑274 (MQ=255)
                                                       tCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                                                                                                                                                                                                 <  2:361483/275‑1 (MQ=255)
                                                         gTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGCCAATTGTCCACTCTATAGCAc                                                                                                                                                                                                                >  2:400910/1‑274 (MQ=255)
                                                                                      gcgcTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCCAAACCGGCGGCTGCGGGCTGGGGCTGGCCATTGTCCACCCTTTAGCACTGGCAATGGGCGGTACGGTTAATTTTGaac                                                                                                                                                                                  >  4:113804/1‑273 (MQ=255)
                                                                                                 aGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGTATTGCCCCAGAAAACCGCGTACATATCTTTGAACCTTTTGTTCGCCTCGCTCCCAGCCGGGCTCGCTCAACCGGCGGCTGCGGTCTGGGGCTGTCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTACTTGTGCCGCCAGCGaatt                                                                                                                                                                        >  1:36342/1‑272 (MQ=255)
                                                                                                   cgcgTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGCTTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCCATGGGCGGTACGGTTAATTGTGACACCAGCGAACTgg                                                                                                                                                                      >  1:73263/1‑274 (MQ=255)
                                                                                                                 tttGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCGCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGACATTGTCCCCTCTATAGCCCTGGCAAGGGGCGGTACGGTTAATTGTTACACCAGAGAACTGGGTGGTGccgccgtgt                                                                                                                                                       >  4:523812/1‑268 (MQ=255)
                                                                                                                      tCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATggg                                                                                                                                                                                                      >  2:88452/1‑223 (MQ=255)
                                                                                                                      tCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTTTGACACCAGGGAACTGGGGGGTTCCCGCTTCCCGttt                                                                                                                                                  >  4:443370/1‑275 (MQ=255)
                                                                                                                          tAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGACCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGt                                                                                                                                                                                                  <  4:88451/223‑1 (MQ=255)
                                                                                                                                        aCTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTc                                                                                                                                                                                                                                                               <  6:906774/148‑1 (MQ=255)
                                                                                                                                             cATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                                                                                                                                                                                                  >  5:599545/1‑140 (MQ=255)
                                                                                                                                             cATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                                                                                                                                                                                                  <  6:599545/140‑1 (MQ=255)
                                                                                                                                               ttCACCGGTTGAATCCAGCCTGCTACTGTCGGGGAATAAAGCGACATTAATTGTTGAGGATGATGGCCAAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATcc                                                                                                                         <  2:284050/275‑1 (MQ=255)
                                                                                                                                                   aCGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGc                                                                                                                                                                                                                                                                 <  7:915394/135‑1 (MQ=255)
                                                                                                                                                        tGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCggg                                                                                                                                                                                                                                               >  8:1133942/1‑148 (MQ=255)
                                                                                                                                                           aaCCAGCCTGCTACTGTCGGGGAATAGAGCGACATGAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACct                                                                                                             <  2:523813/275‑1 (MQ=255)
                                                                                                                                                            aCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGa                                                                                                                                                                                                                                                                     >  6:465881/1‑122 (MQ=255)
                                                                                                                                                            aCCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGa                                                                                                                                                                                                                                                                     <  5:465881/122‑1 (MQ=255)
                                                                                                                                                                 ccTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGcc                                                                                                                                          >  2:29480/1‑240 (MQ=255)
                                                                                                                                                                  cTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCaa                                                                                                                                                                                                                                                             <  8:675783/124‑1 (MQ=255)
                                                                                                                                                                  cTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCaa                                                                                                                                                                                                                                                             >  7:675783/1‑124 (MQ=255)
                                                                                                                                                                    gCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCggg                                                                                                                                                                                                                                                                      <  5:1096413/113‑1 (MQ=255)
                                                                                                                                                                                  aTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggctg                                                                                                                                                                                                                                            >  6:477547/1‑125 (MQ=255)
                                                                                                                                                                                  aTAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtat                                                                                                                                                                                                                     <  7:1133942/148‑1 (MQ=255)
                                                                                                                                                                                    agagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCg                                                                                                                                                                                                                                                                  >  5:554479/1‑101 (MQ=255)
                                                                                                                                                                                     gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTg                                                                                                                                                                                                                               >  4:225307/1‑135 (MQ=255)
                                                                                                                                                                                     gagCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCAAGATgg                                                                                    >  1:94911/1‑274 (MQ=255)
                                                                                                                                                                                         gACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCa                                                                                                                                                                                                                           <  2:225308/135‑1 (MQ=255)
                                                                                                                                                                                           cATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGAGACCAGCGAACTGGGTGTTGCCCGCTTCTCGTTTCGCTGGCCGTTATGGCATAACAGCCCGCAATTTACCTCTGACTGACACTACGCGCACGATGGTCAAGTc                                                                             >  2:239260/1‑275 (MQ=255)
                                                                                                                                                                                                    ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                                                                                                                                                                                                 >  8:900488/1‑134 (MQ=255)
                                                                                                                                                                                                    ttgAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCa                                                                                                                                                                                                                 <  7:900488/134‑1 (MQ=255)
                                                                                                                                                                                                      gAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGAACCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGcc                                                                                                                                          >  2:166031/1‑203 (MQ=255)
                                                                                                                                                                                                        ggATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCgggct                                                                                                                                                                                                                                             >  7:524342/1‑102 (MQ=255)
                                                                                                                                                                                                          atgatgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGAACCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTa                                                                                                                                      <  4:166030/203‑1 (MQ=255)
                                                                                                                                                                                                              tgGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAAc                                                                                                                             >  2:314445/1‑208 (MQ=255)
                                                                                                                                                                                                                 cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTTGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCCTCTCGTTTAGCTGGCCGTCATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACtgt                                                                  >  2:297311/1‑264 (MQ=255)
                                                                                                                                                                                                                 cccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACCCTACGCGCACGATGGTCAAGTCACCACGACTTTGCTATAAAATa                                                       >  2:175241/1‑275 (MQ=255)
                                                                                                                                                                                                                  ccAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATcc                                                                                                                         <  4:314444/208‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                                                                                                                                                                                                    <  7:730601/116‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                                                                                                                                                                                                    >  8:730601/1‑116 (MQ=255)
                                                                                                                                                                                                                    aGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac                                                                                                                                                                                   >  7:676093/1‑148 (MQ=255)
                                                                                                                                                                                                                     gggATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                                                                                                                                          <  7:295811/124‑1 (MQ=255)
                                                                                                                                                                                                                     gggATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCaa                                                                                                                                                                                                          >  8:295811/1‑124 (MQ=255)
                                                                                                                                                                                                                     gggATTGCCCAAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGt                                                    <  2:131191/274‑1 (MQ=255)
                                                                                                                                                                                                                        aTTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAct                                                                                                                                                                                                                         <  5:776530/106‑1 (MQ=255)
                                                                                                                                                                                                                        aTTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCAct                                                                                                                                                                                                                         >  6:776530/1‑106 (MQ=255)
                                                                                                                                                                                                                             cccAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                                                                                                                                                                                                 >  5:705660/1‑93 (MQ=255)
                                                                                                                                                                                                                                aGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                                                                                                                                                                                                    >  8:1120953/1‑103 (MQ=255)
                                                                                                                                                                                                                                aGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCtata                                                                                                                                                                                                                    <  7:1120953/103‑1 (MQ=255)
                                                                                                                                                                                                                                aGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTg                                                                                                                                                                       <  8:580464/148‑1 (MQ=255)
                                                                                                                                                                                                                                 gAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTaa                                       >  4:37736/1‑275 (MQ=255)
                                                                                                                                                                                                                                     aCCGGGAACATACTTTTGACCCTTTTGTGTGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGCCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTaaa                                      <  4:400909/272‑1 (MQ=255)
                                                                                                                                                                                                                                       cgcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGac                                                                                                                                                                                   >  3:81482/1‑129 (MQ=255)
                                                                                                                                                                                                                                       cgcgAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGt                                 <  2:443371/275‑1 (MQ=255)
                                                                                                                                                                                                                                           aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                                                                                                                                                                                                     >  8:626478/1‑107 (MQ=255)
                                                                                                                                                                                                                                           aaCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGc                                                                                                                                                                                                     <  7:626478/107‑1 (MQ=255)
                                                                                                                                                                                                                                            aCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt                                                                                                                                                                                             >  8:211659/1‑114 (MQ=255)
                                                                                                                                                                                                                                            aCATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGt                                                                                                                                                                                             <  7:211659/114‑1 (MQ=255)
                                                                                                                                                                                                                                             cATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCta                                                                                                                                                                                                                      <  6:611532/88‑1 (MQ=255)
                                                                                                                                                                                                                                               tatCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtc                                                                                                                                                        >  5:463683/1‑148 (MQ=255)
                                                                                                                                                                                                                                                atCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAtt                                                                                                                                                                                                                                <  6:729558/75‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGTCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg                                                                                                                                                                                     >  8:498026/1‑117 (MQ=255)
                                                                                                                                                                                                                                                 tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt                                                                                                                                                                                         >  8:630086/1‑113 (MQ=255)
                                                                                                                                                                                                                                                 tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAAt                                                                                                                                                                                         <  7:630086/113‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtgtg                                                                                                                                                                                     <  7:498026/117‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCg                                                                                                                                                     <  8:676093/148‑1 (MQ=255)
                                                                                                                                                                                                                                                        aCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtct                                                                                                                                                       <  6:826618/140‑1 (MQ=255)
                                                                                                                                                                                                                                                        aCCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTtct                                                                                                                                                       >  5:826618/1‑140 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                                                                                                                                                                                            <  5:1039110/70‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                                                                                                                                                                                            >  6:1039110/1‑70 (MQ=255)
                                                                                                                                                                                                                                                          cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAAt                                                                                                                                                                                                                                 >  7:601800/1‑64 (MQ=255)
                                                                                                                                                                                                                                                          cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTgg                                                                                                                                                                      >  6:1070508/1‑123 (MQ=255)
                                                                                                                                                                                                                                                          cTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTg                                                                                                                                             >  5:617560/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           ttttGTTCGCGTGGTACCCGGACGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGTCAATTGTCCACTCTATAGCACTGGCAATGGGTGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTTAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTATGGCGCGTTACg             <  2:110816/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                                                                                                                                                                                            >  5:175636/1‑62 (MQ=255)
                                                                                                                                                                                                                                                                 tCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTcc                                                                                                                                                                                                                            <  6:175636/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtg                                                                                                                                                                                       >  6:998402/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                 tCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATtg                                                                                                                                                                                       <  5:998402/99‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGttt                                                                                                                                                  >  7:379377/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                    ccTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGttt                                                                                                                                                  <  8:379377/133‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCa                                                                                                                                 >  5:241272/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                      ccggtaccAGCCGGGCTCGCTCAACCGGCGGTTGCGGGCTGGTGCTGGACATGGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTATGGCGCGTTACGATCTCGTAGAc  <  4:175240/269‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCgg                                                                                                                                                                                                   >  6:316691/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                       cGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCgg                                                                                                                                                                                                   <  5:316691/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt                                                                                                                                     <  8:685926/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGATCCCGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTAt                                                                                                                                     >  7:685926/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                       cGATCACGGCTGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCgg                                                                                                                                                                                                   >  8:897288/1‑81 (MQ=255)
                                                                                                                                                                                                                                                                        gATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGtt                                                                                                                                       <  8:522229/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa                                                                                                                                                                                    <  6:131434/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGa                                                                                                                                                                                    >  5:131434/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                             cAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATc                                                                                                                          <  6:241272/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAg                                                                                                                                                                              <  8:526824/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAg                                                                                                                                                                              >  7:526824/1‑94 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                       
CATCTGGCAGATATTCAGGCAGTAACGCCCGATAAAACGGTACGGATAAAAACGCTCGTGCAAGGCCATTATGCGGCGTTGGATATGCGCTTAATGGAGCGCGTGCTGGATAATTTGCTCAATAACGCCCTGCGGTACTGCCATTCAACGGTTGAAACCAGCCTGCTACTGTCGGGGAATAGAGCGACATTAATTGTTGAGGATGATGGCCCAGGGATTGCCCCAGAAAACCGCGAACATATCTTTGAACCTTTTGTTCGCCTCGATCCCAGCCGGGATCGCTCAACCGGCGGCTGCGGGCTGGGGCTGGCAATTGTCCACTCTATAGCACTGGCAATGGGCGGTACGGTTAATTGTGACACCAGCGAACTGGGTGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTATGGCGCGTTACGATCTCGTAGAC  >  NC_000913/1683746‑1684282

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