Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,723,504	C→A	10.0%	L20L (CTG→CTT)	htpG ←	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,723,504	0	C	A	10.0%	214.2 / 11.5	70	L20L (CTG→CTT)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (2/5);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGC  >  NC_002947/4723365‑4723651                                                                                                                                              |                                                                                                                                                       tCAGGTCAGCATCGCCCTCGAACAGCTCTGACTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTgcag                                                                                                                                                      >  8:376645/1‑142 (MQ=255) tCAGGTCAGCATCGCCCTCGAACAGCTCTGACTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTgcag                                                                                                                                                      <  7:376645/142‑1 (MQ=255)  cAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGTA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCt                                                                                                                                              >  5:56410/1‑149 (MQ=255)  cAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGc                                                                                                                                               <  5:407957/148‑1 (MQ=255)          cATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCg                                                                                                                                      >  7:211279/1‑149 (MQ=255)          cATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTca                                                                                                                                      <  5:283573/149‑2 (MQ=255)           aTCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGGA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTcat                                                                                                                                     <  6:56410/149‑3 (MQ=255)           aTCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTcat                                                                                                                                     <  2:337536/149‑3 (MQ=255)              gCCCTCGAACCTCGCTTGCTTGGCCATGGCTTGGAAACGCAGGTTTTCGGCGGCGGCGGTATCGTTTGAAATCAGCCCCCGCACCGGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGAAGCAACTGCTGCACCTcatg                                                                                                                                    <  1:122284/149‑4 (MQ=255)                       cAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCa                                                                                                                        <  3:64715/149‑1 (MQ=255)                             tGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACCGCGAGTGGATCATGAGGTGGGGCAGTTGCTTCACCTCGGTCT‑GGACGGCCAG‑GGttt                                                                                                                 >  6:356088/1‑149 (MQ=255)                             tGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGttt                                                                                                                 >  7:224378/1‑149 (MQ=255)                                cTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCt                                                                                                              <  5:126170/149‑1 (MQ=255)                                   ggCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTc                                                                                                                                            >  5:79386/1‑118 (MQ=255)                                   ggCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTc                                                                                                                                            <  6:79386/118‑1 (MQ=255)                                   ggCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCttt                                                                                                            <  3:234288/148‑1 (MQ=255)                                    gCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCtttt                                                                                                           <  3:193053/148‑1 (MQ=255)                                          gCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTc                                                                                                    <  8:224378/149‑1 (MQ=255)                                           cTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGAAAATCAGATCCCGCA‑‑GGAAAATTTCATTGTTCAAGTACAGCAAGTGAATCATGAGTTGAAGCAGTTGATTAACGTCGGTAT‑GGAAGCCCAG‑CGTTTCCTTTTTTTTTTc                                                                                                    >  7:61294/1‑148 (MQ=255)                                               gAACCGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCcaca                                                                                                <  2:280724/148‑1 (MQ=255)                                                aaaCGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTtgt                                                                                                         <  6:413005/138‑1 (MQ=255)                                                aaaCGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTtgt                                                                                                         >  5:413005/1‑138 (MQ=255)                                                    gCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATg                                                                                          >  6:356593/1‑149 (MQ=255)                                                    gCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATg                                                                                          <  8:61294/149‑1 (MQ=255)                                                        cTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCt                                                                                      >  6:148182/1‑149 (MQ=255)                                                           gTCGGCGGCGTCGGAAGCGTTGGAAATCAACTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTc                                                                                    <  1:164712/148‑1 (MQ=255)                                                            tCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGt                                                                                                                                                  <  3:373635/87‑1 (MQ=255)                                                            tCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGt                                                                                                                                                  >  4:373635/1‑87 (MQ=255)                                                                  gcgTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTcc                                                                           <  8:211279/149‑1 (MQ=255)                                                                        gAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCt                                                                      >  2:118608/1‑148 (MQ=255)                                                                           gCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATg                                                                  <  5:148182/149‑1 (MQ=255)                                                                             gTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCtt                                                                                     >  8:115656/1‑129 (MQ=255)                                                                             gTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCtt                                                                                     <  7:115656/129‑1 (MQ=255)                                                                                  aaaTCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTgg                                                           >  6:22211/1‑149 (MQ=255)                                                                                   aaTCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGgc                                                          >  1:326595/1‑149 (MQ=255)                                                                                       aGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCt                                                      <  1:118608/149‑1 (MQ=255)                                                                                        gCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGGTTCACCTCGCTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCGTC‑GATAATCCGATCTGATCGCGTTGGCGCCTg                                                     >  1:45995/1‑149 (MQ=255)                                                                                            ccGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTggcgg                                                 >  1:227503/1‑149 (MQ=255)                                                                                                a‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCgg                                             >  1:155571/1‑149 (MQ=255)                                                                                                    gAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCTGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTc                                                                                                    <  3:107192/93‑1 (MQ=255)                                                                                                    gAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGtt                                           >  8:168750/1‑149 (MQ=255)                                                                                                    gAAGATTTAATTGTTCGAGTACAGAGAGTGAATCATGAGGTGCACCAGTTGCTTCACATCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTCTTTc                                                                                                    >  4:107192/1‑93 (MQ=255)                                                                                                        aTTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGGAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                       <  5:356593/149‑1 (MQ=255)                                                                                                        aTTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                       <  3:105271/149‑1 (MQ=255)                                                                                                         tttCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                       >  7:300426/1‑148 (MQ=255)                                                                                                          ttCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCgg                                             >  6:183919/1‑141 (MQ=255)                                                                                                          ttCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCgg                                             <  5:183919/141‑1 (MQ=255)                                                                                                            ccTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGc                                   >  6:65557/1‑149 (MQ=255)                                                                                                             cTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTggcggc                                                <  3:331434/135‑1 (MQ=255)                                                                                                             cTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTggcggc                                                >  4:331434/1‑135 (MQ=255)                                                                                                             cTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGc                                   <  2:155571/148‑1 (MQ=255)                                                                                                                gTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGtt                                <  5:65557/148‑1 (MQ=255)                                                                                                                gTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGtt                                <  5:22211/148‑1 (MQ=255)                                                                                                                 ttCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                       <  5:305176/140‑1 (MQ=255)                                                                                                                 ttCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGc                                       >  6:305176/1‑140 (MQ=255)                                                                                                                  tCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAAc                             >  8:116103/1‑149 (MQ=255)                                                                                                                     aGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTaa                              >  3:146004/1‑145 (MQ=255)                                                                                                                     aGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTaa                              <  4:146004/145‑1 (MQ=255)                                                                                                                      gTACAGCGAGTCAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GCAACCACAG‑CGTTTCCTTTTGTGTTTACACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGATAAcggcc                         >  2:11238/1‑148 (MQ=255)                                                                                                                        aCAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTcc                                                                                                               <  2:235602/62‑1 (MQ=255)                                                                                                                        aCAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTcc                                                                                                               >  1:235602/1‑62 (MQ=255)                                                                                                                          aGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCtt                                                                                     >  1:141149/1‑86 (MQ=255)                                                                                                                          aGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCtt                                                                                     <  2:141149/86‑1 (MQ=255)                                                                                                                           gCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGAttt                    <  5:152667/149‑3 (MQ=255)                                                                                                                                   ttCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATg                >  4:341126/2‑145 (MQ=255)                                                                                                                                    tCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGAt           <  7:168750/149‑1 (MQ=255)                                                                                                                                      aTGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATg          <  3:341126/148‑1 (MQ=255)                                                                                                                                       tGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATggg        >  6:48977/1‑149 (MQ=255)                                                                                                                                         aGGTGCAGCAGTTGCTTCACCTCGGTCTGGGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑AAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGAGGATGTCATGCAGATgggg       >  1:31126/1‑149 (MQ=255)                                                                                                                                          ggTGAAGCAACTGCTGCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTCAGAAACTCCGATCTGATGGCGTTGGCGCCTggcggc                                                <  1:60435/107‑1 (MQ=255)                                                                                                                                           gTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAGACGTTTCCTTTTGTGCTTCCACAGTCATGGTCTAC‑GAAACTCCGATATGATGGCGTTGGCGCATGGCGGCCGGTTCGGCCAGCGtt                                >  1:67187/1‑122 (MQ=255)                                                                                                                                           gTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGCTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGtt                                <  2:67187/121‑1 (MQ=255)                                                                                                                                             gcagcaGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGc  >  8:233132/1‑149 (MQ=255)                                                                                                                                              cagcagTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTc                                          <  6:71851/108‑1 (MQ=255)                                                                                                                                              cagcagTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTc                                          >  5:71851/1‑108 (MQ=255)                                                                                                                                              |                                                                                                                                                       TCAGGTCAGCATCGCCCTCGAACAGCTCTGGCTTGGCCAGGGCTTCGAAACGCAGCTTGTCGGCGGCGTCGGAAGCGTTGGAAATCAGCTCCCGCA‑‑GGAAGATTTCCTTGTTCGAGTACAGCGAGTGAATCATGAGGTGCAGCAGTTGCTTCACCTCGGTCT‑GGAAGCCCAG‑CGTTTCCTTTTGTGTTTCCACAGTCATGGTCTTC‑GAAACTCCGATCTGATGGCGTTGGCGCCTGGCGGCCGGTTCGGCCAGCGTTAACGGCGGATGTCATGCAGATGGGGGCTGC  >  NC_002947/4723365‑4723651

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