Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,375 A→C 8.2% I120S (ATC→AGC)  betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,3750AC 8.2% 402.2 / 24.6 127I120S (ATC→AGC) betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (60/55);  new base C (7/4);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 130 total aligned reads displayed.

CACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGG  >  NC_000913/328110‑328649
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
cACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCgg                                                                                                                                                                                                                                                                           <  2:127046/275‑1 (MQ=255)
    accGGCAACCCGGTAAAAACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAGCGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgta                                                                                                                                                                                                                                                                       <  4:394671/274‑1 (MQ=255)
                         ccGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGCCGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCAc                                                                                                                                                                                                                                                   >  2:153958/1‑274 (MQ=255)
                           gTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTCTACCCAAACGCCCTTTCCCGCACAGGGATCTGGCCGCCTTCCAGCCCCGGGATCAGCCCCGCCTCATCCTCCCGCCCCTCCGCACCgg                                                                                                                                                                                                                                                >  4:334385/1‑275 (MQ=255)
                              cGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCg                                                                                                                                                                                                                                                                         <  2:219535/247‑1 (MQ=255)
                               ggacggCCCGCTGCGCTGTAAATTTCCNCCGGCTTTGCGGCGGTAAGCGGGGTNACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCACCGCAGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTANTTACATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTTCAGCACGTCCGCACCGGTa                                                                                                                                                                                                                                              <  3:86340/269‑1 (MQ=255)
                                                        cagcagccTTTAACGCGGGAACGGGGGTACATGCGCTCGGTCTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGAGCGGATTTCCACAGGGACACCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGTGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGa                                                                                                                                                                                                                    <  4:153957/267‑1 (MQ=255)
                                                                                   ttgtttggCTGGGTTTGAATTNNNNTTCGNGGGCGGCCGCCAGCGCCGGGGCGTGTGTCCACAGGGCACTCTGGATCGGGTAGTGCCATGCACCAAGCCCGGCCACTACGCCCAGCGGTTCGCGGCGTGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTcc                                                                                                                                                                                        <  3:589/268‑1 (MQ=255)
                                                                                                          aTTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCCCGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCGGGTGGCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTggc                                                                                                                                                                 >  4:200399/1‑273 (MQ=255)
                                                                                                          aTTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGc                                                                                                                                                                  >  2:183231/1‑274 (MQ=255)
                                                                                                              cGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGGCATTGCGTTCa                                                                                                                                                             >  2:461533/1‑275 (MQ=255)
                                                                                                                ttGCCGGCCGCCAGCGCCGGTGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACg                                                                                                                                                           <  2:533684/275‑1 (MQ=255)
                                                                                                                                  ggggCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgt                                                                                                                                                                                                                                                                        <  7:717104/148‑1 (MQ=255)
                                                                                                                                           ttCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCgg                                                                                                                                                                                                                                                                           >  5:1016240/1‑136 (MQ=255)
                                                                                                                                           ttCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCgg                                                                                                                                                                                                                                                                           <  6:1016240/136‑1 (MQ=255)
                                                                                                                                               acaGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGcc                                                                                                                                                                                                                                                                              >  8:14161/1‑129 (MQ=255)
                                                                                                                                               acaGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGcc                                                                                                                                                                                                                                                                              <  7:14161/129‑1 (MQ=255)
                                                                                                                                                   ggCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTCCCAGTTTTGCGAGTTCGTAATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAg                                                                                                                         >  2:470220/1‑274 (MQ=255)
                                                                                                                                                   ggCAAACTGGAGCGGGTAGTCCCATGGGCCAACCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAg                                                                                                                         <  4:183230/274‑1 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCgtagta                                                                                                                                                                                                                                                                    >  8:769571/1‑125 (MQ=255)
                                                                                                                                                             aTCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTCTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCTTTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATCCGCGAGCGCTCCATGGc                                                                                                              >  1:159392/1‑275 (MQ=255)
                                                                                                                                                                 ggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGa                                                                                                                                                                                                                                                                                    <  5:439519/105‑1 (MQ=255)
                                                                                                                                                                 ggTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGa                                                                                                                                                                                                                                                                                    >  6:439519/1‑105 (MQ=255)
                                                                                                                                                                   tAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCGGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       >  8:731707/1‑148 (MQ=255)
                                                                                                                                                                      ttCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTcc                                                                                                                                                                                                                                                       >  2:348311/1‑129 (MQ=255)
                                                                                                                                                                         cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCgg                                                                                                                                                                                                                                                                           >  8:695611/1‑106 (MQ=255)
                                                                                                                                                                         cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCgg                                                                                                                                                                                                                                                                           <  7:695611/106‑1 (MQ=255)
                                                                                                                                                                         cATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGCTGATCCCgg                                                                                                                                                                                                                                                                           >  6:393730/1‑106 (MQ=255)
                                                                                                                                                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCAc                                                                                                                                                                                                                                                   <  4:348310/129‑1 (MQ=255)
                                                                                                                                                                          aTGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGt                                                                                                                                                                                                                                >  8:1004309/1‑148 (MQ=255)
                                                                                                                                                                           tGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCCGAATATCAACGGCCCGACGCACAAGACGCGAGCGCTCAATGGCGGTCATCGACTcc                                                                                                 >  2:189437/1‑274 (MQ=255)
                                                                                                                                                                           tGCGACAATCCCTGCCCCGACGCCCCGCGGGTGGCGGCGGGGTTACCCCAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGGAGTACTCCAGCACGTCCGGACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTGTCCGGTGTCGAGGGTTTCCAGTGTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGcc                                                                                                 <  4:288745/274‑1 (MQ=255)
                                                                                                                                                                            gcgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGccc                                                                                                >  2:323361/1‑274 (MQ=255)
                                                                                                                                                                             cgcCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGa                                                                                                                                                                                                                             <  7:826913/148‑1 (MQ=255)
                                                                                                                                                                               ccAATCCCTGCCCCTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACg                                                                                                                                                                                                                                                          <  7:96922/117‑1 (MQ=255)
                                                                                                                                                                                  aTCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGttt                                                                                                                                                                                                                        <  7:1004309/148‑1 (MQ=255)
                                                                                                                                                                                       tGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGtt                                                                                                                                                               >  4:95972/1‑200 (MQ=255)
                                                                                                                                                                                           ccaccGCCCGGCGGGTCGTGGCGGGTTAACACAAAGGCCTTTTCACGCAACGGGATCTGGCTGCCTTCCCGCGCCGGGATCAGCCCGGCGTAGGACTCCAGCACGTCCGCACCGGTAACGATATCTACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGccc                                                                                   <  4:239951/268‑1 (MQ=255)
                                                                                                                                                                                           aCTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGGTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACg                                                                                                                                                           <  2:95973/200‑1 (MQ=255)
                                                                                                                                                                                            cTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGc                                                                                                                                                                                                              <  7:731707/148‑1 (MQ=255)
                                                                                                                                                                                               cGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                                                                                                                                                       <  5:105023/136‑1 (MQ=255)
                                                                                                                                                                                               cGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTc                                                                                                                                                                                                                       >  6:105023/1‑136 (MQ=255)
                                                                                                                                                                                                    aGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                                                                                                                                             >  5:762797/1‑109 (MQ=255)
                                                                                                                                                                                                    aGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTaa                                                                                                                                                                                                                                             <  6:762797/109‑1 (MQ=255)
                                                                                                                                                                                                     gCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTGGTACTCCAGCACGTCCGCACCGGTAACGata                                                                                                                                                                                                                                        <  8:188907/113‑1 (MQ=255)
                                                                                                                                                                                                     gCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTGGTACTCCAGCACGTCCGCACCGGTAACGata                                                                                                                                                                                                                                        >  7:188907/1‑113 (MQ=255)
                                                                                                                                                                                                        gTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                                                                                                                                                 <  7:20657/117‑1 (MQ=255)
                                                                                                                                                                                                        gTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACgg                                                                                                                                                                                                                                 >  8:20657/1‑117 (MQ=255)
                                                                                                                                                                                                          tCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       <  8:151739/109‑1 (MQ=255)
                                                                                                                                                                                                          tCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGatat                                                                                                                                                                                                                                       >  7:151739/1‑109 (MQ=255)
                                                                                                                                                                                                             cggcggGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         <  7:321041/104‑1 (MQ=255)
                                                                                                                                                                                                             cggcggGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         >  8:321041/1‑104 (MQ=255)
                                                                                                                                                                                                              gacgggtaaaacaaaaGGCCGTTTCACGCACCGGGACCTGGCTGCCTTCCAGCGCCTGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTACGGATATCGACGGTTGAGTTTTCCGAATATGCTTTTCCGGGGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGgcg                                                             <  4:239440/262‑1 (MQ=255)
                                                                                                                                                                                                                  ggTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCCCCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAACGCTCCATGGCGGGACACGCCGCCACGATTTTTTGCCCCTCTTGGGGCCGTTTCCCgggccga                                                          >  4:457557/1‑269 (MQ=255)
                                                                                                                                                                                                                     ataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         <  8:1075338/96‑1 (MQ=255)
                                                                                                                                                                                                                     ataCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGat                                                                                                                                                                                                                                         >  7:1075338/1‑96 (MQ=255)
                                                                                                                                                                                                                      taCACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGCATATCAACGGCCCGACGCAGAATACGCACGCGCTCCATGGCGGTCATCAACGCCCAACTTTTTTGCCCATTCTGGGCGATTTTCACGGCGCGATCGa                                                      >  2:237228/1‑274 (MQ=255)
                                                                                                                                                                                                                       acacAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTAGTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                                                                                                                                               >  7:229712/1‑136 (MQ=255)
                                                                                                                                                                                                                       acacAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTAGTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAg                                                                                                                                                                                               <  8:229712/136‑1 (MQ=255)
                                                                                                                                                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGa                                                                                                                                                                                                                                    >  8:166144/1‑98 (MQ=255)
                                                                                                                                                                                                                        cacaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGa                                                                                                                                                                                                                                    <  7:166144/98‑1 (MQ=255)
                                                                                                                                                                                                                         acaAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 <  8:901711/148‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACt                                                                                                                                                                                                                                                                  <  5:576271/66‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACt                                                                                                                                                                                                                                                                  >  6:576271/1‑66 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAg                                                                                                                                                                                                                                                              <  6:545243/70‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              <  5:554398/134‑1 (MQ=255)
                                                                                                                                                                                                                          caAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTc                                                                                                                                                                                         <  6:1095342/139‑1 (MQ=255)
                                                                                                                                                                                                                            aaGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGCTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTAATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCTGGATCGCAATCCtg                                               >  1:110755/1‑274 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCACCGGCCCGACGCAGAATACGCGAACGCTCCATGGCGGTCATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCCCTTTTCACGGCCCTACCCACATc                                                  >  4:155786/1‑272 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGATGGCTTCATGGCGGTGATCGACGCCCAGGTTTTTTGCCCCCGATGGTCGCTTGTCACGTCGTGATGGGCATCCTc                                               >  2:89078/1‑274 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCt                                                <  4:509344/273‑1 (MQ=255)
                                                                                                                                                                                                                             aGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCGCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGACCGCCATcccc                                               >  4:465194/1‑272 (MQ=255)
                                                                                                                                                                                                                                                     ggCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTNACGCAGAATATCAAc                                                                                                                                               >  3:26986/1‑154 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            >  6:604845/1‑91 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            >  8:618219/1‑91 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            <  5:604845/91‑1 (MQ=255)
                                                                                                                                                                                                                                                       cTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCtt                                                                                                                                                                                                            <  7:618219/91‑1 (MQ=255)
                                                                                                                                                                                                                                                         gCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGcc                                                                                                                                           <  1:26986/154‑1 (MQ=255)
                                                                                                                                                                                                                                                            ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTGTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACg                <  2:201762/274‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGt                                                                                                                                                                                            >  7:124038/1‑101 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGt                                                                                                                                                                                            <  8:124038/101‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                                                                                                                                             >  8:830062/1‑148 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                                                                                                                                              <  8:616173/146‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                                                                                                                                              >  7:616173/1‑146 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGc                                                                                                                                            <  7:830062/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              >  8:399940/2‑95 (MQ=255)
                                                                                                                                                                                                                                                                 acgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAgg                                                                                                                                                                                              <  7:399940/94‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 >  7:876817/2‑108 (MQ=255)
                                                                                                                                                                                                                                                                 acgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTg                                                                                                                                                                                 <  8:876817/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                 acgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     >  5:162621/2‑120 (MQ=255)
                                                                                                                                                                                                                                                                  cgcCGGGCTGATCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGttt                                                                                                                                                                                          >  7:116179/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                    ccGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                                                                                                                                                  >  6:79495/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                    ccGGCATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAat                                                                                                                                                                                                                  <  5:79495/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGa                                                                                                                                                                                                >  5:794005/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                       gggATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGa                                                                                                                                                                                                <  6:794005/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                                                                                                                                                      >  5:293060/1‑80 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCgg                                                                                                                                                                                                      <  6:293060/80‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     <  5:602856/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCAt                                                                                                                                                                     >  6:602856/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                          aTCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTCTCCGAATATGCTTTTCCGGTGTCCAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACTCAGAATCCGCGAGCGCTCCATGGCGGTCATCGCCGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGGTCCACCACCTCCGCCCCGTCCGCCGGCCCGGGCCCCAACACTTTCCCCGTGGcccg  >  4:522444/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                  
CACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGGCCCGCTTCGCTGTAAATTTCAGCCAGCTTTAACGCGGTAAGCGGGGTAACTTCGCTCGGTTTGAAAATCATTGCGTTGCCTGCCGCCAGCGCCGGGGCGGATTTCCACAGGGCAATCTGGATCGGGTAGTTCCATGCGCCAATCCCTGCCACTACGCCCAGCGGTTCGCGGCGGGTATACACAAAGGACGTTTCACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCACGTTACCGTTGGCCGG  >  NC_000913/328110‑328649

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