Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,375,484 A→T 5.4% L160H (CTT→CAT) ‡ menC ← O‑succinylbenzoyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,375,4840AT 5.4% 708.0 / 23.1 202L160H (CTT→CAT) ‡menCO‑succinylbenzoyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base A (83/108);  new base T (5/6);  total (88/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

Only 100 of 207 total aligned reads displayed.

CGCCTGTACCTGCTCGCGTACTTTTTCCAGACTGCCCGTGAGCGTGGGTTTGATAACTACCGCGCGCACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCC  >  NC_000913/2375210‑2375754
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              
      tACCTGCTCGCGTACTTTTTCCAGACTGCCCGTGAGCGTGGGTTTGATAACTACCGCGCGCACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCa                                                                                                                                                                                                                                                                          >  4:253906/1‑275 (MQ=255)
                                   ccggggtggTGGGTTTTATAACTACCGCCCGCACCCCCCCTGCCGCCACAAAGGCAAAATCCGGCTCGCGCGGGCTTTCAACCCAGGCAATGGCAAGTCCGGTTTCCCGGGCAAACGGTCGCGAATCATGGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAAc                                                                                                                                                                                                                                             <  4:211244/267‑1 (MQ=255)
                                   ccggggaggggggTTTGCAAATCACCGCGGGCCGGTCCGCTGCAGCCACAAAGGCAAACGCCGGCTCGCGCAGGGTTCAATCCCAGGGAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGGTGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGGGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAAc                                                                                                                                                                                                                                             <  4:117832/265‑1 (MQ=255)
                                            tGGGTTTGATAACTACCGCGCGCACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACCCGATGCGGTCGCGATAATCCGGGGTAACGTATTTTGCAAACTGCTGACCTTTCCGCGGTTTCCAGGCGCGATTTGCGTCAAAACCCAAATGCAGAACCCGAATTGCCCCCCGCAACAGAATcac                                                                                                                                                                                                                                     >  4:93189/1‑274 (MQ=255)
                                                         tACCGCGCGCACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTcgc                                                                                                                                                                                                                       >  2:97301/1‑275 (MQ=255)
                                                              cgcgcACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCAc                                                                                                                                                                                                                   >  2:24488/1‑274 (MQ=255)
                                                                  cACGCCCTCTTCAGCCACAAAGGCAAAAGCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGAtt                                                                                                                                                                                                                                        <  2:82597/249‑1 (MQ=255)
                                                                                              tCCGGGTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCa                                                                                                                                                                                   <  2:224497/274‑1 (MQ=255)
                                                                                                                      cAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGttt                                                                                                                                                       <  1:11399/278‑1 (MQ=255)
                                                                                                                              ggcaatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                                                                                                                                                                                                <  7:1608229/149‑1 (MQ=255)
                                                                                                                              ggcaatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCaa                                                                                                                                                                                                                                                                                >  5:992290/1‑149 (MQ=255)
                                                                                                                                caatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGa                                                                                                                                                                                                                                                                              <  8:401369/149‑1 (MQ=255)
                                                                                                                                caatgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGa                                                                                                                                                                                                                                                                              <  6:1369898/149‑1 (MQ=255)
                                                                                                                                  atgCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCCCTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                                                                                                                                               >  2:199369/1‑274 (MQ=255)
                                                                                                                                       ggTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATc                                                                                                                                                                                                     <  2:113515/215‑1 (MQ=255)
                                                                                                                                        gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtc                                                                                                                                                                                                                                                                      <  8:1115998/149‑2 (MQ=255)
                                                                                                                                        gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCTTGACGCAAAtc                                                                                                                                                                                                                                                                      <  6:711846/149‑2 (MQ=255)
                                                                                                                                        gTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATg                                                                                                                                                                                                                                                                      <  7:787225/149‑1 (MQ=255)
                                                                                                                                               gggCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATcc                                                                                                                                                                                                                                                               <  6:316758/149‑1 (MQ=255)
                                                                                                                                                ggCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCg                                                                                                                                                                                                                                                              <  6:223750/149‑1 (MQ=255)
                                                                                                                                                   aaaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGaa                                                                                                                                                                                                                                                           <  7:1360629/149‑1 (MQ=255)
                                                                                                                                                   aaaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAAc                                                                                                                                                                                                                                             <  1:29892/163‑1 (MQ=255)
                                                                                                                                                    aaCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAAt                                                                                                                                                                                                                                                          <  8:257362/149‑1 (MQ=255)
                                                                                                                                                     aCGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACa                                                                                                                                                                                                                                            >  2:267880/1‑162 (MQ=255)
                                                                                                                                                         tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGAtt                                                                                                                                                                                                                                        <  4:267877/162‑1 (MQ=255)
                                                                                                                                                         tCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGCTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                                                                                                                                                     >  8:531242/1‑149 (MQ=255)
                                                                                                                                                           gcgAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTcc                                                                                                                                                                                                                                                   >  5:1321837/1‑149 (MQ=255)
                                                                                                                                                             gAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                                                                                                                                                                                                 <  7:322122/149‑1 (MQ=255)
                                                                                                                                                             gAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAg                                                                                                                                                                                                                                                 <  7:1072632/149‑1 (MQ=255)
                                                                                                                                                              aaTCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCTGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGc                                                                                                                                                                                                                                                >  5:761109/1‑149 (MQ=255)
                                                                                                                                                                tcatcGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCaa                                                                                                                                                                                                                                              <  8:1620779/149‑1 (MQ=255)
                                                                                                                                                                  atcGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACa                                                                                                                                                                                                                                            >  6:145249/1‑149 (MQ=255)
                                                                                                                                                                              gCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCAccacca                                                                                                                                                                                                                                <  7:351697/149‑1 (MQ=255)
                                                                                                                                                                                aCGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATg                                                                                                                                                                                                                              >  6:527379/1‑149 (MQ=255)
                                                                                                                                                                                 cGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGc                                                                                                                                                                                                                             >  5:918606/1‑149 (MQ=255)
                                                                                                                                                                                     tctTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCTGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTc                                                                                                                                                                                                                         <  5:1177485/149‑1 (MQ=255)
                                                                                                                                                                                       ttCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTcgc                                                                                                                                                                                                                       <  5:145249/149‑1 (MQ=255)
                                                                                                                                                                                       ttCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTcgc                                                                                                                                                                                                                       >  5:892414/1‑149 (MQ=255)
                                                                                                                                                                                          gagaAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCa                                                                                                                                                                                                                    <  8:330355/149‑1 (MQ=255)
                                                                                                                                                                                                   gATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGt                                                                                                                                                                                                           <  8:312583/149‑1 (MQ=255)
                                                                                                                                                                                                       cgtgctcgATAATCCGGTTTAAAGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGGGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTc                                                                      <  1:17445/272‑1 (MQ=255)
                                                                                                                                                                                                         gTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCt                                                                                                                                                                                                                                                     >  6:47118/1‑101 (MQ=255)
                                                                                                                                                                                                         gTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGt                                                                        >  2:314890/1‑274 (MQ=255)
                                                                                                                                                                                                            gcgATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCa                                                                                                                                                                                                  <  5:526741/149‑1 (MQ=255)
                                                                                                                                                                                                                  aTCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGa                                                                                                                                                                                            <  8:73708/149‑1 (MQ=255)
                                                                                                                                                                                                                   tCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGc                                                                                                                                                                                                                                                       >  5:3987/1‑89 (MQ=255)
                                                                                                                                                                                                                    ccGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAcc                                                                                                                                                                                          >  5:579760/1‑149 (MQ=255)
                                                                                                                                                                                                                    ccGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAcc                                                                                                                                                                                          >  5:1507874/1‑149 (MQ=255)
                                                                                                                                                                                                                    ccGGGCTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAcc                                                                                                                                                                                          <  7:531242/149‑1 (MQ=255)
                                                                                                                                                                                                                       ggTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTc                                                                                                                                                                                       <  7:1576856/149‑1 (MQ=255)
                                                                                                                                                                                                                            aCGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTcaccac                                                                                                                                                                                                                                  >  5:1148287/1‑101 (MQ=255)
                                                                                                                                                                                                                               tATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGGTCCGGAATTGCCGCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTTACCTTCGCCACttt                                                                                                                                                                               >  6:429327/1‑149 (MQ=255)
                                                                                                                                                                                                                                aTTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTa                                                                                                                                                                                                          <  6:1554086/121‑1 (MQ=255)
                                                                                                                                                                                                                                aTTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTa                                                                                                                                                                                                          >  5:1554086/1‑121 (MQ=255)
                                                                                                                                                                                                                                 tttGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACtttat                                                                                                                                                                             >  7:1010450/1‑147 (MQ=255)
                                                                                                                                                                                                                                    ggCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACttt                                                                                                                                                                                              >  7:908222/1‑129 (MQ=255)
                                                                                                                                                                                                                                      cAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGAACCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCt                                                                                                                                                                        <  8:1425570/149‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTg                                                                                                                                                                       >  8:186789/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTg                                                                                                                                                                       <  7:1598381/149‑1 (MQ=255)
                                                                                                                                                                                                                                         aCTGCTGACCTTTCAGCGGTGTCCATGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGCAATTGCCTACAGCAACAGATTCACCACCATGCCGTGTCACACCGCTTCGTACACTCCCACTTTCACCTTCGCCACATTCCCGCCTGGCATACCTGCACCTCTGAGGATCTGATCGTCCGGATCGCCATCCCACAGCCCGGCCGCACGGTAGTTCGCTCCGTGCGGCAACGGATCTCGCCGTGCAGACAAGGCACACCTTACCCACAAGGcccag                                       >  4:323850/1‑272 (MQ=255)
                                                                                                                                                                                                                                         aCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCaa                                                                                                                                                           >  4:445167/1‑159 (MQ=255)
                                                                                                                                                                                                                                         aCTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTCGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGCTTGAGGATCAg                                                                                                                                              >  1:26161/1‑172 (MQ=255)
                                                                                                                                                                                                                                             ctgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGttt                                                                                                                                                       <  2:445171/159‑1 (MQ=255)
                                                                                                                                                                                                                                              tgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATc                                                                                                                                                                >  6:955970/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tgACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATc                                                                                                                                                                >  6:6635/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCt                                                                                                                                                               >  7:589030/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCt                                                                                                                                                               >  8:1398276/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aCCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTg                                                                                                                                                              >  7:885419/1‑149 (MQ=255)
                                                                                                                                                                                                                                                aCCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTg                                                                                                                                                              >  7:509701/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGcc                                                                                                                                                                                                                                                      >  8:439868/1‑59 (MQ=255)
                                                                                                                                                                                                                                                  cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGcc                                                                                                                                                                                                                                                      <  7:439868/59‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCa                                                                                                                                                            >  7:494164/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   tttCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGa                                                                                                                                                                                                                                          <  8:1598921/70‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCaa                                                                                                                                                           <  6:793177/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACTCCAAATGCCACGGAAGGCATCTGCgg                         >  1:206687/1‑277 (MQ=255)
                                                                                                                                                                                                                                                       aGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGttt                                                                                                                                                       <  5:241904/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAg                                                                                                                                                    <  5:1609713/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                                                                                                                                               <  5:642913/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               agtgctgcgAGTTGCTTCAAGCCGTAACTGCAGATCCGGAATACCCTCCATCACCGGACTGACCGCGATGCCGTCGCGCGCCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGGGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGc                 <  4:143416/269‑1 (MQ=255)
                                                                                                                                                                                                                                                                cAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTtctc                                                                                                                                                                            <  7:1533141/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                cAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTtctc                                                                                                                                                                            >  8:1533141/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                 aGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGa                                                                                                                                             <  8:268183/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cgATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCCAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACATCGGTGCCTCACGGTCATTGGCTGCTTTCGGCAACGTATCTGTCAGTTCTGCCTATGCACAGATTACGCCCAACGCCACGGGAGGCCTCTGCGGTAGCCCGctaatc            >  4:616733/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                       gATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGc                                                                                                                                                                                     >  5:504391/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                       gATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCCCCTTCCGTCGCGCACCGCTTCGGACAATCCCACTTTGACCTTCGCCACTTTCTCCCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg                                                                                                                                      >  7:81333/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCa                                                                                                                                                                                                  <  8:1018524/89‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCa                                                                                                                                                                                                  >  7:1018524/1‑89 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGAc                                                                                                                                                                                           >  7:261862/1‑96 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCACTTTTAAGGATCAGACCGCCCg                                                                                                                                      >  6:1291731/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                                                                                                                                               <  6:111873/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                                                                                                                                               >  5:111873/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                        aTTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCg                                                                                                                                      <  6:1059477/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tttGCGTCTTGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCa                                                                                                                                               >  7:1398327/1‑139 (MQ=255)
                                                                                                                                                                                                                                                                         tttGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGa                                                                                                                                                     >  6:1154687/1‑133 (MQ=255)
                                                                                                                                                                                                                                                                             cGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATTGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTgccagcc  <  2:606483/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGc                                                                                                                                <  6:884764/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGc                                                                                                                                >  7:1252269/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              gTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGc                                                                                                                                >  7:729136/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               tCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGcc                                                                                                                               <  5:709213/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               tCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGcc                                                                                                                               <  5:840108/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              
CGCCTGTACCTGCTCGCGTACTTTTTCCAGACTGCCCGTGAGCGTGGGTTTGATAACTACCGCGCGCACGCCCTCTTCAGCCACAAAGGCAAAATCCGGCTCGCGCAGGCTTTCATCCCAGGCAATGGCAATGCCGGTTTCACGGGCAAACGCTCGCGAATCATCGCGGGTTTTGCACGGCTCTTCGAGAAACGCGATGCGGTCGCGATAATCCGGGTTAACGTATTTGGCAAACTGCTGACCTTTCAGCGGTGTCCAGGCGCGATTTGCGTCAAGACGCAAATGCAGATCCGGAATTGCCTCCAGCAACAGATTCACCACCATGCCGTCGCGCACCGCTTCGTACAATCCCACTTTGACCTTCGCCACTTTCTCGCCTGGCATATCTGCAAGTTTGAGGATCAGATCGTCCGGATCGCCATTACACAGCGGTGCCGCACGGTAGTTGGCTGCTTGCGGCAACGTATCTGTCAGTTCTGCCAATGCACAGCTTACGCCAAAGGCCACGGAAGGCATCTGCGGTAGCTCGCAATCGCCTGCCAGCC  >  NC_000913/2375210‑2375754

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