Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,077 G→A 5.0% intergenic (‑113/+15) mfd ← / ← ycfT transcription‑repair coupling factor/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,0770GA 5.0% 1040.4 / 27.0 313intergenic (‑113/+15)mfd/ycfTtranscription‑repair coupling factor/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (160/137);  new base A (9/7);  total (169/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

Only 100 of 318 total aligned reads displayed.

CGCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAA  >  NC_000913/1173804‑1174348
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
cGCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCgg                                                                                                                                                                                                                                                                                 >  2:41131/1‑275 (MQ=255)
                 atCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAg                                                                                                                                                                                                                                                                >  2:606747/1‑275 (MQ=255)
                                                          aaattaccgcgcgcAGCGTCGAAAAGGCTCGCCGGGTTAACTTGCCAAGACGCTGCTGTCCACCCGCTTGGATGGGCAGCGAATAACGATATTGTTCAGGGATAACGGTTGTCAGATTCTCGTTAGGATATGCCTAAACATATGAGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAAGCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaa                                                                                                                                                                                                                       <  4:113524/262‑1 (MQ=255)
                                                                                                                cGCTGCTCACACGTGTGGACGGGCGGCTGATAAGGATCTGGTTCAGGCAGAACGGTTGTCAGATTGTCGTTAGGAGATGCATCAACATAGGGGGGCATATCTCTGATAGGCAATGTTGTTCATTATCCGTTGTCGTCAGGGTTTAGCAAAGTGAATCGCCCCGGGTGGGGCGATGGCTTAGCGGGAGTAACGGACGGAAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGaaa                                                                                                                                                                  <  4:276802/275‑1 (MQ=255)
                                                                                                                                   cGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                                                                                                                                             <  5:201255/149‑1 (MQ=255)
                                                                                                                                    gggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATCGTTTAGCAAATTGAATCGCCCCGGCt                                                                                                                                                                                                                                                                               >  8:759747/1‑146 (MQ=255)
                                                                                                                                      gCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGc                                                                                                                                                                                                                                                                          >  8:30225/1‑149 (MQ=255)
                                                                                                                                        aGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAc                                                                                                                                                                                >  2:529285/1‑237 (MQ=255)
                                                                                                                                           gTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                                                                                                                                                     <  5:1369095/149‑1 (MQ=255)
                                                                                                                                               aaCGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgggg                                                                                                                                                                                                                                                                           >  7:782636/1‑139 (MQ=255)
                                                                                                                                               aaCGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTa                                                                                                                                                                                                                                                                 <  5:609622/149‑1 (MQ=255)
                                                                                                                                                 cGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGc                                                                                                                                                                                                                                                               <  8:608886/149‑1 (MQ=255)
                                                                                                                                                  gATATTGGTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATGTGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACGTGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCg                                                                                                                                                                                                                                                              >  8:234460/1‑149 (MQ=255)
                                                                                                                                                   atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCgg                                                                                                                                                                                                                                                             <  7:234460/149‑1 (MQ=255)
                                                                                                                                                     atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGa                                                                                                                                                                                                                                                           >  5:1329121/1‑149 (MQ=255)
                                                                                                                                                     atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGa                                                                                                                                                                                                                                                           <  7:1224407/149‑1 (MQ=255)
                                                                                                                                                        gTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTA                                                                                                                                                                                                                                                        <  5:980575‑M1/149‑1 (MQ=255)
                                                                                                                                                         ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAA                                                                                                                                                                                                                                                       >  5:477677‑M1/1‑149 (MQ=255)
                                                                                                                                                           cAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                                                                                                                                                     <  5:946156/133‑1 (MQ=255)
                                                                                                                                                                aTAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACgg                                                                                                                                                                                                                                                >  7:594566/1‑149 (MQ=255)
                                                                                                                                                                  aaCGGTTGTCAGAATCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCCTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCTGGGGCGATGGCTTAGCGGGAGTAA                                                                                                                                                                                                                                                       <  6:450546‑M1/140‑1 (MQ=255)
                                                                                                                                                                   aCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAg                                                                                                                                                                                                                                             <  8:1225935/149‑1 (MQ=255)
                                                                                                                                                                     ggTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCg                                                                                                                                                                                                                                           >  6:554890/1‑149 (MQ=255)
                                                                                                                                                                         gTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGc                                                                                                                                                                        >  2:334315/1‑212 (MQ=255)
                                                                                                                                                                          tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAg                                                                                                                                                                                                                                      >  6:8759/1‑149 (MQ=255)
                                                                                                                                                                          tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAg                                                                                                                                                                                                                                      >  7:977880/1‑149 (MQ=255)
                                                                                                                                                                           cAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGt                                                                                                                                                                                                                                     >  5:127483/1‑149 (MQ=255)
                                                                                                                                                                            aGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTg                                                                                                                                                                                                                                    >  5:1473617/1‑149 (MQ=255)
                                                                                                                                                                             gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTg                                                                                                                                                                                                                                                                              >  8:1401127/1‑106 (MQ=255)
                                                                                                                                                                             gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGc                                                                                                                                                                    <  4:334311/212‑1 (MQ=255)
                                                                                                                                                                              aTTCTCGGTAGGATATGCCTCAACATATGGGGGATTATTACTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTGAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTgcggc                                                                                                    <  4:606743/275‑1 (MQ=255)
                                                                                                                                                                               ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                                                                                                                                                 >  8:667227/1‑149 (MQ=255)
                                                                                                                                                                               ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTCTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAg                                                                                                                                                                                                                                 <  5:295554/149‑1 (MQ=255)
                                                                                                                                                                                    gTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            >  8:1467118/1‑149 (MQ=255)
                                                                                                                                                                                    gTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAACCCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAg                                                                                              >  2:175774/1‑275 (MQ=255)
                                                                                                                                                                                    gTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGACGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            >  6:1136340/1‑149 (MQ=255)
                                                                                                                                                                                     ttAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgggg                                                                                                                                                                                                                           <  5:1170919/149‑1 (MQ=255)
                                                                                                                                                                                            atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                                                                                                                                                          <  6:365036/143‑1 (MQ=255)
                                                                                                                                                                                            atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaa                                                                                                                                                                                                                    <  7:1055633/149‑1 (MQ=255)
                                                                                                                                                                                               ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAg                                                                                                                                                                                                                 >  5:1721674/1‑149 (MQ=255)
                                                                                                                                                                                               ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAg                                                                                                                                                                                                                 >  7:1375312/1‑149 (MQ=255)
                                                                                                                                                                                               ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAg                                                                                                                                                                                                                 >  6:1692021/1‑149 (MQ=255)
                                                                                                                                                                                               ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAg                                                                                                                                                                                                                 >  5:1610646/1‑149 (MQ=255)
                                                                                                                                                                                                cTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGa                                                                                                                                                                                                                >  8:1036801/1‑149 (MQ=255)
                                                                                                                                                                                                 tCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                                                                                                                                               <  7:667227/149‑1 (MQ=255)
                                                                                                                                                                                                   aaCATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                                                                                                                                                     <  7:1287919/93‑1 (MQ=255)
                                                                                                                                                                                                   aaCATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                                                                                                                                                     >  8:1287919/1‑93 (MQ=255)
                                                                                                                                                                                                     cATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCt                                                                                                                                                                                                           <  8:352592/149‑1 (MQ=255)
                                                                                                                                                                                                     cATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGAGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCCGGGAGAACAACAGATCGCt                                                                                                                                                                                                           <  8:1431967/149‑1 (MQ=255)
                                                                                                                                                                                                          gggggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATg                                                                                                                                                                                                      >  7:747758/1‑149 (MQ=255)
                                                                                                                                                                                                          gggggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATg                                                                                                                                                                                                      >  5:74256/1‑149 (MQ=255)
                                                                                                                                                                                                            gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            <  7:42695/125‑1 (MQ=255)
                                                                                                                                                                                                            gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            >  8:42695/1‑125 (MQ=255)
                                                                                                                                                                                                             ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCAc                                                                                                                                                                                                   >  7:1167200/1‑149 (MQ=255)
                                                                                                                                                                                                               cATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCaga                                                                   >  2:589812/1‑275 (MQ=255)
                                                                                                                                                                                                                     tctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGaa                                                                                                                                                                                           >  6:512624/1‑149 (MQ=255)
                                                                                                                                                                                                                      ctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAa                                                                                                                                                                                          >  8:206695‑M2/1‑148 (MQ=255)
                                                                                                                                                                                                                       tGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgggg                                                                                                                                                                                                                                                                           <  5:546819/67‑1 (MQ=255)
                                                                                                                                                                                                                          tAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGttt                                                                                                                                                                                      >  8:1007337/1‑149 (MQ=255)
                                                                                                                                                                                                                           aaGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                                                                                                                                                     >  6:1658551/1‑149 (MQ=255)
                                                                                                                                                                                                                           aaGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                                                                                                                                                     <  5:366804/149‑1 (MQ=255)
                                                                                                                                                                                                                            aGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTc                                                                                                                                                                                    >  8:1602091/1‑149 (MQ=255)
                                                                                                                                                                                                                             gCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCggg                                                                                                                                                                                                                                                                            <  5:1674730/60‑1 (MQ=255)
                                                                                                                                                                                                                               aaTGCTGCTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATAGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGTACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                                                                                                                                                 >  8:217275/1‑149 (MQ=255)
                                                                                                                                                                                                                                     gTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACTGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGCGAAAGTTTTCTTAc                                                                                                                                                                                <  8:974592/144‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            >  3:203167/1‑99 (MQ=255)
                                                                                                                                                                                                                                      ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            <  1:203167/99‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCaa                                                                                                                                                                          <  7:450630/149‑1 (MQ=255)
                                                                                                                                                                                                                                        cATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGa                                          <  2:322448/274‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                                                                                                                                                       <  7:1593274/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                                                                                                                                                       >  6:1307302/1‑149 (MQ=255)
                                                                                                                                                                                                                                          ttATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGccc                                                                                                                                                                      >  7:1254431/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aTCCGTAATGGACATGGTTTACCAAATTGATGGGCCCCGGCTGGGGCGATGGCTTGGCGGGAGTAACTGACGGCAGGGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCg                                     <  4:502607/275‑1 (MQ=255)
                                                                                                                                                                                                                                             tCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCaa                                                                                                                                                                          >  5:431072/1‑142 (MQ=255)
                                                                                                                                                                                                                                              ccGTAATCGACATGTTTTAGCAAAGTGAATCGCCCCGGGTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATggg                                <  1:103131/278‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                                                                                                                                                <  7:179468/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTTAATATTTCAACCAGAATGcattgg                                 >  2:493913/1‑270 (MQ=255)
                                                                                                                                                                                                                                                gTAATCGAAATGGTTTAGCAAATTGAATAGCCCCGGCTGGGGCGCTGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAAATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCATCAACACCGCAACAGTCCAGATGAACAAACTGACAATGGGAAACAACAGCCGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCCTTTGCGCCGACCGAGTTAATGTGAATATTACAACCAGAATGctcatg                                 >  4:437491/1‑270 (MQ=255)
                                                                                                                                                                                                                                                 tAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAg                                                                                                                                                               <  5:869442/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                                                                                                                                             >  5:1012655/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aCATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                                                                                                                                                          >  8:941344/1‑84 (MQ=255)
                                                                                                                                                                                                                                                       aCATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                                                                                                                                                          <  7:941344/84‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCCGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCACAACAGTACAGa                                                                                                                                                         <  8:1058106/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aCATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCAGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                                                                                                                                                <  7:342737/142‑1 (MQ=255)
                                                                                                                                                                                                                                                        cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                            <  6:142768/81‑1 (MQ=255)
                                                                                                                                                                                                                                                        cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGAt                                                                                                                                                        <  7:1449043/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGAt                                                                                                                                                        <  6:954667/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTGGTTTAGCAAATTGACTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAAAAACAGATCGATGAATGCACGCGGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAATTGACGAAGGGATAAACCAACAGGAGTGTCAGTTAGACTTGCGTCGACAAGCGTGCTGCGTCCATTTGCGCAAGCAAAGGTAAGCTGAATGTCTCAACCCACATGCGCTGGGTGGAGTAg                        >  4:227471/1‑275 (MQ=255)
                                                                                                                                                                                                                                                            gTTTAGCAAATTGAATCGACCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATAGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCaa                 <  1:281138/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                gtaaaagtaaaaggCCCCGGATTGGGGGAGGGTTTGGGGGGAGTACCTGACGGAAGCGGGCAGGGAAGGCGGGGAGAACAACAGATCGATGAATGAACGCTGTGAAAGTTTTCTTACAAGCAAGCACGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCaa                 <  4:589808/262‑1 (MQ=255)
                                                                                                                                                                                                                                                                aGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGATTAACTGATGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGTTCGCTGAATGCACGCTGTGAAAGTTTTCTAACCAGCAAGCCCGCCACCGTACAGATTAACAAACTGACAAAGGGAAAAACCAGCATGCGTGTCCGTTCGACTTACGGCTAACACCGTGCCACGTGCCTTTGCGCAACATGAGTCAATCGCATACTTTCACCCAGAATGCGCTGGGTGGAGAAGATAACaa                 >  4:561062/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                 gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACt                                                                                                                                               >  6:1028570/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTc                                                                                                                                                                                    <  6:250327/110‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACaa                                                                                                                                          >  6:179189/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACaaa                                                                                                                                         <  8:11779/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aTCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAg                                                                                                                                                          >  6:158102/1‑129 (MQ=255)
                                                                                                                                                                                                                                                                           tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa                                                                                                                                     >  6:1237450/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             gCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATa                                                                                                                                   >  5:905426/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                   >  6:929122/1‑68 (MQ=255)
                                                                                                                                                                                                                                                                                 cGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCaa  <  2:692151/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                
CGCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAA  >  NC_000913/1173804‑1174348

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