Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012971 550,663 G→A pseudogene (114/364 nt) ECD_RS02680 → prophage tail fiber N‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_012971550,6630GA100.0% 218.1 / NA 66pseudogene (114/364 nt)ECD_RS02680prophage tail fiber N‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (46/20);  total (46/20)

TACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATG‑‑‑GCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCGGCACGCAACGCCGAAGCCGCCTCTCAGAGCGCAGCGGCGGCAAAGAAATCCGAAACTGCAGC  >  NC_012971/550507‑550950
                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
atgacagcccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTAc                                                                                                                                                     <  1:768578/271‑1 (MQ=25)
       cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGc                                                                                                                                                   <  1:609414/273‑1 (MQ=25)
       cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCACTACCGTCTATGAAGGTTCCAGACCTGGTACGc                                                                                                                                                   <  1:203931/273‑1 (MQ=25)
       cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCAAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGc                                                                                                                                                   <  1:77083/273‑1 (MQ=25)
                                       ggatatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAGCTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGAAGGaa                                                                                                                      >  2:71287/23‑292 (MQ=37)
                                       ggatatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCNAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGCCGGAAGATGATGt                                                                                                              >  2:254064/23‑300 (MQ=37)
                                       ggatatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt                                                                                                              >  1:193918/23‑300 (MQ=37)
                                       gaatatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCNGATGAAGCAGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGt                                                                                                              >  2:586926/23‑300 (MQ=255)
                                       gaatatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACTATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt                                                                                                              <  2:732903/278‑1 (MQ=37)
                                         atatggcagcagtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACTATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGt                                                                                                              <  2:755783/278‑1 (MQ=37)
                                                   aggaaaaatctCAGGTGTGTTTAAAGATGGTGCTGGAAAAACAATACGGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAACACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGAGGAAGCCGGGCGTTACAGCATGGATGTCGATTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCGTCACATGCCGGAACCATTACCTTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGc                                                                                                  <  2:112426/290‑1 (MQ=255)
                                                    gtaaaaatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCAACCTGCTGGTTGAAGGTTTTCCGCCTTAACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCAGGAGGCa                                                                                                 >  2:446561/10‑300 (MQ=37)
                                                         aatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATGCAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGc                                                                                             >  1:451228/5‑299 (MQ=255)
                                                         aatctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGc                                                                                             >  1:451220/5‑299 (MQ=255)
                                                           tctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGACAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGcgtc                                                                                          >  1:432627/3‑300 (MQ=255)
                                                           tctCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGcgtc                                                                                          >  1:393875/3‑300 (MQ=255)
                                                                   gtgCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGCAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGa                                                                                   >  2:46860/1‑299 (MQ=255)
                                                                     gCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGCCGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGc                                                                                >  1:278721/1‑300 (MQ=255)
                                                                     gCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGc                                                                                <  1:593905/300‑1 (MQ=255)
                                                                     gCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGc                                                                                >  1:279421/1‑300 (MQ=255)
                                                                     gCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCAGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGc                                                                                <  1:79912/300‑1 (MQ=255)
                                                                     gCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCAGGATGAAGCCGGGCGTTAAAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGc                                                                                <  1:79721/300‑1 (MQ=255)
                                                                      cTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCa                                                                               >  2:134601/1‑300 (MQ=255)
                                                                       tGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACGGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGNAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAgg                                                                                 <  2:620137/297‑1 (MQ=255)
                                                                                gTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCNGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGgaa                                                                      >  2:614452/1‑299 (MQ=255)
                                                                                gTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCNGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGACGGAGGTATTGCGTCGTTCTGAGGCAACGGTGgaa                                                                      >  2:614445/1‑299 (MQ=37)
                                                                                    gggAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGc                                                                 >  1:771610/1‑300 (MQ=255)
                                                                                         aaCCAATACAGAACTGCACTATTCAACTTAAGGCAAAGCGTAACAGCACCACGGTACTGTTGAACACGGTAGCCTCTGAAAATCCCGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGGTGTTTGAATGTTTTCCGCCTCCACATGCCGGAACCATTATCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCAc                                                            <  2:539894/300‑1 (MQ=255)
                                                                                           ccAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGACTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGc                                                          >  1:255044/1‑300 (MQ=255)
                                                                                              aTACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTNGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAAc                                                       >  1:563246/1‑300 (MQ=255)
                                                                                                   gAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGctga                                                  >  2:200576/1‑297 (MQ=255)
                                                                                                   gAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGctga                                                  >  1:616924/1‑297 (MQ=255)
                                                                                                   gAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGctga                                                  >  1:17607/1‑297 (MQ=255)
                                                                                                             aTTCAACTGAAGGCAAATCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCctct                                        >  1:339316/1‑300 (MQ=255)
                                                                                                             aTTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAGATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATNGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTNTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCctat                                        >  2:607326/1‑298 (MQ=255)
                                                                                                             aTTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTCCGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCatct                                        >  2:54468/1‑300 (MQ=38)
                                                                                                             aTTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCctct                                        >  2:305923/1‑300 (MQ=255)
                                                                                                                 aaCTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCag                                     <  1:633018/299‑1 (MQ=255)
                                                                                                                 aaCTGAAGGCAAAGCGNAACAGCACCACGGTACTGNTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGNTGAAGCCGCCTCTCag                                     <  1:669284/299‑1 (MQ=38)
                                                                                                                    tGAAGGCAACGCGTAACAGCACCACGGTACTGGTGAACACGGTATCCTCGGAAAATCCGGATGAAGCCGGGCGTTACAGAATGAATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGATTTTCCTCCTTCACATGCCGGAACCATTACCGTCTANGAAGGTTCCAGACCTGGTACGCTGAATGATTTNCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAgc                                  <  2:606705/299‑1 (MQ=255)
                                                                                                                     gAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGTGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTGGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAgcgc                                <  2:493818/300‑1 (MQ=255)
                                                                                                                     gAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAGGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAgcgc                                <  2:493821/300‑1 (MQ=255)
                                                                                                                     gAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAgcgc                                <  2:494166/300‑1 (MQ=255)
                                                                                                                        ggCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGTAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAgc                             <  1:779696/300‑1 (MQ=255)
                                                                                                                        ggCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAgc                             <  1:779707/300‑1 (MQ=255)
                                                                                                                            aaGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAgcggca                         >  1:684145/1‑299 (MQ=255)
                                                                                                                              gCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCAGCCTATCAGAGCGCAgc                             >  2:317022/1‑294 (MQ=38)
                                                                                                                              gCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCAGCCTATCAGAGCGCAgcggcagc                       >  2:718311/1‑297 (MQ=38)
                                                                                                                                    cAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGaa                 >  1:267999/1‑300 (MQ=255)
                                                                                                                                    cAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGaa                 >  2:561387/1‑300 (MQ=255)
                                                                                                                                    cAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCAACCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTATCAGAGCGCAGCGGCAGCAAAGaa                 >  2:669428/1‑300 (MQ=38)
                                                                                                                                      gCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAAt               >  1:450001/1‑300 (MQ=255)
                                                                                                                                      gCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCAGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACNCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGTCATGCAGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTATCAGAGCGCAGTGGCAGCAAAGAAAt               >  2:250917/1‑300 (MQ=38)
                                                                                                                                         ccacGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATca             >  1:129310/1‑298 (MQ=255)
                                                                                                                                         ccacGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTAAGAGCGCAGCGGCAGCAAAGAAATca             >  2:160066/1‑298 (MQ=38)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCTGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAAca       >  1:720753/1‑299 (MQ=255)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAAca       >  2:477730/1‑299 (MQ=38)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCAGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAAca       >  2:196336/1‑299 (MQ=38)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAGATCAGAAAca       >  2:703213/1‑299 (MQ=38)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCAGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAAAGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAAca       >  2:159308/1‑299 (MQ=38)
                                                                                                                                              gTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGAAGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGCGCGCAGCGGCAGCAGAGAAATCAGAAAca       >  2:193370/1‑299 (MQ=38)
                                                                                                                                                 cTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCAGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCCGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGAAGGAAGATGATGTAATGAAGGAGGCATTGCGTCGTTTTGAGGAAATGGTGGCAGAAGCAGCACGCAACGCTGCAGCCGCATCTAAGAGCGCAgcagcagcagaaaaatcaaaaaaaaca    >  2:185059/1‑275 (MQ=255)
                                                                                                                                                   ggttaaCACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAACagcagc  >  1:556910/5‑300 (MQ=38)
                                                                                                                                                   ggTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAACagcagc  >  1:653683/1‑300 (MQ=38)
                                                                                                                                                   ggTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCAGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCAGAGCGCAGCGGCAGCAAAGAAATCAGAAAcagc     >  2:478587/1‑294 (MQ=38)
                                                                                                                                                   ggTGAACACGGTAGCCNCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGCTGAATGATTTTCTAGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCAGCACGCAACGCTGAAGCCGCCTCTCNGAGCGCAgttgcagtaaangaaatcagaaacgggcc  >  1:659572/1‑272 (MQ=37)
                                                                                                                                                               |                                                                                                                                                                                                                                                                                               
TACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATG‑‑‑GCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAAGATGATGTCATGCCGGAGGCATTGCGTCGTTTTGAGGCAATGGTGGAAGAAGCGGCACGCAACGCCGAAGCCGCCTCTCAGAGCGCAGCGGCGGCAAAGAAATCCGAAACTGCAGC  >  NC_012971/550507‑550950

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