Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	576,137	T→G	pseudogene (689/1040 nt)	nmpC ←	DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	576,137	0	T	G	100.0%	301.3 / NA	83	pseudogene (689/1040 nt)	nmpC	DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (39/44);  total (39/44)

TCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACAC  >  NC_000913/575908‑576368                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        tCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTTTGCATTTTTACCGGAAGCGACTACTGCAgt                                                                                                                                                                                                                                     >  2:562972/1‑233 (MQ=255)      gATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAAc                                                                                                                                                                                                                               >  1:26912/1‑235 (MQ=255)      gATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAAc                                                                                                                                                                                                                               >  1:423258/1‑235 (MQ=255)          cTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTcc                                                                                                                                                                                                                          <  1:376397/236‑1 (MQ=255)            tACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCt                                                                                                                                                                                                                        <  2:651373/236‑1 (MQ=255)               cTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                                                                                                                                                                     <  2:736743/236‑1 (MQ=255)                 ttAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                                                                                                                                                                     >  1:558474/1‑234 (MQ=255)                      ttGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACtt                                                                                                                                                                                                              <  2:51087/236‑1 (MQ=255)                         cAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAgg                                                                                                                                                                                                                                     >  1:554751/1‑210 (MQ=255)                         cAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAgg                                                                                                                                                                                                                                     <  2:554751/210‑1 (MQ=255)                         cAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAAc                                                                                                                                                                                                                >  1:95424/1‑231 (MQ=255)                         cAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAg                                                                                                                                                                                                           >  2:334980/1‑236 (MQ=255)                                         ggacgCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGa                                                                                                                                                                                           <  2:253827/236‑1 (MQ=255)                                                gACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCATCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTGCCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGc                                                                                                                                                                                     <  1:461187/235‑1 (MQ=255)                                                     aaaTCGAACTGATATTGTTCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGt                                                                                                                                                                               <  1:334980/236‑1 (MQ=255)                                                                    tgtgCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGa                                                                                                                                                                 <  1:321926/235‑1 (MQ=255)                                                                      tgCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAAt                                                                                                                                                               >  2:92768/1‑235 (MQ=255)                                                                        cAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATACATTGCACCGATACCGAATACt                                                                                                                                                            >  2:337640/1‑236 (MQ=255)                                                                          aCAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCtt                                                                                                                                                           >  2:344617/1‑235 (MQ=255)                                                                          aCAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCtt                                                                                                                                                           <  1:344617/235‑1 (MQ=255)                                                                          aCAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTc                                                                                                                                                          >  1:222661/1‑236 (MQ=255)                                                                           cAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                                                                                                                                                                     <  1:108851/176‑1 (MQ=255)                                                                           cAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCa                                                                                                                                                                                                                     >  2:108851/1‑176 (MQ=255)                                                                                tCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTc                                                                                                                                                    >  2:666137/1‑236 (MQ=255)                                                                                tCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAAATTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTc                                                                                                                                                    >  1:711596/1‑236 (MQ=255)                                                                                  gAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                                                                                                                                                         <  2:615216/229‑1 (MQ=255)                                                                                   aaGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCAt                                                                                                                                                  <  2:95424/235‑1 (MQ=255)                                                                                       tttGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGt                                                                                                                                              <  2:526546/235‑1 (MQ=255)                                                                                       tttGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGt                                                                                                                                              <  1:337640/235‑1 (MQ=255)                                                                                                  ttaGCAACGAAGTGATCAGCAAATACAGTCATATTCAGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAAcc                                                                                                                                  <  2:241933/236‑1 (MQ=255)                                                                                                        aCGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCt                                                                                                                                                            <  2:347166/204‑1 (MQ=255)                                                                                                        aCGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCt                                                                                                                                                            >  1:347166/1‑204 (MQ=255)                                                                                                            aGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCa                                                                                                                        <  2:610547/236‑1 (MQ=255)                                                                                                            aGTGATCAGCAAATACAGTAATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCa                                                                                                                        <  2:632107/236‑1 (MQ=255)                                                                                                                 gcAGCAAATACGGTAATATTCTGGGTTTCAGAATAGGGAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGCCCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGtt                                                                                                                   <  2:232299/235‑1 (MQ=255)                                                                                                                    gCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTAcc                                                                                                                <  1:398044/236‑1 (MQ=255)                                                                                                                               cATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAg                                                                                                      >  2:108652/1‑235 (MQ=255)                                                                                                                                   ttCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCaga                                                                                                                                        >  1:149723/1‑197 (MQ=255)                                                                                                                                   ttCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCaga                                                                                                                                        <  2:149723/197‑1 (MQ=255)                                                                                                                                   ttCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATc                                                                                                 >  1:101858/1‑236 (MQ=255)                                                                                                                                   ttCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATc                                                                                                 >  1:694061/1‑236 (MQ=255)                                                                                                                                     cTGGGTTTCAGAATAGGGAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGa                                                                                               <  2:664144/236‑1 (MQ=255)                                                                                                                                        ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGaa                                                                                              >  1:685234/1‑234 (MQ=255)                                                                                                                                        ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAAt                                                                                            >  2:124815/1‑236 (MQ=255)                                                                                                                                        ggTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAAt                                                                                            <  1:135724/236‑1 (MQ=255)                                                                                                                                          tttCAGAATAGGTAGTGGCCAGGTATATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATc                                                                                           >  2:15611/1‑235 (MQ=255)                                                                                                                                            tCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTTAGTGTAGTTACTGAAATCGcc                                                                                        >  2:270348/1‑235 (MQ=255)                                                                                                                                                  tAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGGATTCATAGGTAGCAGAGAAACCGAAGCCATAACCGTTACCTTCAGTGTAGTTATCGCAATCGCTACGATc                                                                                  >  1:765767/1‑236 (MQ=255)                                                                                                                                                         tGGCCAGGTATATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                             <  1:15611/234‑1 (MQ=255)                                                                                                                                                         tGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                             <  1:124815/234‑1 (MQ=255)                                                                                                                                                          ggCCAGTTAAATGTTGTTAGCGTCATATTGCAGTCCTGCTGCCCAATCTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAATAACTTTCCCTTCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTg                                                                           <  2:711596/235‑1 (MQ=255)                                                                                                                                                          ggCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCACTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGc                                                                          >  2:702404/1‑236 (MQ=255)                                                                                                                                                                    aTGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTa                                                                   <  2:558474/233‑1 (MQ=255)                                                                                                                                                                     tgttgtTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTg                                                                <  1:666137/235‑1 (MQ=255)                                                                                                                                                                     tgttgtTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGAACGTTTTTGCATTGGTACTGa                                                               >  2:432317/1‑236 (MQ=255)                                                                                                                                                                         gttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcag                                                           >  1:567950/1‑236 (MQ=255)                                                                                                                                                                          ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGtt                                                                                                                   >  1:698004/1‑179 (MQ=255)                                                                                                                                                                          ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGtt                                                                                                                   <  2:698004/179‑1 (MQ=255)                                                                                                                                                                          ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                          >  1:93323/1‑236 (MQ=255)                                                                                                                                                                          ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                          >  1:120748/1‑236 (MQ=255)                                                                                                                                                                          ttAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGagcagc                                                          >  1:215022/1‑236 (MQ=255)                                                                                                                                                                           tAGCGTCATATTTCAGACCTGGCGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAgcagca                                                         <  2:603253/236‑1 (MQ=255)                                                                                                                                                                           tAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAgcagca                                                         <  2:120748/236‑1 (MQ=255)                                                                                                                                                                            aGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCaa                                                        >  1:447960/1‑236 (MQ=255)                                                                                                                                                                                  atatTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCa                                                  <  1:270348/236‑1 (MQ=255)                                                                                                                                                                                    atTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCAGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGa                                                >  1:141326/1‑236 (MQ=255)                                                                                                                                                                                    atTACAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGa                                                <  1:739242/236‑1 (MQ=255)                                                                                                                                                                                      ttCAGTCCTGCCGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGAcc                                              <  2:737349/236‑1 (MQ=255)                                                                                                                                                                                       tCAGGCCTGCGGCCCAAACTTCTTCATTTTTACCGGAAGCAAATACTGCAGGAATAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCa                                             <  2:585448/236‑1 (MQ=255)                                                                                                                                                                                        cAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTAc                                                                                                                 <  1:611072/167‑1 (MQ=255)                                                                                                                                                                                        cAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTAc                                                                                                                 >  2:611072/1‑167 (MQ=255)                                                                                                                                                                                        cAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCa                                             >  1:117062/1‑235 (MQ=255)                                                                                                                                                                                                cGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAAcc                                       >  2:57332/1‑233 (MQ=255)                                                                                                                                                                                                cGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAAcc                                       <  1:57332/233‑1 (MQ=255)                                                                                                                                                                                                          ttCTGCATTTTTACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGt                          <  2:117062/236‑1 (MQ=255)                                                                                                                                                                                                                  ttttACCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTa                  <  1:92768/236‑1 (MQ=255)                                                                                                                                                                                                                      aCCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                                                                                                                                                         <  2:681951/97‑1 (MQ=255)                                                                                                                                                                                                                      aCCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCg                                                                                                                                                         >  1:681951/1‑97 (MQ=255)                                                                                                                                                                                                                      aCCGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGAt              >  1:762158/1‑236 (MQ=255)                                                                                                                                                                                                                       ccGGAAGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATa             <  2:762158/236‑1 (MQ=255)                                                                                                                                                                                                                            aGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                             <  2:321082/167‑1 (MQ=255)                                                                                                                                                                                                                            aGCAAATACTGCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGtttt                                                                             >  1:321082/1‑167 (MQ=255)                                                                                                                                                                                                                                  tACTGCAGGAATGACTTTCCCTGCATTAACTTGAGTGTCGGTTCGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAacac  <  2:609899/236‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        TCCAAGATCCTTACCTTTAGATTGCAGGTAAGCAACGGACGGACGCAGACCGAAATCGAACTGATATTGTGCAACAGCTTCGAAGTTTTGGGCTTTATTAGCAACGAAGTGATCAGCAAATACAGTCATATTCTGGGTTTCAGAATAGGTAGTGGCCAGGTAAATGTTGTTAGCGTCATATTTCAGACCTGCGGCCCAAACTTCTGCATTTTTACCGGAAGCAAATACTTCAGGAAGAACTTTCCCTGCATTAACTTGAGTGTCGGTACGATCAGATTTCGCATAAGTTGCACCGATACCGAATCCTTCGTATTCATAGGTAGCAGAGAAACCGAAGCCATCACCGTTACCTTCAGTGTAGTTATCGAAATCGCTACGATCGTTTTTGCCTTGGTACTGAGCAGCAAAGTTCAGACCATCAACCAGACCAAAGAAGTCGTTGTTACGATAGGTTGCAACAC  >  NC_000913/575908‑576368

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