Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,912,813	A→T	10.6%	L629Q (CTG→CAG)	ccmF ←	holocytochrome c synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,912,813	0	A	T	10.6%	177.0 / 9.5	66	L629Q (CTG→CAG)	ccmF	holocytochrome c synthetase
Reads supporting (aligned to +/- strand):  ref base A (31/27);  new base T (4/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AGGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGAC  >  NC_002947/4912673‑4912956                                                                                                                                              |                                                                                                                                                agGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCg                                                                                                                                         >  5:86039/1‑149 (MQ=255)    aaCACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCACGCAACCCGCCCAGGGCGGTCTGCTGAcc                                                                                                                                          <  1:202042/146‑1 (MQ=255)    aaCACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGGCGGTCTGCTGAcc                                                                                                                                          >  2:202042/1‑145 (MQ=255)        cTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                   >  2:21335/1‑148 (MQ=255)         tGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                   >  8:118432/1‑147 (MQ=255)         tGCAAGGGGCACTACCAAGATCCAACGCTTCATGCAGCTGCTCCAGACACGCACAGGACATCCCGCACCCGGGTCTTGACCTAGCCGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCCACCCGCCCAGGCCGGTCCGCAGACCGCccagcca                                                                                                                                 >  1:350695/1‑149 (MQ=255)            aaGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCccagc                                                                                                                                   >  2:342019/1‑144 (MQ=255)              ggggCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                             >  6:222287/1‑148 (MQ=255)                   aCTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                             >  4:128377/1‑143 (MQ=255)                   aCTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                             <  3:128377/143‑1 (MQ=255)                      aCCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACg                                                                                                                    >  1:126674/1‑149 (MQ=255)                       ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACg                                                                                                                    <  1:180074/148‑1 (MQ=255)                       ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGt                                                                                                                   <  2:222942/149‑1 (MQ=255)                       ccATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGt                                                                                                                   <  6:86039/149‑1 (MQ=255)                                    cTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTcagca                                                                                                                                              >  4:144780/1‑109 (MQ=255)                                         tGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAcgc                                                                                                 >  4:341947/1‑149 (MQ=255)                                                gctCCAGACACGCCTAGGACATCAGGCCCCCGGGGCTTGACCTTTACGCGGTATCGCCGGTCGAGGGCCGCCA‑GCAACCCCCCCAGCCGGGGCAGCTGACCGCCCAGCCAGATCCACCGGACGTAAGGCTTAATATGTACGCGCACAGc                                                                                          <  2:350695/149‑1 (MQ=255)                                                 ctCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAc                                                                                                   <  5:220030/139‑1 (MQ=255)                                                 ctCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTAc                                                                                                   >  6:220030/1‑139 (MQ=255)                                                    cAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAg                                                                                      >  3:2750/1‑149 (MQ=255)                                                       acacGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcg                                                                                   >  1:72016/1‑149 (MQ=255)                                                        cacGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcg                                                                                   <  4:2750/148‑1 (MQ=255)                                                         acGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGgcgc                                                                                  >  7:206786/1‑148 (MQ=255)                                                           gCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACCGCCCAGGCGCCAt                                                                               >  6:9716/1‑149 (MQ=255)                                                               aGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATtctc                                                                           >  3:277235/1‑149 (MQ=255)                                                                  acctCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGGCAGCAGACCGCCCAGCCCGATCCAGCGGACGTAAGGCTTAGTATGTACGCGCACAGCCCAGGCGCCATTCTCCAg                                                                        <  8:238831/146‑1 (MQ=255)                                                                     tCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTgg                                                                     >  1:41449/1‑149 (MQ=255)                                                                      cACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCTCCCAGCCAGATc                                                                                                                             <  3:3187/92‑1 (MQ=255)                                                                      cACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATc                                                                                                                             >  4:3187/1‑92 (MQ=255)                                                                        cGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTc                                                                  >  6:341017/1‑149 (MQ=255)                                                                        cGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTc                                                                  <  5:222287/149‑1 (MQ=255)                                                                          cACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                <  1:342019/149‑1 (MQ=255)                                                                           aCCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                >  2:325430/1‑148 (MQ=255)                                                                           aCCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGc                                                                >  8:114634/1‑148 (MQ=255)                                                                              cGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAg                                                            <  2:72016/149‑1 (MQ=255)                                                                                     tgaccttgacGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGATCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACg                                                     >  8:36006/1‑149 (MQ=255)                                                                                        ccttgacGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTAc                                                  <  8:206786/149‑1 (MQ=255)                                                                                                 ggTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                         <  1:21335/149‑1 (MQ=255)                                                                                                 ggTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                         <  5:47777/149‑1 (MQ=255)                                                                                                 ggTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGATCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCgg                                         <  7:36006/149‑1 (MQ=255)                                                                                                  gTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                      <  8:86759/119‑1 (MQ=255)                                                                                                  gTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTg                                                                      >  7:86759/1‑119 (MQ=255)                                                                                                   tAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGt                                       <  1:325430/149‑1 (MQ=255)                                                                                                     gcgcCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACCTGTACGTcgcgcttg                                      <  2:126674/148‑5 (MQ=255)                                                                                                       gcCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGCCGTACa                                                 >  2:252734/1‑135 (MQ=255)                                                                                                       gcCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                 <  1:247658/135‑1 (MQ=255)                                                                                                       gcCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                 <  1:252734/135‑1 (MQ=255)                                                                                                       gcCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACa                                                 >  2:247658/1‑135 (MQ=255)                                                                                                        ccGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGc                                  >  5:250674/1‑149 (MQ=255)                                                                                                         cGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGcc                                 <  3:341947/149‑1 (MQ=255)                                                                                                           gTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAgc                                                           <  4:17546/121‑1 (MQ=255)                                                                                                           gTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAgc                                                           >  3:17546/1‑121 (MQ=255)                                                                                                            tCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTcgcg                                          >  6:121579/1‑137 (MQ=255)                                                                                                            tCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTcgcg                                          <  5:121579/137‑1 (MQ=255)                                                                                                                gggCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                            <  1:349646/147‑1 (MQ=255)                                                                                                                gggCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGt                            >  2:349646/1‑147 (MQ=255)                                                                                                                 ggCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGa                         >  7:197769/1‑149 (MQ=255)                                                                                                                      ccA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGccc                    <  7:114634/149‑1 (MQ=255)                                                                                                                       cA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCg                   >  3:84749/1‑149 (MQ=255)                                                                                                                             aCCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCg              >  6:206114/1‑149 (MQ=255)                                                                                                                                 gCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACCGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtca          >  7:251948/1‑149 (MQ=255)                                                                                                                                    ctgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatca       >  8:43734/6‑149 (MQ=255)                                                                                                                                     tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc        >  2:92948/5‑147 (MQ=255)                                                                                                                                     tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc        <  1:97566/143‑1 (MQ=255)                                                                                                                                     tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc        <  1:92948/143‑1 (MQ=255)                                                                                                                                     tgggcGGTCTGCTGACCGGCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGtcatc        >  2:97566/5‑147 (MQ=255)                                                                                                                                         cGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGAc  <  4:277235/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                AGGAACACTGCAAGGGGCACTACCATGATCCAACGCTTCATGCAGCTGCTCCAGACACGCCCAGGACATCACGCACCCGGGTCTTGACCTTGACGCGGTAGCGCCGGTCGAGGGCCGCCA‑GCAACCCGCCCAGGCCGGTCAGCAGACCGCCCAGCCAGATCCAGCGGACGTAAGGCTTAATATGTACGCGCACAGCCCAGGCGCCATTCTCCAGTGGCTCGCCAAGCGCGACGTACAGGTCGCGGGTAAAGCCGGCGTCGATGCCCGCTTCGGTCATCATCGAC  >  NC_002947/4912673‑4912956

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: ‑