Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,191,432	C→T	100%	G72S (GGC→AGC)	yehC ←	putative periplasmic pilin chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,191,432	0	C	T	100.0%	296.7 / ‑6.7	87	G72S (GGC→AGC)	yehC	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/56);  total (31/56)

TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCG  >  NC_000913/2191336‑2191527                                                                                                 |                                                                                                tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                              <  2:1251407/100‑1 (MQ=255) tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                              >  1:2209612/1‑100 (MQ=255) tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                              <  1:2454436/100‑1 (MQ=255)  ttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatca                                                                                             <  1:330160/100‑1 (MQ=255)    ttACTTGCTGCCCGGAATTTGCCCATATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                           <  2:457340/100‑1 (MQ=255)       cTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAAt                                                                                        <  2:1449636/100‑1 (MQ=255)       cTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAAt                                                                                        <  1:2226841/100‑1 (MQ=255)       cTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAAt                                                                                        <  1:252820/100‑1 (MQ=255)         tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                      >  1:1949663/1‑100 (MQ=255)         tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                      >  2:131103/1‑100 (MQ=255)          gctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCa                                                                                     <  1:35701/100‑1 (MQ=255)          gctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCa                                                                                     >  1:2208955/1‑100 (MQ=255)           ctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCAc                                                                                    <  2:226523/100‑1 (MQ=255)             ncccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                  >  2:664611/2‑100 (MQ=255)             gcCCGGAATTTGCCCCTATTTTTTTTTTTTTTTTCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCTATCCACGc                                                                                  >  2:799081/1‑100 (MQ=255)               ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                                >  2:723202/1‑100 (MQ=255)               ccGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCt                                                                                >  1:1793825/1‑100 (MQ=255)                cGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTg                                                                               >  1:274009/1‑100 (MQ=255)                 ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                              <  2:1938683/100‑1 (MQ=255)                  gAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgca                                                                             <  1:924138/100‑1 (MQ=255)                   aaTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcag                                                                            <  2:740675/100‑1 (MQ=255)                    aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                           <  1:1816018/100‑1 (MQ=255)                      ttGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGcagcag                                                                         <  1:1902608/100‑1 (MQ=255)                       tGCCGCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGACCCTGAAGTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                        <  2:2096670/100‑1 (MQ=255)                       ctccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                        <  2:1608941/98‑1 (MQ=255)                        gtccttATTTTTTCCACTTATGGCATTAACATGAAAGGAACCTTAATTTTTTCTGGTGGTAATGAAGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                       <  2:1255037/95‑1 (MQ=255)                         ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagcag                                                                      >  1:1608952/1‑98 (MQ=255)                         ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                      <  1:2016617/100‑1 (MQ=255)                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                    >  2:1180617/1‑100 (MQ=255)                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                    >  1:1617527/1‑100 (MQ=255)                           ccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaag                                                                    <  2:5594/100‑1 (MQ=255)                              aTTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAAc                                                                 <  2:1949649/100‑1 (MQ=255)                               tttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACg                                                                <  2:961216/100‑1 (MQ=255)                                ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                               <  1:415843/100‑1 (MQ=255)                                ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                               <  2:670408/100‑1 (MQ=255)                                   tGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTAACGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGttt                                                            <  2:922402/100‑1 (MQ=255)                                    gCCACTGGTGGTGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGATAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                           >  2:2182172/1‑100 (MQ=255)                                    gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                           >  1:836192/1‑100 (MQ=255)                                     ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                          >  1:1114927/1‑100 (MQ=255)                                      cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                         <  1:1821286/100‑1 (MQ=255)                                      cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                         >  1:2015377/1‑100 (MQ=255)                                       aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                        <  1:1453956/100‑1 (MQ=255)                                          ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                     >  2:1873557/1‑100 (MQ=255)                                           gtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                    <  1:1392729/100‑1 (MQ=255)                                            tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                   >  1:102838/1‑100 (MQ=255)                                             ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                  <  1:1292407/100‑1 (MQ=255)                                               cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCACTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                <  2:92823/100‑1 (MQ=255)                                               cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                >  2:2007803/1‑100 (MQ=255)                                                gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                               >  2:2303101/1‑100 (MQ=255)                                                  tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                             <  2:1926404/100‑1 (MQ=255)                                                   aaCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                            >  1:278375/1‑100 (MQ=255)                                                    aCATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTg                                           <  2:1325351/100‑1 (MQ=255)                                                     cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                          >  2:1682376/1‑100 (MQ=255)                                                     cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                          >  2:364777/1‑100 (MQ=255)                                                      aTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCa                                         <  1:2131307/100‑1 (MQ=255)                                                       tGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAc                                        <  2:675335/100‑1 (MQ=255)                                                        gAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAcc                                       >  1:1205227/1‑100 (MQ=255)                                                         aaaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCa                                      <  1:1873571/100‑1 (MQ=255)                                                           aGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATc                                    <  1:2551377/100‑1 (MQ=255)                                                             gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                  >  1:2510380/1‑100 (MQ=255)                                                             gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                  <  2:1028605/100‑1 (MQ=255)                                                             gAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCAc                                  <  1:1184298/100‑1 (MQ=255)                                                               aCCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActt                                >  1:2295958/1‑100 (MQ=255)                                                                ccTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                               <  1:1852764/100‑1 (MQ=255)                                                                ccTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                               >  2:1123317/1‑100 (MQ=255)                                                                 cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                              <  2:1563013/100‑1 (MQ=255)                                                                 cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                              <  1:977745/100‑1 (MQ=255)                                                                  tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                             <  1:635310/100‑1 (MQ=255)                                                                  tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                             <  2:1661825/100‑1 (MQ=255)                                                                  tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                             <  1:2053310/100‑1 (MQ=255)                                                                    aaTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTa                           <  2:1114919/100‑1 (MQ=255)                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                        <  1:75106/100‑1 (MQ=255)                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                        <  1:223493/100‑1 (MQ=255)                                                                       tttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAttt                        <  1:728561/100‑1 (MQ=255)                                                                          ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                     <  1:138717/100‑1 (MQ=255)                                                                           tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                    <  1:756758/100‑1 (MQ=255)                                                                           tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                    <  2:1328370/100‑1 (MQ=255)                                                                              ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                 >  2:1342289/1‑100 (MQ=255)                                                                              ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg                 >  1:1353943/1‑100 (MQ=255)                                                                               gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                <  1:2167158/100‑1 (MQ=255)                                                                               gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                >  2:533463/1‑100 (MQ=255)                                                                                tggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg               <  1:1716141/100‑1 (MQ=255)                                                                                 ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt              >  1:1478515/1‑100 (MQ=255)                                                                                 ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt              <  1:2546314/100‑1 (MQ=255)                                                                                  gtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa             <  1:2020038/100‑1 (MQ=255)                                                                                     aTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat          <  1:1180625/100‑1 (MQ=255)                                                                                       gACGTATCGCTATAATCAATCCACGCCTGCAGGAGCGAAGAACGGTTTCCCTGGTTCATTACCTGAACCCACACTTCTTTATTTTATGCCGGGTCaataa        >  2:2541930/1‑100 (MQ=255)                                                                                         cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg      <  1:716612/100‑1 (MQ=255)                                                                                            aTCGCTATCATCATTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcg    <  2:250444/99‑1 (MQ=255)                                                                                             tCGCTATCATCAATCCTCGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg  <  2:337582/100‑1 (MQ=255)                                                                                             tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg  <  1:2284197/100‑1 (MQ=255)                                                                                                 |                                                                                                TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCG  >  NC_000913/2191336‑2191527

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