Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,949 A→C 5.5% I64L (ATC→CTC)  dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9490AC 5.5% 1785.6 / 76.1 523I64L (ATC→CTC) dinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base A (279/212);  new base C (16/13);  total (295/228)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-05

Only 100 of 536 total aligned reads displayed.

ACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGG  >  NC_000913/3817682‑3818223
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
aCGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGa                                                                                                                                                                                                                                                                             >  4:265866/1‑275 (MQ=255)
                                                 tttGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATcct                                                                                                                                                                                                                             >  2:510126/1‑272 (MQ=255)
                                                       aTTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGat                                                                                                                                                                                                                                                                            <  4:64618/221‑1 (MQ=255)
                                                          gATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                                                                                                                                                                                                         >  2:202328/1‑253 (MQ=255)
                                                                          ggCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTCATGCAAACGGAGCAGAACAATGGTCAGCAAGACACCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCAAAAGAAGCCTGTGAGAACAGTGGTCACAAAAGTGATGACCATTTCGAGAAGATCCTCGAAATGGTCAAAATGGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAACCAACCGCTATACATGGTACCTgg                                                                                                                                                                                                   >  4:517808/1‑275 (MQ=255)
                                                                          ggCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATTACTACTTCGAGGAGATCCTCGATATGCTCAAAATTGGCTAAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCGCTTACCTgg                                                                                                                                                                                                   >  4:109723/1‑275 (MQ=255)
                                                                          ggCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAACGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                                                                                                                                                   >  2:6927/1‑275 (MQ=255)
                                                                                                           gaTAAAACTGATTAATGCAAACGGAGCAGAACCATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTAACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTAAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTATCCCGCTATGCCTGGTACATGGTAGTAAAAACCGGCGACCCTGGGAACCCGGTCa                                                                                                                                                                  >  4:321437/1‑275 (MQ=255)
                                                                                                           gaTAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCa                                                                                                                                                                                                                                                   >  2:554452/1‑194 (MQ=255)
                                                                                                                     aTTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggaagg                                                                                                                                                        >  4:218440/1‑275 (MQ=255)
                                                                                                                                   aGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCg                                                                                                                                                                                                                              <  2:23939/191‑1 (MQ=255)
                                                                                                                                     cAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGt                                                                                                                                                                                                                                                                       <  7:995037/148‑1 (MQ=255)
                                                                                                                                         aCAATGGTCAGCAAGACAACTTGGGACACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaaa                                                                                                                                                                                                                                                                   <  5:489737/148‑1 (MQ=255)
                                                                                                                                          aaaTGGTCAGCAAGAAAAGTTTGGACAATNCTTGGTTTTTCCGAGTTGAGTCCCTNTATACCTTTATTAACGCGCGACAAAGACGCCTGTGAGAACAGTGGTCACACAAGTGATGACCATTTCGAGGAGATCCTCGNTATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCa                                                                                                                                   <  3:94245/274‑1 (MQ=255)
                                                                                                                                                 cAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTa                                                                                                                                                                                                                            <  2:308863/179‑1 (MQ=255)
                                                                                                                                                  aGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATg                                                                                                                                                                                                               >  2:570298/1‑191 (MQ=255)
                                                                                                                                                      aGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                                                                                                                                               >  6:979885/1‑139 (MQ=255)
                                                                                                                                                           aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTc                                                                                                                                                                                                                                                                               <  8:1044376/118‑1 (MQ=255)
                                                                                                                                                           aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTc                                                                                                                                                                                                                                                                      >  8:637914/1‑127 (MQ=255)
                                                                                                                                                           aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACGTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTc                                                                                                                                                                                                                                                                               >  7:1044376/1‑118 (MQ=255)
                                                                                                                                                            cTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaaaa                                                                                                                                                                                                                                                >  7:247593/1‑148 (MQ=255)
                                                                                                                                                                       actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGa                                                                                                                                                                                                                                                                             <  5:328514/108‑1 (MQ=255)
                                                                                                                                                                        ctGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGa                                                                                                     <  2:175223/275‑1 (MQ=255)
                                                                                                                                                                         tGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGa                                                                                                                                                                                                                                   <  8:610825/148‑1 (MQ=255)
                                                                                                                                                                          gggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCg                                                                                                                                                                                                                              <  4:558313/152‑1 (MQ=255)
                                                                                                                                                                                   agagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGAc                                                                                                                                                                                                                                  <  5:125084/139‑1 (MQ=255)
                                                                                                                                                                                    gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                      >  8:1080127/1‑134 (MQ=255)
                                                                                                                                                                                     agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCATGGCAGACTTATTTTGCTATCCAGACCCTACGTCAGGCGCTTGCTTACGACGAAGCATTTACGCaa                                                                                        >  1:59241/1‑275 (MQ=255)
                                                                                                                                                                                     agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCaa                                                                                        >  2:47428/1‑275 (MQ=255)
                                                                                                                                                                                     agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTATCCAGCCATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAACCCGNTCATTGCGGCAGGGCAGACTAATTTTGCCATCTAGAACCGACGGACGGAGCTTGCTCACGATGAAGCATTTCACCaa                                                                                        >  3:122064/1‑275 (MQ=255)
                                                                                                                                                                                            tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                      <  6:122607/126‑1 (MQ=255)
                                                                                                                                                                                            tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTAt                                                                                                                                                                                                                >  8:592821/1‑148 (MQ=255)
                                                                                                                                                                                             cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCt                                                                                                                                                                                                                                                                                <  7:611463/83‑1 (MQ=255)
                                                                                                                                                                                              acctttaaCCTGTATTAACGCGCGCCAAAGAAGCCTGTTAGAACAGTGGTCACAGAATTTATGACCATTTCGAGGAGAGCCTCGATATGGTCAACATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCCGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGaa                                                                               <  4:147266/268‑1 (MQ=255)
                                                                                                                                                                                              aCTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCg                                                                                                                                                                                                                              <  4:349958/132‑1 (MQ=255)
                                                                                                                                                                                              aCTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                                                                                                                                                   <  2:138631/159‑1 (MQ=255)
                                                                                                                                                                                                 ttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCg                                                                                                                                                                                                                                                                              >  5:990663/1‑81 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGt                                                                                                                                                                                                                             <  7:53183/123‑1 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                    >  8:127555/1‑148 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                    >  7:599892/1‑148 (MQ=255)
                                                                                                                                                                                                        tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                    <  7:127555/148‑1 (MQ=255)
                                                                                                                                                                                                          tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGTAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTTAAGACGAAAAGCGt                                                                   >  2:390672/1‑275 (MQ=255)
                                                                                                                                                                                                          tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGAATTTAAGCAACTTCGTGAAGACGAAAAGCGt                                                                   >  2:357232/1‑275 (MQ=255)
                                                                                                                                                                                                           aTTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGa                                                                                                                                                                                                                                                                             <  8:528210/72‑1 (MQ=255)
                                                                                                                                                                                                              aaCGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACAt                                                                                                                                                                                                                                <  5:428429/114‑1 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTTCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                         >  4:479809/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGa                                                                                                                                                   >  2:213445/1‑185 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGCCTTATTTTGCTATCCAGACCCGACGGACGGAGCTTGCTGACGATGAAGCATTGAAGCAACTTCGTGTAGACGAAAAGCGCCTGTTTCGga                                                         >  4:363807/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAg                                                                                                                                  >  2:442022/1‑202 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCATGAGCTTGCTTACGGTGAAGCATTTAAGCCACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                         >  2:412441/1‑275 (MQ=255)
                                                                                                                                                                                                                    cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTga                                                         >  2:259283/1‑275 (MQ=255)
                                                                                                                                                                                                                       cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTcc                                                                                                                                                                                                                      <  6:403200/115‑1 (MQ=255)
                                                                                                                                                                                                                        aaaGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGAccc                                                                                                                              <  4:442021/202‑1 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAACTTGGCt                                                                                                                                                                                                                                                            <  8:716028/75‑1 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                                                                                                                                                              <  7:360757/73‑1 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                                                                                                                                                              >  8:360757/1‑73 (MQ=255)
                                                                                                                                                                                                                         aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                                                                                                                                                                                                            >  7:716028/1‑75 (MQ=255)
                                                                                                                                                                                                                         aaaaagCCTGTGGGAAATGGGGTGAAACGATTGATGACCATTTCGGAGAGATCCTCGTGTGGGTCAAAAATGGATCAAATGCCAAAAGAGCATTAAAAGGCAACGTACTCTCCCGCTATGCCTGTTACCCGGTAGGACAAAACGGCGACCCTGGGAAAGCGGTAATTTGGGCAGGGCAGACTTATTGTGCTATCCAGGCCCGACGGCAGGAGCTTGCGGACGATGAAGCATTGAAGCAACGTCGTTAAGACGAAAAGCGTCTGTTTCTGAGAAAt                                                    <  4:2977/272‑1 (MQ=255)
                                                                                                                                                                                                                          agaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAc                                                                                                                                                                                                       <  7:844112/127‑1 (MQ=255)
                                                                                                                                                                                                                                  gtgAAAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta                                                                                                                                                                                              >  6:31895/1‑128 (MQ=255)
                                                                                                                                                                                                                                   tgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGATTTGAAAg                                          >  1:94364/1‑275 (MQ=255)
                                                                                                                                                                                                                                    gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                                                                                                                                               <  7:906548/61‑1 (MQ=255)
                                                                                                                                                                                                                                    gagaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                                                                                                                                                                                                               >  8:906548/1‑61 (MQ=255)
                                                                                                                                                                                                                                     agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGa                                         <  2:278769/274‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                                                                                                                                                                                                    >  7:644985/1‑117 (MQ=255)
                                                                                                                                                                                                                                           gTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCCCAAATGCCAAAagagta                                                                                                                                                                                                                                           <  8:84251/74‑3 (MQ=255)
                                                                                                                                                                                                                                            tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATaaca                                 >  2:336115/1‑275 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                      <  7:756841/77‑1 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaa                                                                                                                                                                                                                                      >  8:756841/1‑77 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATaacaa                                >  4:415346/1‑275 (MQ=255)
                                                                                                                                                                                                                                             ggTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTAAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGCGCGTGCTGACGATGAAGCATTTAAGCAACTTCGTGACGAAGCAACGAGACTGTCTCTGAGAATAGAGTTAAACGAACATaaca                                 >  4:209799/1‑274 (MQ=255)
                                                                                                                                                                                                                                              gacacaaaATGGATGAACATTGCGAGGAGATCCTCGATAGAGTAAAAATTGGCACAAATGCCAAAAGAGCATTAAAAGACATCGTACTCGCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATaacaa                                <  1:153100/272‑1 (MQ=255)
                                                                                                                                                                                                                                               tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCg                                                                                                                                                             <  6:1097390/148‑1 (MQ=255)
                                                                                                                                                                                                                                                cacacaATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                                                                                                                                                                                                          <  6:698270/38‑1 (MQ=255)
                                                                                                                                                                                                                                                cacacaATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                                                                                                                                                                                                          >  5:698270/1‑38 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                                                                                                                                                                                                          >  6:541064/1‑100 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                                                                                                                                                                                                          <  5:541064/100‑1 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                                                                                                                                                                                                  >  5:60930/1‑108 (MQ=255)
                                                                                                                                                                                                                                                  cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggca                                                                                                                                                          <  8:247593/148‑1 (MQ=255)
                                                                                                                                                                                                                                                    caATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                                                                                                                                                     >  4:396053/1‑71 (MQ=255)
                                                                                                                                                                                                                                                      atggatgaCCATTTCGATGAGATCCACGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGGCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGt                       <  4:243538/272‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTTATGACCATTACGGGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGt                       <  4:195346/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGt                       >  2:366729/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           tgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                                                                                                                                                                                                             >  6:487395/1‑88 (MQ=255)
                                                                                                                                                                                                                                                            gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTa                                                                                                                                                                                                                            >  6:940039/1‑72 (MQ=255)
                                                                                                                                                                                                                                                            gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTa                                                                                                                                                                                                                            <  5:940039/72‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggca                                                                                                                                                          >  6:645180/1‑137 (MQ=255)
                                                                                                                                                                                                                                                             aCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTa                                                                                                                                                                                                        <  6:821558/91‑1 (MQ=38)
                                                                                                                                                                                                                                                              ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGaaa                                                                                                                                                                         >  5:1103867/1‑121 (MQ=255)
                                                                                                                                                                                                                                                              ccATTTCGAGAAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                                                                                                                                                                                                          >  5:836857/1‑88 (MQ=255)
                                                                                                                                                                                                                                                               cATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcag                                                                                                                                                         >  7:1106778/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                 tttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGNGGCAGGGCAGACTTATTTTGCTATCCAGACACGACGGCAGGAGCTTGCTGACGATGAAACATTTAAGCANATTCGTTAAGACGAAAAGCGTCGGTTTATGAGAAATGAGGTGAAANANACTAACAAACAAATGGTTAGGgcaaaaa            >  3:35035/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                  ttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTc                                                                                                                                                                   >  8:461636/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                  ttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTc                                                                                                                                                                   <  7:461636/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagt                                                                                                                                                                                               >  5:1010767/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                    cGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt                                                                                                                                             <  7:932731/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cGAGGATCTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt                                                                                                                                             >  8:932731/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                        gagaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                                                                                                                                                                                                  <  6:215451/86‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gagaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                                                                                                                                                                                                  >  5:215451/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                           aTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGAACAGCAAGCTgg  >  1:113756/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
ACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGACCCGACGGCAGGAGCTTGCTGACGATGAAGCATTTAAGCAACTTCGTGAAGACGAAAAGCGTCTGTTTCTGAGAAATGAGTTGAAAGAACATAACAAACAATTGGTTGAGGCAGCACAGCAAGCTGG  >  NC_000913/3817682‑3818223

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