Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,704,306 A→T 6.2% intergenic (‑243/+29) rnc ← / ← lepB RNase III/leader peptidase (signal peptidase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,704,3060AT 6.2% 665.3 / 23.0 190intergenic (‑243/+29)rnc/lepBRNase III/leader peptidase (signal peptidase I)
Reads supporting (aligned to +/- strand):  ref base A (95/83);  new base T (3/9);  total (98/92)
Fisher's exact test for biased strand distribution p-value = 7.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

Only 100 of 195 total aligned reads displayed.

CCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATCC  >  NC_000913/2704040‑2704564
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                  
                                      taaatGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGAc                                                                                                                                                                                                                                                            <  4:240540/237‑1 (MQ=255)
                                              tCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTTAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                                                                                                                                                                                              >  2:65496/1‑275 (MQ=255)
                                                   agagacagacacGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                                                                                                                                                                                          <  2:146128/266‑1 (MQ=255)
                                                          cacaCGAAACAGACCTGAGATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCa                                                                                                                                                                                                   <  2:562084/274‑1 (MQ=255)
                                                                   cAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                                                                                                                                         <  4:65496/275‑1 (MQ=255)
                                                                                cATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTAGTCGCCATAACGGTGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAATTAAGCGCAGACCCGCCGGCCATTCGCCTTCTTGcgt                                                                                                                                                                            >  4:165965/1‑273 (MQ=255)
                                                                                    gTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTTTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGa                                                                                                                                                                                                                                                 <  4:21553/202‑1 (MQ=255)
                                                                                    gTTCTGCTAACTGCTTCGCTGCAGTTTAAGTTTACCGGTATATGGGGACCAACGCTGGTTCGTGTGCCGTGGCAAACTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                                                                                                                                         <  2:18715/274‑1 (MQ=255)
                                                                                         gCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACt                                                                                                                                                                                                                <  2:369703/230‑1 (MQ=255)
                                                                                          cTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCAATTCGCCTTCTTGCTTATCGAAGCta                                                                                                                                                                  >  1:53818/1‑274 (MQ=255)
                                                                                           tAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACg                                                                                                                                                                                                                                                  <  2:344954/194‑1 (MQ=255)
                                                                                            aaCTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTc                                                                                                                                                                                                <  2:43045/243‑1 (MQ=255)
                                                                                                              tcagtttACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACCACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGGGACTTAAGCACAGACCAGTCGGCCATTCGCCTTATTTGTTAACGAAGCCCAACCAAATTCCCGTTGccc                                                                                                                                               >  4:394535/1‑274 (MQ=255)
                                                                                                                      ccGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATa                                                                                                                                                                                                                                                                    >  5:174205/1‑149 (MQ=255)
                                                                                                                       cGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCCATACACAACGAAGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATaa                                                                                                                                                                                                                                                                   >  6:785395/1‑149 (MQ=255)
                                                                                                                          tataTGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACgg                                                                                                                                                                                                                                                                <  6:420964/149‑1 (MQ=255)
                                                                                                                           atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                                                                                                                               <  6:134185/149‑1 (MQ=255)
                                                                                                                           atatGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                                                                                                                               <  7:470552/149‑1 (MQ=255)
                                                                                                                                gggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACCCAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACaa                                                                                                                                                                                                                                                          <  8:707477/149‑1 (MQ=255)
                                                                                                                                 ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                                                                                                                                         <  6:740160/149‑1 (MQ=255)
                                                                                                                                 ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                                                                                                                                         <  7:228694/149‑1 (MQ=255)
                                                                                                                                 ggACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAAc                                                                                                                                                                                                                                                         <  6:944463/149‑1 (MQ=255)
                                                                                                                                  gACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACg                                                                                                                                                                                                                                                        >  7:584850/1‑149 (MQ=255)
                                                                                                                                   aCCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGt                                                                                                                                                                                                                                                       <  5:797897/149‑1 (MQ=255)
                                                                                                                                      aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGaa                                                                                                                                                                                                                                                    >  7:72857/1‑149 (MQ=255)
                                                                                                                                      aaCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAACGTGaa                                                                                                                                                                                                                                                    >  8:766843/1‑148 (MQ=255)
                                                                                                                                       aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                                                                                                                                   <  5:76180/149‑1 (MQ=255)
                                                                                                                                       aCGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGAGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGACACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAAc                                                                                                                                                                                                                                                   <  7:708283/149‑1 (MQ=255)
                                                                                                                                         gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                                                                                                                               <  6:60790/135‑1 (MQ=255)
                                                                                                                                         gCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGc                                                                                                                                                                                                                                                               >  5:60790/1‑135 (MQ=255)
                                                                                                                                           tGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAg                                                                                                                                                                                                                                               >  5:339146/1‑149 (MQ=255)
                                                                                                                                                       cGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAAt                                                                                                                                                                                                                                   <  7:1080938/149‑1 (MQ=255)
                                                                                                                                                        gTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATg                                                                                                                                                                                                                                  >  8:598216/1‑149 (MQ=255)
                                                                                                                                                         tGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATgg                                                                                                                                                                                                                                 >  7:650043/1‑149 (MQ=255)
                                                                                                                                                           gCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGTGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACACCGTGAACGAATATGGCTATTAATGGAt                                                                                                                                                                                                                               >  6:1064457/1‑149 (MQ=255)
                                                                                                                                                             aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTaa                                                                                                                                                                                                                                    <  2:397018/142‑1 (MQ=255)
                                                                                                                                                             aaCCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                                                                                                                                             >  5:966227/1‑149 (MQ=255)
                                                                                                                                                             aaCCTGGAGGTGCCAACCTTAAACTTCAGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                                                                                                                                                                             >  7:722193/1‑149 (MQ=255)
                                                                                                                                                                cTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgc                                                                                                                                                                                                                          >  5:361184/1‑149 (MQ=255)
                                                                                                                                                                  ggAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCa                                                                                                                                                                                                                        >  5:532711/1‑149 (MQ=255)
                                                                                                                                                                   gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                                                                                                                                       >  8:902627/1‑149 (MQ=255)
                                                                                                                                                                   gAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCaa                                                                                                                                                                                                                       >  8:207212/1‑149 (MQ=255)
                                                                                                                                                                      gTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                                                                                                                                    >  5:455160/1‑149 (MQ=255)
                                                                                                                                                                         ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                                                                                                                                 >  8:763339/1‑149 (MQ=255)
                                                                                                                                                                         ccAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGAc                                                                                                                                                                                                                 >  8:167201/1‑149 (MQ=255)
                                                                                                                                                                            aCCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTa                                                                                                                                                                                                              <  6:695036/149‑1 (MQ=255)
                                                                                                                                                                              cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                                                                                                                                            >  7:950283/1‑149 (MQ=255)
                                                                                                                                                                              cTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAg                                                                                                                                                                                                            >  8:277534/1‑149 (MQ=255)
                                                                                                                                                                                 aaaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAACGCCCCAGTAACGGCCGTCGCCGCTGTTGTCGCGGtttgcc                                                                           >  4:552271/1‑271 (MQ=255)
                                                                                                                                                                                  aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCGTTATGGCGACAAc                                                                                                                                                                                                                                                         >  5:777111/1‑100 (MQ=255)
                                                                                                                                                                                  aaCTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                                                                                                                                                                                        <  7:207212/149‑1 (MQ=255)
                                                                                                                                                                                    cTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                                                                                                                                                                                      <  8:584850/149‑1 (MQ=255)
                                                                                                                                                                                       cGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                                                                                                                                         >  2:697022/1‑175 (MQ=255)
                                                                                                                                                                                        gggggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCg                                                                                                                                                                                                                   >  4:194054/1‑132 (MQ=255)
                                                                                                                                                                                          gggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                                                                                                                                                                                    <  6:970074/129‑1 (MQ=255)
                                                                                                                                                                                          gggAATATGCTAAACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATg                                                                  <  2:242362/275‑1 (MQ=255)
                                                                                                                                                                                           ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAg                                                                                                                                                                     <  4:697022/175‑1 (MQ=255)
                                                                                                                                                                                           ggAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCcg                                                                                               <  4:365023/245‑1 (MQ=255)
                                                                                                                                                                                            gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACtt                                                                                                                                                                                                               <  2:194054/132‑1 (MQ=255)
                                                                                                                                                                                            gAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                                                                                                                                                                              <  6:361184/149‑1 (MQ=255)
                                                                                                                                                                                                        cacaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttc                                                                                                                                                                                  >  5:206154/1‑149 (MQ=255)
                                                                                                                                                                                                         acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                                                                                                                                 >  5:685361/1‑149 (MQ=255)
                                                                                                                                                                                                         acaACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                                                                                                                                                                                 >  5:233894/1‑149 (MQ=255)
                                                                                                                                                                                                           aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTg                                                                                                                                                                               >  8:334404/1‑149 (MQ=255)
                                                                                                                                                                                                           aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaa                                                                                                                    >  1:144542/1‑208 (MQ=255)
                                                                                                                                                                                                           aaCGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCggagg                                                 >  2:606677/1‑275 (MQ=255)
                                                                                                                                                                                                            acgacgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCCCGTTGTCGCCATAGCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCAATTCGCCTTCTTGc                                                                                                                                                                              <  6:156069/149‑1 (MQ=255)
                                                                                                                                                                                                                cgGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCACGACCGACCAGATTCGCTTCCGGCACAAAGCCCAAGTAACGGCTGTCAGCGCTGGTGCCGGGGTCTGCGCCCCACATTGTGTATGGTCCCGGAGGAACag                                            >  4:73025/1‑274 (MQ=255)
                                                                                                                                                                                                                 gggggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc                                                                                                                                                                         >  7:302688/1‑149 (MQ=255)
                                                                                                                                                                                                                    ggATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGTTGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa                                                                                                                                                                      <  5:945321/149‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCTCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGTTTATCGAATCTCATCCAGATAGCCGTTGCCCGACCGACCTGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTTGCCCTTCATGAAGTATTTTCCTGGCGGAACAATCTAGGTTGc                                  >  1:140943/1‑275 (MQ=255)
                                                                                                                                                                                                                                ccACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGAt                                                                                                                                                          <  7:1023942/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATa                                                                                                                                                         >  6:123153/1‑149 (MQ=255)
                                                                                                                                                                                                                                 accGGGAGATTTATCTCATAAATAATTCCCGTTGTCGCCATAACGGCAACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAgttg                            <  2:73025/272‑1 (MQ=255)
                                                                                                                                                                                                                                      gagaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                                                                                                                                    >  6:584725/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gagaTTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACCTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGt                                                                                                                                                    <  8:356018/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ttATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGccc                                                                                                                                               <  5:659817/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGgctg                <  2:223522/275‑1 (MQ=255)
                                                                                                                                                                                                                                              tctcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgac                                                                                                                                            >  7:612509/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                                                                                                                                          >  7:375774/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tcATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCccgaccg                                                                                                                                          >  5:524034/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                                                                                                                                       <  8:43876/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   taaataaTTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCgaccgacc                                                                                                                                       <  6:886255/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      ataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGa                                                                                                                                    >  5:654984/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       taatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGAt                                                                                                                                   <  5:1064457/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttacagTTGGCGCCAAATCGGCGAAAACGTGAACGATGATGGCTATTAATGGATGCCCCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGCTAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATcc  <  4:254303/270‑1 (MQ=255)
                                                                                                                                                                                                                                                          ttCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATcc  <  4:435918/275‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGc                                                                                                                                                                    >  5:289564/1‑110 (MQ=255)
                                                                                                                                                                                                                                                             aCGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCtt                                                                                                                             >  8:914075/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTc                                                                                                                            <  5:123153/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTCAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCCAAGCTCTTCCAGATAGCCGTTGCCCGACCGACCAGATTCCCTTc                                                                                                                            <  6:238973/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gTTGTCGCCATAACGGCGACAACGTGAACGAAGCTGACTATTAATGGATGCCACCAATGCGACTTAGGCGCAAACCGGTTGGCCATTCTTCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCTGATTCGCTTcc                                                                                                                           >  8:569782/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                                                                                                                                                                             >  1:174978/1‑83 (MQ=255)
                                                                                                                                                                                                                                                               gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc                                                                                                                           <  6:191739/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTcc                                                                                                                           <  6:524034/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGc                                                                                                                        >  7:893416/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   tCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                                                                                                                                                                                         <  3:174978/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                     gCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGcaca                                                                                                                     >  8:1047776/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      ccATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaa                                                                                                                    >  6:570357/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       cATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACaaa                                                                                                                   <  7:764022/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                  
CCGATTAATTACGATGGGGTTCATGCGATACCAATAAATAAATGCGTCAACAATTCAGCACACGAAACAGACCTGATATACATGGTTCTGCTAACTGCTTCGCTGCAGTTTCAGTTTACCGGTATATGGGGACCAACGCTGTTTCGTGTGCCGTGGCAACCTGGAGGTGCCAACCTTAAACTTCGGGGGAATATTCTATACACAACGACGGGGGATGTCGTTAGCCACGGGAGATTTATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATCCAGATAGCCGTTGCCCGACCGACCAGATTCGCTTCCGGCACAAAGCCCCAGTAACGGCTGTCCGCGCTGTTGTCGCGGTTGTCGCCCATCATGAAGTATTGTCCCGGAGGAACAATCCAGGTTGCCAGTTGTTGCCCTGGCTGCTGGTAATACATCC  >  NC_000913/2704040‑2704564

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