Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913904,9610TG25.3% 186.2 / 14.7 75C17G (TGT→GGT) amiDN‑acetylmuramoyl‑L‑alanine amidase D
Reads supporting (aligned to +/- strand):  ref base T (12/44);  new base G (19/0);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 4.92e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAACTGGCCTGG  >  NC_000913/904711‑905206
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
gcgATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATTAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGgt                                                                                                                                                                                                                                                           <  1:341947/251‑1 (MQ=255)
     atCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGGTGAGAAGATTCTTTTGGGTGGTCGCTGCCGCCCTGTTATTTGCAGGGGGTGcc                                                                                                                                                                                                                                                      >  1:553518/1‑250 (MQ=255)
           ggCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGGTCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGGAGGATGATGATGGTTAGCGTCTGCGGTACGGCGGTGAAAACGCGGCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGaa                                                                                                                                                                                                                                               <  1:75316/251‑1 (MQ=255)
                  gTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCa                                                                                                                                                                                                                                        <  2:482190/251‑1 (MQ=255)
                    tCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGGAGGGGGTGCAGGCGAAAAAGGGAtt                                                                                                                                                                                                                                     >  1:145287/1‑251 (MQ=255)
                     ccGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGGCGGTGAGAAGATTTTTCTGGCTGGTCGCTGCCGCTCTGTTATTTGGCGGGGGTGCAGGCGAAAAAGGCAttt                                                                                                                                                                                                                                    >  1:137794/1‑250 (MQ=255)
                               tGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGGGAAAACACGTCGATGAGAAGGTTCTTTTTGCTGGGCGCTGGCCCTCGGTTTTTGGGAGGGGGTGGAGACGATAAAGGGAtttttgagcaaa                                                                                                                                                                                                                           >  2:353912/1‑240 (MQ=255)
                                  aaaaGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTTGCTGGTTGCTGCCGCCCTGTTATTTGCCGGGGGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAggg                                                                                                                                                                                                                       >  1:280803/1‑251 (MQ=255)
                                     aGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGGCGATGAGAAGAATCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGGGGGGAGGCGAAAAAGGGATTGTCGGGAAAGAGGGat                                                                                                                                                                                                                     >  1:164784/1‑250 (MQ=255)
                                           gCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGc                                                                                                                                                                                                               <  1:462124/250‑1 (MQ=255)
                                                aaGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGAATCTTCTTGGTGGTCGCGGCCGCCCTGTTATTGGGAGGGGGTGGAGGCGAAAAAGGCAATTGTGAGAAAAAGGGATATCAGCtttat                                                                                                                                                                                                          >  1:179841/1‑247 (MQ=255)
                                                  gCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAAACACGCGATTCGAAGCTTCTTTTGGCAGGTTGCCGCCGCTCTTTTATTGGAAGGGGGTGCAGGGGAAAAAAGGATTGAGCAGAAAGAGGGaaaacaaattgataac                                                                                                                                                                                                       >  1:465769/1‑236 (MQ=255)
                                                      gggAGATCGCCTGTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGGCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTTGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGAc                                                                                                                                                                                                  <  2:279688/251‑1 (MQ=255)
                                                         agaTCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGcc                                                                                                                                                                                               <  1:340682/251‑1 (MQ=255)
                                                          gaTCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCCGCGGTGAGAAGATTCTTCTGGCTGGGCGGTGCCGCTCTGTTTTTGGCAGGGGGTGGAGGGGGAAAAGGCAATGTTGGGGAAGAggggtgacagcttggaaacccaccccc                                                                                                                                                                                               >  2:469101/1‑227 (MQ=255)
                                                           aTCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGGCGGTGAGAAGATTCTTCTGGCTGGGCGCTGCCGCCCTGTTATTGGCCGTGGGTGCAGGCGAAAAAAGGATTTTCGAGAAAAAGGGATATCAGCTTGATAcccccccccc                                                                                                                                                                                               >  2:211473/1‑243 (MQ=255)
                                                              gCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGGAAACCCGGCGATGAGAAGATTCTTCTGGCTGGTTGCTGCCGCGCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGGAGGAGATCAGCATGGTAccctccggcaagcg                                                                                                                                                                                           >  1:546425/1‑240 (MQ=255)
                                                              gCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGGGGAAACGCGTCGGTGAGAAGATTCTTCTGGCTGGGCGCTGCCGCCCTTTTATTTGCCGGGGGTGGGGGCGGAAAAGGCATTGTTGCGAAAGAGGGAGAACAACGTGATACCCGACGCcaggcg                                                                                                                                                                                          >  2:335040/1‑251 (MQ=255)
                                                                cTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAAAGCGGCGGTGAGGAGATTCTTTTTGCTGGTCGCTGCCGCCCTGTTATTTGCCGGGGGGGTAGGGGGAAAAGGCATTTTCGTGAAAAAAGGATATCCGCTTGATACCCgcagcccagcggc                                                                                                                                                                                        >  2:492527/1‑239 (MQ=255)
                                                                    gAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCcaggcgcaggcg                                                                                                                                                                                    <  1:173341/251‑1 (MQ=255)
                                                                     aGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGGTGAGAAGATTCTTTTGGGTGGTTGCCGGCGCCCTGTTTTTGGCGGGGGGGGGCGGGGAAAAAGGCATTTTTGGGAAAAAGGGGTATAAGCCTGATAACCTACGCcaggggacggggg                                                                                                                                                                                   >  1:362588/1‑242 (MQ=255)
                                                                            aGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATcc                                                                                                                                                                            <  2:347805/251‑1 (MQ=255)
                                                                            aGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATcc                                                                                                                                                                            <  1:267726/251‑1 (MQ=255)
                                                                                       cgcgGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAg                                                                                                                                                                <  1:422462/251‑1 (MQ=255)
                                                                                          gggCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGACGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGGGTCGATGAGAAGATTTTTCTGGCTGGTCGCTGCCCCCCTGTTTTTGGCAGGGGGTGCAGGCGAAAAAGGCAGTGTCGGGGAAGAGGGATATCAGCTTGAAACCCCAACGCCAGCCCAGGCGGGGGATACGCGCATTTAAgagc                                                                                                                                                             >  1:434820/1‑248 (MQ=255)
                                                                                          gggCGCTGGGGGCCGATGCCGTCGGCGGTATTTATATCCACTACGAAACGGGCGGGCAAAACGGCAGTATTCTGATGGTTAGCGTCAGCCGTACGGCGGTGAAAACGCGTCGATGAGAGGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGc                                                                                                                                                             <  2:354211/251‑1 (MQ=255)
                                                                                                 gggggCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTg                                                                                                                                                        <  1:166639/249‑1 (MQ=255)
                                                                                                     gCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCAc                                                                                                                                                  <  2:220629/251‑1 (MQ=255)
                                                                                                        gATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAAACCGTCGATGAAAAGATTCTTCTGGCTGGGCGCTGCCGCTCTGTTTTTGGGAGGGGGTGCAGGGGGAAAAAGGATTGTCGGGGAAAAGGGGTATCAGATTGATACCCCACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAATGGTGGGGATCCACTac                                                                                                                                               >  1:525314/1‑251 (MQ=255)
                                                                                                                cgacgGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAg                                                                                                                                        <  2:362588/247‑1 (MQ=255)
                                                                                                                 gtcgGTATTGATATCGACTACGAAACGGTCGAGCAAAACGGAAGTATGCTGATGGTTGGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCGGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATGAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATACACTACGCCGCAgtt                                                                                                                                      <  2:211117/251‑3 (MQ=255)
                                                                                                                  tcgGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAga                                                                                                                                       <  1:362921/249‑1 (MQ=255)
                                                                                                                   cgGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGAACGGCGGAGAAAACGCCTCGATGGGAAGATTCTTCTGGATGGTCGCTGGCGCCCTGTTATTGGCAGGGTGTGAGGGCGAAAAAAGCAAAGTCGGGAAAGGGGGGGATCAGCTTTTTTCCAGACGCCAGGCGCAGGCGGGGTAGCCGCGACATTAAAGGCTgggtgtacactacacccgagatgg                                                                                                                                    >  1:389543/1‑229 (MQ=255)
                                                                                                                          gATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGAt                                                                                                                             <  1:27640/251‑1 (MQ=255)
                                                                                                                          gATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGAt                                                                                                                             <  2:300249/251‑1 (MQ=255)
                                                                                                                               cGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCt                                                                                                                         <  1:457922/250‑1 (MQ=255)
                                                                                                                                    aCGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCt                                                                                                                     <  1:198425/249‑1 (MQ=255)
                                                                                                                                     cGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGc                                                                                                                  <  1:514319/251‑1 (MQ=255)
                                                                                                                                        aaCGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGac                                                                                                               <  2:145287/251‑1 (MQ=255)
                                                                                                                                              cGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGACAAAGCGTCGATGAGAAGATTCTTCTTGCTGGGCGCTGCCGCCCAGTTTTTGGGAGGGGGGGCCGGCGGAAAAGGGATTGTCGAGAAAGAGGGGTATCAGCTTGGTAACCCACGCCAGGGGCCGGGGGCGGATCCCCGGTTTAAAGAGGTGGGGATTCCCTACACCGCCGATGATTTTTATAGATCACTGGGGGCACTGa                                                                                                          >  2:207966/1‑250 (MQ=255)
                                                                                                                                               gggCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGAcc                                                                                                        <  2:553518/251‑1 (MQ=255)
                                                                                                                                                 gCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCCCTGTTATTGGCAGGGTGTGGAGGCGAAAAAAGCATTGTCGGGAAAGAGGGATATCAGCTTGATAACCGCCGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAAGGCTGGGGATCCACTACACCGCAGATGATTTTGATAACACGCTTGCGACACTGACCg                                                                                                       >  1:18812/1‑250 (MQ=255)
                                                                                                                                                      aCGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTTGCTGCCGCTCTGTTATTGGCAGGGGGTGCAGGCGAAAAAAGCCTTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGGATCCCCGCATTAAAATGCTTGGGATCCACCACACCGCAGATGATTTTTATAGCTCGCTGGCGACACTGACCGGTaaac                                                                                                 >  1:272940/1‑249 (MQ=255)
                                                                                                                                                          cAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGt                                                                                             <  1:527262/251‑1 (MQ=255)
                                                                                                                                                            gTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCa                                                                                           <  2:179841/251‑1 (MQ=255)
                                                                                                                                                                   gATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGTATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGc                                                                                     <  1:228897/250‑1 (MQ=255)
                                                                                                                                                                       gTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTAt                                                                                <  1:349747/251‑1 (MQ=255)
                                                                                                                                                                          aGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCg                                                                                                                                                                                                    <  1:288366/133‑1 (MQ=255)
                                                                                                                                                                          aGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCg                                                                                                                                                                                                    >  1:371963/1‑133 (MQ=255)
                                                                                                                                                                          aGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCg                                                                                                                                                                                                    <  2:371963/133‑1 (MQ=255)
                                                                                                                                                                          aGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTg                                                                             <  1:323076/251‑1 (MQ=255)
                                                                                                                                                                                 gCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCcag                                                                                                                                                                                             <  2:49960/133‑1 (MQ=255)
                                                                                                                                                                                 gCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGGGTGCAGGCGAAAAAGGCATTGTCGAGAAAAAGGGGTATCAGCTTGATACCCGACGcccg                                                                                                                                                                                             >  1:49960/1‑131 (MQ=255)
                                                                                                                                                                                    gTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCgg                                                                   <  1:420211/251‑1 (MQ=255)
                                                                                                                                                                                       cggcggTGAAAACGCGTCGAGGAGAAGATGCTTCGGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGGAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGTCGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACGACACCGCAGATGATTGTGATAGCTCGCTGGCGTCACTGACCGAGAAGTAGGTCAGCTAGCATTATCTGTGCCCTGCGGTac                                                                <  2:465769/251‑1 (MQ=255)
                                                                                                                                                                                         gcggTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCg                                                                                                                                           <  2:101656/174‑1 (MQ=255)
                                                                                                                                                                                         gcggTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGGGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGGGCTGGTGATCCACTACACCg                                                                                                                                           >  1:101656/1‑174 (MQ=255)
                                                                                                                                                                                           ggTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTaccacc                                                            <  1:469101/251‑1 (MQ=255)
                                                                                                                                                                                            gTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACcg                                                           >  1:173278/1‑251 (MQ=255)
                                                                                                                                                                                              gAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACcgcg                                                         >  2:86372/1‑251 (MQ=255)
                                                                                                                                                                                                      cgTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAAc                                                  <  1:41248/251‑1 (MQ=255)
                                                                                                                                                                                                           aTGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCATTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCATGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTa                                              <  2:323793/250‑1 (MQ=255)
                                                                                                                                                                                                              agaagaTTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAAcc                                          <  1:455847/251‑1 (MQ=255)
                                                                                                                                                                                                              agaagaTTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAAcc                                          >  1:446926/1‑251 (MQ=255)
                                                                                                                                                                                                                    ttcttcTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCAt                                    >  1:482380/1‑251 (MQ=255)
                                                                                                                                                                                                                        tctGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCt                                  <  1:247255/249‑1 (MQ=255)
                                                                                                                                                                                                                               ggTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGAt                                                                                                     <  2:25373/175‑1 (MQ=255)
                                                                                                                                                                                                                               ggTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGAt                                                                                                     >  1:25373/1‑175 (MQ=255)
                                                                                                                                                                                                                               ggTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACt                           <  1:86372/249‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAAc                >  1:512606/1‑250 (MQ=255)
                                                                                                                                                                                                                                             tGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGa            >  2:132794/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATGAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAAc          >  2:208453/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAAct         <  1:361083/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAActg        <  2:482380/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAActggcc     <  2:137794/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       cAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAActggcctgg  <  2:280803/251‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
GCGATATCGTTGGCGGACGTTCCGGTGCGTATGAAAAAGAACTGCGTAAAGCACGGGAGATCGCCTTTGAGGAATTAGGCTCCCAGGCGCGGGCGCTGGGGGCCGATGCCGTCGTCGGTATTGATATCGACTACGAAACGGTCGGGCAAAACGGCAGTATGCTGATGGTTAGCGTCAGCGGTACGGCGGTGAAAACGCGTCGATGAGAAGATTCTTCTGGCTGGTCGCTGCCGCTCTGTTATTGGCAGGGTGTGCAGGCGAAAAAGGCATTGTCGAGAAAGAGGGATATCAGCTTGATACCCGACGCCAGGCGCAGGCGGCGTATCCGCGCATTAAAGTGCTGGTGATCCACTACACCGCAGATGATTTTGATAGCTCGCTGGCGACACTGACCGATAAGCAGGTCAGCTCGCATTATCTGGTCCCTGCGGTACCACCGCGATACAACGGTAAACCGCGCATCTGGCAACTGGTGCCGGAACAAGAACTGGCCTGG  >  NC_000913/904711‑905206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: