Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,061,916 C→T A356V (GCC→GTC)  scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,061,9160CT96.0% 234.2 / NA 76A356V (GCC→GTC) scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (40/31);  minor base G (0/3);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGA  >  NC_000913/3061670‑3062156
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
cGTTCTTCTTCGGCATCGGCATGTATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATTCCAACGTCtt                                                                                                                                                                                                                                                <  2:298308/250‑1 (MQ=255)
      tcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACg                                                                                                                                                                                                                                          <  1:360245/249‑1 (MQ=255)
      tcttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCAGCCCGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGATGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGCCTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACACCTTTGACGa                                                                                                                                                                                                                                         >  1:262399/1‑250 (MQ=255)
                 ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCATTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCGCAGGCTGGAGCCTTACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGGGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGttt                                                                                                                                                                                                                             <  1:219228/251‑1 (MQ=255)
                      ggATCTGTTTATGACCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGCGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTTCCt                                                                                                                                                                                                                         >  2:286899/1‑250 (MQ=255)
                      ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCt                                                                                                                                                                                                                         >  1:258462/1‑250 (MQ=255)
                       gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTAc                                                                                                                                                                                                                       >  1:172473/1‑251 (MQ=255)
                           tGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGAt                                                                                                                                                                                                                   >  1:106071/1‑251 (MQ=255)
                                 tGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGGAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTtctc                                                                                                                                                                                                              >  2:335733/1‑250 (MQ=255)
                                  gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGTCGAAGCGCTTGGTTTGCCTACCGATTTCTCAg                                                                                                                                                                                                            <  1:90230/251‑1 (MQ=255)
                                   aaCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCGCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGc                                                                                                                                                                                                           <  2:98994/251‑1 (MQ=255)
                                     aggcgcaaTGTTGCGAGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGCTTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGCAGCGCTTGGTTTGACGACCGATTTCTCAGccc                                                                                                                                                                                                         >  2:233945/8‑249 (MQ=255)
                                      gTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAAGCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTTCCTACCGATTTCGCAGCACg                                                                                                                                                                                                        <  2:276925/251‑1 (MQ=255)
                                        cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGACGCGCTTGGTTTGCCTACCGATTTCTCAGCACGc                                                                                                                                                                                                       >  2:217875/1‑250 (MQ=255)
                                         gCCATGTCGGGTGCGGCAAGTTATTTATGGAGCGCAGCGGTCAGTGGGTTTGGCGCACAGGACCCGAACTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGCACAGGATCCGTATAACAAGGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGCCGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAt                                                                                                                                                                                                     >  2:281449/1‑251 (MQ=255)
                                          ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAtt                                                                                                                                                                                                    >  1:43960/1‑251 (MQ=255)
                                            aTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGATCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGACAGACCTCCGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGCCGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCCGCACGCATTGc                                                                                                                                                                                                  >  2:86227/1‑251 (MQ=255)
                                              gTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGGGGAGTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGGCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGAGTGCGCAGCACGCATTGccc                                                                                                                                                                                                <  2:189045/251‑1 (MQ=255)
                                               ttGCGGGCGGCACGTTATTTATGAAGCGAAGCGGTCAGTGGTTTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACAACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCATCACGCATTGCCCg                                                                                                                                                                                               >  2:36369/1‑251 (MQ=255)
                                                       ggCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTAAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACAccc                                                                                                                                                                                       >  2:116434/1‑251 (MQ=255)
                                                          aCGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGGCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACTCCCAGa                                                                                                                                                                                    <  2:408509/251‑1 (MQ=255)
                                                           cGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGat                                                                                                                                                                                   <  1:96252/251‑1 (MQ=255)
                                                                 ttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCGCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcat                                                                                                                                                                             >  1:11392/1‑251 (MQ=255)
                                                                 ttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGTAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGAAGAAGCGCTTGGTTTGCCTACCGAGTTCTCAGCACGCATTGCCCGCAACACCCAGatcatca                                                                                                                                                                              <  1:315588/250‑1 (MQ=255)
                                                                  tatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcat                                                                                                                                                                             <  1:408752/250‑1 (MQ=255)
                                                                  tatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatcatcat                                                                                                                                                                             <  2:225264/250‑1 (MQ=255)
                                                                   atGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATcc                                                                                                                                                                           >  1:367504/1‑251 (MQ=255)
                                                                      gAGCGAAGCGGTCAGTGGATTGGGCGCACAGGACCCGAAAGCACGGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGGACTCAGTCACTGCATACCAACGTCTGTGACGAAGCGCTTGGTTTGCTTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAgg                                                                                                                                                                        <  1:281449/251‑1 (MQ=255)
                                                                         cGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCAAGCATTGCCCGCAACACCCAGATCATCCTCCAGga                                                                                                                                                                       >  2:357419/1‑249 (MQ=255)
                                                                               ggTCAGTGGATTTGGCTCACCGGACACGAAATCACTGCCGCTGCGTACCCACTGCCCGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTCTAACAACGATATCCGCCCCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGCATTCTCAGCCCGCATTGCCCGCAACACCCAGATCATCATCCAGGCAGAATc                                                                                                                                                                 >  2:361252/1‑249 (MQ=255)
                                                                               ggTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATc                                                                                                                                                                 >  2:278125/1‑249 (MQ=255)
                                                                                  cAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAAc                                                                                                                                                            <  2:143098/251‑1 (MQ=255)
                                                                                               gcACAGGACCCGAAATAACTGGCGCTGCGTACCCACTGCCCGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATCCCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCCTTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCCCCGt                                                                                                                                               >  2:153355/1‑251 (MQ=255)
                                                                                                 acaGGACACGAAATCACTGGCGCTGCGTACCAACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCCTGGATTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGACGAATCAGAACTCTGCCGCACCGTCg                                                                                                                                             >  1:88094/1‑251 (MQ=255)
                                                                                                         cGAAATCACTGGCGCTGCGTTCCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTg                                                                                                                                     <  2:79251/251‑1 (MQ=255)
                                                                                                         cGAAATCACTGGCGCTGCGTACCCACTGCCAGAACTCAGGCTGGAGCCTGACTGAACAGGATCCGTATCACAACGTTATCCGCACCACCATTGACGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACACAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTg                                                                                                                                     >  2:214682/1‑251 (MQ=255)
                                                                                                             aTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCg                                                                                                                                 <  2:26074/251‑1 (MQ=255)
                                                                                                               cACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCgg                                                                                                                                >  1:328989/1‑250 (MQ=255)
                                                                                                                aCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGAt                                                                                                                              >  2:41261/1‑251 (MQ=255)
                                                                                                                   ggtgctgccgaCCCACCGCCAGACCTAAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTTTCCGTACCACCATTTAAGCGATGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGGCTTTGACGAAGCGCTTGGGGTGCCTCACGATTTCTCAGCACGCATTGCCCGCAACCCCCAGATCATCATCCAGGAAGAAGCAGAACTCTGCCTCAACGGCGATCCACTGGCCggttc                                                                                                                             <  2:230496/239‑4 (MQ=255)
                                                                                                                   ggCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATc                                                                                                                             <  2:106071/249‑1 (MQ=255)
                                                                                                                          cGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACAt                                                                                                                    >  2:1293/1‑251 (MQ=255)
                                                                                                                                 ccTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCAGCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCTGATCCTATTACAGTGAGTCg                                                                                                             <  2:164724/250‑1 (MQ=255)
                                                                                                                                       aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAAAGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGCACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGAcc                                                                                                       >  1:160162/1‑251 (MQ=255)
                                                                                                                                        gACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCg                                                                                                      <  2:7219/251‑1 (MQ=255)
                                                                                                                                          ccTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGAt                                                                                                    <  1:1293/251‑1 (MQ=255)
                                                                                                                                           cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATc                                                                                                   >  1:330649/1‑251 (MQ=255)
                                                                                                                                            ccAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCa                                                                                                  <  2:72754/250‑1 (MQ=255)
                                                                                                                                              aGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCaa                                                                                                 >  1:379523/1‑250 (MQ=255)
                                                                                                                                                gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATGGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATc                                                                                              <  1:286899/251‑1 (MQ=255)
                                                                                                                                                  tGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGt                                                                                            >  1:286568/1‑251 (MQ=255)
                                                                                                                                                    gAGCCTGACTGAACAGGATCCGTATAACAACGTTATCAGCACCACCATTGACGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCCGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTc                                                                                           >  1:54526/1‑250 (MQ=255)
                                                                                                                                                      gCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCaaa                                                                                        <  2:420491/251‑1 (MQ=255)
                                                                                                                                                          gactgaACAGGATCCGTATAACAACGTTATCCGCCCCACCAGTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAg                                                                                    <  2:172473/251‑1 (MQ=255)
                                                                                                                                                          gactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCTGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACa                                                                                      <  1:357419/249‑1 (MQ=255)
                                                                                                                                                          gactgaACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGATGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACa                                                                                      >  2:339101/1‑249 (MQ=255)
                                                                                                                                                          gacggaACAGGATCCGTATAACAACGTTATCCGCACCAACATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTGGACGAAGCGCTTGGGTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCGGACCGATCAGATCGTCAAACAAg                                                                                    <  2:288411/247‑1 (MQ=255)
                                                                                                                                                                 cAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCACCACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGc                                                                             >  1:410033/1‑251 (MQ=255)
                                                                                                                                                                        cGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATc                                                                      >  1:399035/1‑251 (MQ=255)
                                                                                                                                                                             aacaacGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACa                                                                 <  1:156495/251‑1 (MQ=255)
                                                                                                                                                                                 acGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGGCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATc                                                             <  1:294623/251‑1 (MQ=255)
                                                                                                                                                                                  cGTTATCCGCTCCTCCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCTCTGCATACCAACGTCTTTGACGAAGCGCTTGGTGTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcg                                                            <  2:293141/251‑1 (MQ=255)
                                                                                                                                                                                      aTCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcgacg                                                         >  1:330083/1‑250 (MQ=255)
                                                                                                                                                                                       tCCGCACCACCATTGAAGAGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGCCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATcgacac                                                        >  2:419075/1‑248 (MQ=255)
                                                                                                                                                                                        ccGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAg                                                      <  2:328989/251‑1 (MQ=255)
                                                                                                                                                                                         cGCACCACCCTTGAATCGCTTGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTGGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGc                                                     <  2:296300/251‑1 (MQ=255)
                                                                                                                                                                                             ccaccaTTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCggg                                                 >  2:92731/1‑250 (MQ=255)
                                                                                                                                                                                              caccaTTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTg                                                >  2:4559/1‑251 (MQ=255)
                                                                                                                                                                                                      aagcgatcgccgcgacgaccggCGGTACTCAGTCCCTGCATACCAACGTCTTTACAGCAGCGCTTGGTTAGCCTACCGATTTCTCCAAACGCATTACCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCCCTGGCCGGCTCCTATTACATTGAGTCGCTGACCGCTCAAATCGTCAAACAAGCCAGAGCTATTATCCACCAGATCGACGAAGCCGGTGGCATGGCg                                        >  2:288112/21‑251 (MQ=255)
                                                                                                                                                                                                      aaGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCg                                        <  2:366977/251‑1 (MQ=255)
                                                                                                                                                                                                            tggctgCGACGCTGGGCGGTAATCACTCACCGCATACCAACGTCTTTAACGTAGCGCTTGGCTTGCCTACCGACTTCTCAGCACGCCCTGCCCGCACCCCCCCGATCTTCATCCAGGAAGAATCAGAACTCTGCCGCACCGTACATCCACTGGCCGGCTCCTCTTACATTGAGTCGCTGACCGATCAAATAGTCAAACACGCCAGAGCTATTATCCAACAGATCGCGGAAGCCGGTGGCATGGCGAAAGCg                                  >  2:178157/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGGGCGGTACTCAGTCACTGCATACCAACGTCTTTGACGAAGCGCTTGGTATGCCTACAGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATACACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCAGGTGGCATgacgttgacggtgggggcgggc                      >  1:72716/1‑230 (MQ=255)
                                                                                                                                                                                                                                   cAGTCACTGCATACCAACGTCTTTGACTAAGCTCTTGGTTCGCCTACCGATTTCTCAGCACGCATTGCCCGCACCACCCAGATCATCATCCATGAAGAATCAGAACTCTGCCGCACCGTCGATCCTCTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAATCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAATCGATCTTAGCATGTCTGCCAAAACg           <  2:399534/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cTGCATACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAACAGGAATGAt     >  2:84322/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa  >  2:364107/1‑250 (MQ=255)
                                                                                                                                                                                                                                             aTACCAACGTCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTACCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGa  >  1:435611/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
CGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGCTATTATCCAACAGATCGACGAAGCCGGTGGCATGGCGAAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGA  >  NC_000913/3061670‑3062156

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