Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,061,817 G→C 6.3% W323S (TGG→TCG) ‡ scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,061,8170GC 6.3% 643.4 / 24.8 188W323S (TGG→TCG) ‡scpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (79/96);  new base C (3/9);  total (83/105)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 189 total aligned reads displayed.

TACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGC  >  NC_000913/3061543‑3062081
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
                                          gCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcac                                                                                                                                                                                                                                <  2:195709/275‑1 (MQ=255)
                                                     cgctGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      >  4:246365/1‑274 (MQ=255)
                                                                  gAGTACATCAAAGCAGCAATCTCTGCCGGACGGACAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCCTCGGCATGGATCTGTTTATGAACGTCGCCCTGTTGCGTTCTGCACGTTATTTATGGAGCTAACCGGTCAGTGGATTTGTTGCACCGGACCCTAAATCACTGGCGCTGCGTACCCACTTCCCTACCCACAGCTGGAGCCTGACTTAACACGATCCCTCTAACAACGTTTCCTTCACTACCATTCAACCccttgctccca                                                                                                                                                                                                         >  4:419202/1‑270 (MQ=255)
                                                                             agcagcAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGc                                                                                                                                                                                                                                                                      >  4:552862/1‑202 (MQ=255)
                                                                                 gcaATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                                                                                                                                                                                                                                                  <  2:552862/202‑1 (MQ=255)
                                                                                         tGCCGGACTGAAAATTGAGTACTTCGCTCCTCGCCGGTCGTTCTTCTTCGGCATCGGCATGGATCGGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTGGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGc                                                                                                                                                                                 <  1:143253/275‑1 (MQ=255)
                                                                                              gACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCCGTCACTGCATAcc                                                                                                                                                                            >  2:101240/1‑275 (MQ=255)
                                                                                                          attccTTCGCCCCTCGCATGTCGTTCTTCTTCCGCTTCGTCATGGATCTGTTCAGGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGCGGATTTGGCGAACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACCACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGAc                                                                                                                                                                <  1:24958/271‑1 (MQ=255)
                                                                                                                         gCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTg                                                                                                                                                 >  2:321646/1‑275 (MQ=255)
                                                                                                                             gTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTa                                                                                                                                             <  2:246365/275‑1 (MQ=255)
                                                                                                                              tCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGCACGTCGCCATGTTGCGTGCGGCCCGTTACTTATGGAGCGAAGCGGTCAGTGGATTTCGCGCACAGGACCCGAAATCACCTGCTCTGCGTACCCACTGCCAGACCTCAGGCTTGCGCCTGACTGACCAGGATCCCTATAACAACGTTATCCGCCCCACCATTGAATCGCTGGCTGCGCCTCTGGGCGGTCCCCAGTCACTGCCTACCAACGCCTTCGACGAGGCGCGTGGTTCGCCTa                                                                                                                                             >  4:136349/1‑274 (MQ=255)
                                                                                                                                gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                                                                                                                                                                                                                                                        <  5:443317/149‑1 (MQ=255)
                                                                                                                                gTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                                                                                                                                                                                                                                                        <  7:262068/149‑1 (MQ=255)
                                                                                                                                      cttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTTCCCACTGCCAGACCTCAGGCTGGAGCctga                                                                                                                                                                                                                                                                  <  5:604136/149‑1 (MQ=255)
                                                                                                                                       ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgac                                                                                                                                                                                                                                                                 <  7:920672/149‑1 (MQ=255)
                                                                                                                                       ttcGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 <  7:119763/149‑2 (MQ=255)
                                                                                                                                        tcggcatcggcaTGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCTAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgagg                                                                                                                                                                                                                                                                <  6:943851/149‑3 (MQ=255)
                                                                                                                                              tcggcatGGATCTGTTTATGAACGTTGCCATGCTGCGTGCTGCACGTTATTTATGGAGCGAAGCGGTCCGTGGATTTGGCGCACAGGACCCGAGATCACTCGTGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACCGAACAg                                                                                                                                                                                                                                                          <  6:1089897/149‑1 (MQ=255)
                                                                                                                                                ggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGAACCCGAAATCACTGGCGATGCGTACCCACTGCCAGACCTCGGGCTGGAGCCTGACTGAACAGGa                                                                                                                                                                                                                                                        >  8:643937/1‑149 (MQ=255)
                                                                                                                                                 gcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       >  7:831902/1‑149 (MQ=255)
                                                                                                                                                   atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                                                                                                                                                                                                                                                     >  7:218629/1‑149 (MQ=255)
                                                                                                                                                   atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATcc                                                                                                                                                                                                                                                     <  7:547801/149‑1 (MQ=255)
                                                                                                                                                   atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgag                                                                                                                                                                                                                                                                 >  8:106095/1‑136 (MQ=255)
                                                                                                                                                    tGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                                                                                                                                                                                                                                                    <  8:273426/149‑1 (MQ=255)
                                                                                                                                                      gATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGta                                                                                                                                                                                                                                                  >  8:210858/1‑149 (MQ=255)
                                                                                                                                                       aTCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                                                                                                                                                 >  7:481916/1‑149 (MQ=255)
                                                                                                                                                           gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCCCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                                                                                                                                             >  8:13939/1‑149 (MQ=255)
                                                                                                                                                                tGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGTCCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCa                                                                                                                                                                                     <  2:45488/200‑1 (MQ=255)
                                                                                                                                                                   aCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATc                                                                                                                                                                                                                                     >  7:690034/1‑149 (MQ=255)
                                                                                                                                                                    cGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATcc                                                                                                                                                                                                                                    >  6:837635/1‑149 (MQ=255)
                                                                                                                                                                     ggcGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGAGTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGGCCTCAGGCTGGAGCCTGACTGAACAGAATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGatca                                                                                                      <  1:190615/272‑1 (MQ=255)
                                                                                                                                                                      tCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTTCCAGACCTCAGCCTGGAGCCTGACTGAACAGGATCCGtat                                                                                                                                                                                                                                                 >  1:62483/1‑134 (MQ=255)
                                                                                                                                                                      tCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCCCCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCACCGCCTTTGACGAAGAGCTTGGTTTGCATACCGATTTCCCATCACGCATTGCCCGCAACAccccaatcct                                                                                                     >  4:4840/1‑267 (MQ=255)
                                                                                                                                                                      tCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAg                                                                                                                                                                                        >  2:57946/1‑191 (MQ=255)
                                                                                                                                                                      tCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAc                                                                                                                                                                         >  4:659804/1‑206 (MQ=255)
                                                                                                                                                                      tCGCCATGTTGCGTGCGGCACGTTATTTACGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCCACGCTTGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACACAGCGTCTGGTTTGCCTACCGCATTCCCAGCCGGCATTGCCCCCACCACACAGatccat                                                                                                    >  4:230829/1‑272 (MQ=255)
                                                                                                                                                                        gCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCa                                                                                                                                                                           <  4:517080/202‑1 (MQ=255)
                                                                                                                                                                         ccATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacc                                                                                                                                                                                                                               <  7:577896/149‑1 (MQ=255)
                                                                                                                                                                          cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcacca                                                                                                                                                                                                                              <  5:982907/149‑1 (MQ=255)
                                                                                                                                                                          cATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTAAGTGGATTTGGAGGACAGGACCCGAAATCACTGGCGCAGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGAACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCa                                                                                                           <  2:230829/264‑1 (MQ=255)
                                                                                                                                                                                  gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  8:791838/149‑1 (MQ=255)
                                                                                                                                                                                  gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  7:202136/149‑1 (MQ=255)
                                                                                                                                                                                  gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  5:544209/149‑1 (MQ=255)
                                                                                                                                                                                  gtgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                                                                                                                                                      <  6:621705/149‑1 (MQ=255)
                                                                                                                                                                                   tgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACa                                                                                                                                                                                                                                                           >  2:589167/1‑111 (MQ=255)
                                                                                                                                                                                       gCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       <  4:589167/111‑1 (MQ=255)
                                                                                                                                                                                           gTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggct                                                                                                                                                                                                             >  8:1040257/1‑149 (MQ=255)
                                                                                                                                                                                            ttattaatGGAGCGAAGCGGTCAGTGGATTTTGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgctggctg                                                                                                                                                                                                            >  6:570624/1‑149 (MQ=255)
                                                                                                                                                                                              atttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCg                                                                                                                                                                                                          >  5:108083/1‑149 (MQ=255)
                                                                                                                                                                                                       gCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACt                                                                                                                                                                                           >  3:74672/1‑155 (MQ=255)
                                                                                                                                                                                                         gAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATcc                                                                                                                                                                                                                                    >  4:185591/1‑112 (MQ=255)
                                                                                                                                                                                                         gAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCgg                                                                                                                                                                                               <  7:402631/149‑1 (MQ=255)
                                                                                                                                                                                                            gCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAg                                                                                                                                                                                                                                                                       <  2:702339/74‑1 (MQ=255)
                                                                                                                                                                                                             cGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGcac                                                                                                                                                                                                                                <  2:185591/112‑1 (MQ=255)
                                                                                                                                                                                                                    tGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAg                                                                                                                                                                                        <  8:359848/145‑1 (MQ=255)
                                                                                                                                                                                                                     ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                                                                                                                   <  7:5004/149‑1 (MQ=255)
                                                                                                                                                                                                                     ggATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACt                                                                                                                                                                                   >  8:466789/1‑149 (MQ=255)
                                                                                                                                                                                                                       aTTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGa                                                    <  2:278140/274‑1 (MQ=255)
                                                                                                                                                                                                                          tGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACa                                                                                                                                                                                                                                                           >  2:338192/1‑72 (MQ=255)
                                                                                                                                                                                                                          tGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGgaagaa                                                                                        >  2:604507/1‑235 (MQ=255)
                                                                                                                                                                                                                           ggCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATAc                                                                                                                                                                             <  6:463163/149‑1 (MQ=255)
                                                                                                                                                                                                                           ggCCCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTggg                                                                                                                                                                                                  >  5:854948/1‑128 (MQ=255)
                                                                                                                                                                                                                           ggCCCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTggg                                                                                                                                                                                                  <  6:854948/128‑1 (MQ=255)
                                                                                                                                                                                                                              gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGAt                                                                                                                                                                                                                                                       <  4:338192/72‑1 (MQ=255)
                                                                                                                                                                                                                              gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCaa                                                                                                                                                                          >  5:301739/1‑149 (MQ=255)
                                                                                                                                                                                                                              gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAg                                                                                    <  4:604507/235‑1 (MQ=255)
                                                                                                                                                                                                                                 caGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACTCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCTATCAACTGGCTGGATCCTCTTACa                                          >  1:208802/1‑274 (MQ=255)
                                                                                                                                                                                                                                  aGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgct                                                                                                                                                                                                                 >  8:517970/1‑106 (MQ=255)
                                                                                                                                                                                                                                  aGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGcc                                                                                                                                                                      >  5:248577/1‑149 (MQ=255)
                                                                                                                                                                                                                                    gACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCtt                                                                                                                                                                    <  8:23666/149‑1 (MQ=255)
                                                                                                                                                                                                                                     aCCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCGGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                                                                                                                                                   <  7:643937/149‑1 (MQ=255)
                                                                                                                                                                                                                                     aCCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                                                                                                                                                   <  7:442425/149‑1 (MQ=255)
                                                                                                                                                                                                                                      cccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTg                                                                                                                                                                  <  8:739593/149‑1 (MQ=255)
                                                                                                                                                                                                                                      cccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTg                                                                                                                                                                  >  6:238380/1‑149 (MQ=255)
                                                                                                                                                                                                                                       ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTAAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGa                                                                                                                                                                 >  5:661687/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aaaTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTa                                                                                                                                             >  1:164483/1‑166 (MQ=255)
                                                                                                                                                                                                                                              cACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGa                                                                                                                                         <  3:164483/166‑1 (MQ=255)
                                                                                                                                                                                                                                               aCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACa                                          >  4:359566/1‑260 (MQ=255)
                                                                                                                                                                                                                                               aCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACa                                          <  2:359566/260‑1 (MQ=255)
                                                                                                                                                                                                                                                       tGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTg                                                                                                                                                 <  7:966088/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        gCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGc                                                                                                                                                <  8:546746/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTAc                                                                                                                                            >  7:1041826/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                          >  7:229425/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                                                                                                                                          >  7:1036357/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGa                                                                                                                                         <  8:837233/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAt                                                                                                                                        >  8:594012/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAttt                                                                                                                                      >  7:996674/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ccAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGt                                                                                                                                                                                                                                                   >  8:181098/1‑38 (MQ=255)
                                                                                                                                                                                                                                                                      ttcagtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGaa                                                                                                                                                             >  6:660937/6‑122 (MQ=255)
                                                                                                                                                                                                                                                                      ttcagtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGaa                                                                                                                                                             <  5:660937/117‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                  <  8:1041826/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCa                                                                                                                                  <  8:416983/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccTCAGGCTGGAGCCTGACTGAAAAGGATCCGTATAACAACGTTATCCGGACCACCATTGAAGCGCTGGCTCCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCa                                                                                                                               >  7:217188/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaac                                                                                                                                                                                                                                              >  4:291111/1‑37 (MQ=255)
                                                                                                                                                                                                                                                                           tCAGTCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTAACTCCATACCAACGCCTTTGACGAATCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                             >  8:1036343/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           tCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                             >  5:713360/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           tCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACg                                                                                                                             >  7:245440/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                               gCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCAtt                                                                                                                         >  7:537042/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                cTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCa                                                                                                                  <  2:657229/155‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cTGGAGCCTGACGGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACTCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACAtt                                        <  2:144828/229‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                        
TACCACATGGGTGAAGCGGGTGCCAACTGCGTGCAGCAGGTAGCATTTACGCTCGCTGATGGGATTGAGTACATCAAAGCAGCAATCTCTGCCGGACTGAAAATTGATGACTTCGCTCCTCGCCTGTCGTTCTTCTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGCATTGCCCGCAACACCCAGATCATCATCCAGGAAGAATCAGAACTCTGCCGCACCGTCGATCCACTGGCCGGATCCTATTACATTGAGTCGCTGACCGATCAAATCGTCAAACAAGCCAGAGC  >  NC_000913/3061543‑3062081

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