Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,992,624 T→G 48.6% L491R (CTG→CGG)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,6240TG48.6% 14.1 / 303.3 165L491R (CTG→CGG) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (39/46);  new base G (38/42);  total (77/88)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 177 total aligned reads displayed.

AATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTG  >  NC_000913/3992362‑3992886
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                       
  tacccTTAGCGCTCCCGCCAGTCCGCCGGATATCTGCGTCCTTCGGCCTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAACGTTACGCTGGTAAACCCGCAGCTTTCACCCGCTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCTCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCa                                                                                                                                                                                                                                                            <  4:181450/271‑1 (MQ=255)
               tcCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGTCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAAACGCGCGCCAACTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGAGGCGAGGGCAGGGTTTAAAGGcc                                                                                                                                                                                                                                               >  2:214333/1‑274 (MQ=255)
                             ggATATCGGCGGGCTGACGCGTAAACTGTAGGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGGTGGTAAACCCGCCGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCa                                                                                                                                                                                                                                <  2:30572/275‑1 (MQ=255)
                                       gTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGGTGGTAAACCCGCAGATTTCAACCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGt                                                                                                                                                                                                                       <  2:111976/274‑1 (MQ=255)
                                                                            tACCAGGTAAAGTGACGCTGGTAAAACCGCAGCTGTCGCCCGAGCTCTCGGAACCGAATGTGACCTTTATTTATGGGCCGCCGGGCCGGGCTAACCGTTCCGGTTGGTAGCTGTATAACCGCGCGCCAAATATGGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATAAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTg                                                                                                                                                                                <  4:339277/275‑1 (MQ=255)
                                                                                            gCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTACAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGcag                                                                                                                                                                >  2:33974/1‑273 (MQ=255)
                                                                                                              ttCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCg                                                                                                                                                                                                                                                                  <  2:378764/160‑1 (MQ=255)
                                                                                                                  cccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACt                                                                                                                                                                                                                                                                         >  5:217982/1‑149 (MQ=255)
                                                                                                                   ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtg                                                                                                                                                                                                                                                                        >  6:952584/1‑149 (MQ=255)
                                                                                                                   ccGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtg                                                                                                                                                                                                                                                                        >  8:1259926/1‑149 (MQ=255)
                                                                                                                    cGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtgg                                                                                                                                                                                                                                                                       >  7:1441935/1‑149 (MQ=255)
                                                                                                                     gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                                                                                                                                                                                      <  5:1173604/149‑1 (MQ=255)
                                                                                                                     gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                                                                                                                                                                                      <  7:687988/149‑1 (MQ=255)
                                                                                                                     gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACtggt                                                                                                                                                                                                                                                                      <  7:836134/149‑1 (MQ=255)
                                                                                                                     gATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggt                                                                                                                                                                                                                                                                      <  7:331581/149‑1 (MQ=255)
                                                                                                                      aTCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACgggtg                                                                                                                                                                                                                                                                     >  5:434192/1‑149 (MQ=255)
                                                                                                                          ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGt                                                                                                                                                                                                                                                                 >  7:155039/1‑149 (MQ=255)
                                                                                                                          ctcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGt                                                                                                                                                                                                                                                                 <  5:558145/149‑1 (MQ=255)
                                                                                                                           tcGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTg                                                                                                                                                                                                                                                                >  5:346277/1‑149 (MQ=255)
                                                                                                                            cGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTgg                                                                                                                                                                                                                                                               <  7:208350/149‑1 (MQ=255)
                                                                                                                              gAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACt                                                                                                                                                                                                                                                                         >  6:1052327/1‑137 (MQ=255)
                                                                                                                              gAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACt                                                                                                                                                                                                                                                                         <  5:1052327/137‑1 (MQ=255)
                                                                                                                                aCCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCAt                                                                                                                                                                                                                                                           <  6:1460764/149‑1 (MQ=255)
                                                                                                                                 ccGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATg                                                                                                                                                                                                                                                          <  8:1373383/149‑1 (MQ=255)
                                                                                                                                 ccGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                                                                                                                                                          >  7:460228/1‑149 (MQ=255)
                                                                                                                                 ccGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATg                                                                                                                                                                                                                                                          <  5:160612/149‑1 (MQ=255)
                                                                                                                                  cGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATgg                                                                                                                                                                                                                                                         >  5:1312089/1‑149 (MQ=255)
                                                                                                                                   gAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGt                                                                                                                                                                                                                                                        >  8:260508/1‑149 (MQ=255)
                                                                                                                                    aaTCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGAGCTCGCGCCCCCGTTTGTATATTAAAGGTAAAGGCATTGTCGATTTGCCTAAAGTTG‑AGGGTCGGGTCGCCGACTTGTTGCCCCATCTCCCGTCGCGCATACCGGccacgac                                                                                                                        >  4:81382/1‑275 (MQ=255)
                                                                                                                                       cTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTaa                                                                                                                                                                                                                                                    <  5:520232/149‑1 (MQ=255)
                                                                                                                                       cTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTaa                                                                                                                                                                                                                                                    <  7:785143/149‑1 (MQ=255)
                                                                                                                                       cTGACCTTTATTTATGTGCCGCCGGGCCGAGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTaa                                                                                                                                                                                                                                                    <  5:891267/149‑1 (MQ=255)
                                                                                                                                        tGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCACGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAAc                                                                                                                                                                                                                                                   <  8:1002637/149‑1 (MQ=255)
                                                                                                                                         gACTTTTATTTATGTGTCGCCGGGCCGTGCGAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgacaccga                                                                                                                   <  4:214333/275‑1 (MQ=255)
                                                                                                                                              ttatttatGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgc                                                                                                               <  2:189217/274‑1 (MQ=255)
                                                                                                                                                       tgCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGACGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCGCCGACCCCGACGGCACTCTACAg                                                                                                      >  2:6039/1‑274 (MQ=255)
                                                                                                                                                         cgcccggaccgggCTAACGTTTAAGGTNNTTATCTATATAACCAGTCGCCAAATATTGGNTCGGTCATCAGCCATCAGCCTCTGGAATATAACCGTTACCTGAGTAGACTGGTGGCGTGGGCATGGTGTAACGGNCTGCTGACCTCGCGCACCCGTTTGTGTATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGc                                                                                                     <  3:44588/270‑1 (MQ=255)
                                                                                                                                                                gccgggcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGt                                                                                                                                                                                                                           <  5:387892/149‑1 (MQ=255)
                                                                                                                                                                     gcTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGta                                                                                                                                                                                                                      <  7:1261526/149‑1 (MQ=255)
                                                                                                                                                                      cTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGtat                                                                                                                                                                                                                     <  5:784405/149‑1 (MQ=255)
                                                                                                                                                                       tAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGtata                                                                                                                                                                                                                    <  7:195146/149‑1 (MQ=255)
                                                                                                                                                                        aaCCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTatat                                                                                                                                                                                                                   >  5:1482045/1‑149 (MQ=255)
                                                                                                                                                                            gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                                                                                                                                               <  8:1441935/149‑1 (MQ=255)
                                                                                                                                                                            gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                                                                                                                                               <  8:1115693/149‑1 (MQ=255)
                                                                                                                                                                            gTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTaaa                                                                                                                                                                                                               >  6:148364/1‑149 (MQ=255)
                                                                                                                                                                             ttCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAg                                                                                                                                                                                                              >  5:1472326/1‑149 (MQ=255)
                                                                                                                                                                                  gTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAAc                                                                                                                                                                                                         >  7:1195816/1‑149 (MQ=255)
                                                                                                                                                                                     ggTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGc                                                                                                                                                                                                      <  6:1417261/149‑1 (MQ=255)
                                                                                                                                                                                         tCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTg                                                                                                                                                                                                  >  5:850151/1‑149 (MQ=255)
                                                                                                                                                                                             tataACCGCGCGCCAAATATTGAGTCGATGATCAGCCAGCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGt                                                                                                                                                                                                 <  8:1189232/146‑1 (MQ=255)
                                                                                                                                                                                             tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGa                                                                                                                                                                                              <  8:155039/149‑1 (MQ=255)
                                                                                                                                                                                             tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGa                                                                                                                                                                                              <  6:1109917/149‑1 (MQ=255)
                                                                                                                                                                                             tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTgg                                                                                                                                                                                                                                                               >  5:513861/1‑84 (MQ=255)
                                                                                                                                                                                             tataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTgg                                                                                                                                                                                                                                                               <  6:513861/84‑1 (MQ=255)
                                                                                                                                                                                              ataACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCAtt                                                                                                                                                                                                   >  8:1279825/1‑143 (MQ=255)
                                                                                                                                                                                                aaCCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGAttt                                                                                                                                                                                           >  7:149545/1‑149 (MQ=255)
                                                                                                                                                                                                   cgcgcgCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGcc                                                                                                                                                                                        >  6:797720/1‑149 (MQ=255)
                                                                                                                                                                                                      gcgcCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAg                                                                                                                                                                                                              <  5:1152186/124‑1 (MQ=255)
                                                                                                                                                                                                        gcCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGt                                                                                                                                                                                  <  7:260508/149‑1 (MQ=255)
                                                                                                                                                                                                          cAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑aa                                                                                                                                                                                    >  6:917314/1‑145 (MQ=255)
                                                                                                                                                                                                          cAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑aa                                                                                                                                                                                    <  5:917314/145‑1 (MQ=255)
                                                                                                                                                                                                           aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTa                                                                                                                                                                                                                                                     <  7:1123604/80‑1 (MQ=255)
                                                                                                                                                                                                           aaaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTa                                                                                                                                                                                                                                                     >  8:1123604/1‑80 (MQ=255)
                                                                                                                                                                                                            aaTATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCa                                                                                                                                                                              <  7:914717/149‑1 (MQ=255)
                                                                                                                                                                                                             atatTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAg                                                                                                                                                                             >  8:598568/1‑149 (MQ=255)
                                                                                                                                                                                                              tatTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAgg                                                                                                                                                                            >  6:684048/1‑149 (MQ=255)
                                                                                                                                                                                                              tatTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTaa                                                                                                                                                                                                                                                    <  6:349803/78‑1 (MQ=255)
                                                                                                                                                                                                                ttGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGgag                                                                                                                                                                          >  5:1069103/1‑149 (MQ=255)
                                                                                                                                                                                                                  gAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGgcgc                                                                                                             >  2:222148/1‑208 (MQ=255)
                                                                                                                                                                                                                   aGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGgag                                                                                                                                                                          >  8:905265/1‑146 (MQ=255)
                                                                                                                                                                                                                   aGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATg                                                                                                                                                                       <  5:59384/149‑1 (MQ=255)
                                                                                                                                                                                                                    gTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATgg                                                                                                                                                                      >  6:1208454/1‑149 (MQ=255)
                                                                                                                                                                                                                       gATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCg                                                                                                                                                                   >  8:292335/1‑149 (MQ=255)
                                                                                                                                                                                                                          catcagccatcagcCGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTTGCCg                                                                                                                                                                >  8:879688/1‑149 (MQ=255)
                                                                                                                                                                                                                             cagccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACg                                                                                                                                                             >  8:195523/1‑149 (MQ=255)
                                                                                                                                                                                                                              agccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACgt                                                                                                                                                            >  7:1093922/1‑149 (MQ=255)
                                                                                                                                                                                                                                ccatcagccGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACgtgt                                                                                                                                                          <  5:952584/149‑1 (MQ=255)
                                                                                                                                                                                                                                ccatcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACgtgt                                                                                                                                                          <  5:932227/149‑1 (MQ=255)
                                                                                                                                                                                                                                   tcagccGCTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAggtg                          >  2:67968/1‑274 (MQ=255)
                                                                                                                                                                                                                                    aggcggttGGAAAAAAACGGTTACTGGAATAAACTGGTGGGGTGGGAAGGGTTTAACGGCTTGATGACCTCGCGCACCCGGTTGTATATTAAAGGGAACGGCATTGTCGATTTGCCT‑AAGTTGCAAGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAggtggt                        <  4:6039/268‑1 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCAttt                                                                                                                                                <  7:613870/149‑1 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATATAACCGTTACCTGAATAAACTGGTCGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCTCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCCACATCTGGCGATTACCGTTAACCTGGAATATGACCAGACAGTCGCCTTCCGCAACCAGCTGGTGCAtt                   >  4:126154/1‑274 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCAttt                                                                                                                                                >  8:787960/1‑149 (MQ=255)
                                                                                                                                                                                                                                          cTGGAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACAACGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATTGTTAACCTGGCATATGACCCTACAGAGGGGTTCCGCACTCAGCTGGTGCAt                    >  1:61119/1‑273 (MQ=255)
                                                                                                                                                                                                                                           tGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTc                                                                                                                                               >  6:120882/1‑149 (MQ=255)
                                                                                                                                                                                                                                            ggAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTcc                                                                                                                                              >  6:478398/1‑149 (MQ=255)
                                                                                                                                                                                                                                             gAATATAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTccc                                                                                                                                             <  8:533404/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gtataATCGGTAGCTGAATAGCCTGTTGGGGTGGTCGTGGGTTAGCGGCATGCTGACGTCGCGCACCCGGATGTATCTTACAGGTAACGTCATTGTCGATTTGCCT‑ATGTTGCAGGAGATGGTCGCCGACTTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACAACGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAt              <  4:4942/273‑1 (MQ=255)
                                                                                                                                                                                                                                               atatAACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGc                                                                                                                                           <  6:131134/149‑1 (MQ=255)
                                                                                                                                                                                                                                                tataACCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCt                                                                                                                                          <  8:460228/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaCCGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTgcg                                                                                                                                       >  6:1472557/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ccGTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTa                                                                         >  2:139765/1‑209 (MQ=255)
                                                                                                                                                                                                                                                       gTTACCTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTa                                                                                                                                   <  7:1304953/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          aCCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGACGTCGCTCTACAGCCCGTGTGAGGTCCCCCATCTGGCGATTATCGTTAACCTGGACTCTGACCCGACACCGGCGTTCCGCAATCCGGTTGTCCGTTTCAATTTCCTTAAGcgg  >  4:289287/1‑273 (MQ=255)
                                                                                                                                                                                                                                                            cTGAATAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGc                                                                                                                              >  5:20573/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcac                                                                                                                            >  5:843547/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcac                                                                                                                            <  7:1259926/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaTAAACGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcacc                                                                                                                           >  8:570188/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaCTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgac                                                                                                                        <  6:217982/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaCGGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGcaccgaca                                                                                                                       >  5:1125734/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                       
AATAACCTTAGCGTCTCCGCCAGTCCGCAGGATATCGGCGTGCTGACGCGTAAGCTGTATGCCGCGTTTGAAGCATTACCAGGTAAAGTGACGCTGGTAAACCCGCAGATTTCACCCGATCTCTCGGAACCGAATCTGACCTTTATTTATGTGCCGCCGGGCCGGGCTAACCGTTCAGGTTGGTATCTGTATAACCGCGCGCCAAATATTGAGTCGATCATCAGCCATCAGCCGCTGGAATATAACCGTTACCTGAATAAACTGGTGGCGTGGGCATGGTTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCT‑AAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTG  >  NC_000913/3992362‑3992886

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