Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,891,210	A→C	7.7%	R136R (CGT→CGG)	plpB ←	NLPA lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,891,210	0	A	C	7.7%	192.8 / 6.4	65	R136R (CGT→CGG)	plpB	NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base C (2/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CAGCACACGGGGCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTCGTT  >  NC_002947/5891069‑5891356                                                                                                                                              |                                                                                                                                                   cAGCACACGGGGCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGc                                                                                                                                             >  6:74530/1‑149 (MQ=255) cAGCACACGGGGCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGc                                                                                                                                             >  6:80928/1‑149 (MQ=255)     acacGGGGCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGtt                                                                                                                                         >  2:127765/1‑149 (MQ=255)          gggCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAACAGCGGCCGTGCGCtgctgctgc                                                                                                                                    <  5:185986/149‑1 (MQ=255)           ggCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCAGCTGTTGCTGcc                                                                                                                                   >  7:93268/1‑149 (MQ=255)            gcagcaTGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAACAGCGGCCGTGCGCtgctgctgctg                                                                                                                                  <  1:230159/149‑3 (MQ=255)             cagcaTGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCt                                                                                                                                 >  7:212860/1‑149 (MQ=255)                    gCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAg                                                                                                                           <  1:127765/148‑1 (MQ=255)                    gCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTAGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAgg                                                                                                                          <  8:284469/149‑1 (MQ=255)                                  gcgGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTg                                                                                                            >  2:210502/1‑149 (MQ=255)                                    ggAACTTCAGCTTCCTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTgcg                                                                                                          >  4:244512/1‑149 (MQ=255)                                           cAGCTTCTTCGGGTTCTCGGGGATATCCTTGGAGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCcttctt                                                                                                   >  8:155125/1‑149 (MQ=255)                                            aGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCttt                                                                                                  <  8:212860/149‑1 (MQ=255)                                                 cttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCCGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAg                                                                                               <  5:80928/147‑1 (MQ=255)                                                 cttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGGTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAgc                                                                                              <  3:348007/148‑1 (MQ=255)                                                 cttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAgc                                                                                              <  5:74530/148‑1 (MQ=255)                                                 cttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAg                                                                                               <  1:210502/147‑1 (MQ=255)                                                  ttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGcc                                                                                                                                               >  7:85240/1‑98 (MQ=255)                                                  ttcGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGcc                                                                                                                                               <  8:85240/98‑1 (MQ=255)                                                    cGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAgctcgc                                                                                          >  7:320201/1‑149 (MQ=255)                                                        ttCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAACAGCGGCCGTGCGCtgctgc                                                                                                                                       >  1:131968/1‑100 (MQ=38)                                                        ttCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAACAGCGGCCGTGCGCtgctgc                                                                                                                                       <  2:131968/100‑1 (MQ=38)                                                                           gggTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAg                                                                   <  7:155125/149‑1 (MQ=255)                                                                              tGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCT‑‑‑GCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCa                                                             >  3:146334/1‑149 (MQ=255)                                                                                 ccAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCttt                                                                                                  >  6:139963/1‑112 (MQ=255)                                                                                 ccAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCttt                                                                                                  <  5:139963/112‑1 (MQ=255)                                                                                   aGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCAcc                                                           >  3:14615/1‑149 (MQ=255)                                                                                     ggCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCttt                                                                                                  >  8:176210/1‑108 (MQ=255)                                                                                     ggCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCttt                                                                                                  <  7:176210/108‑1 (MQ=255)                                                                                      gCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGaa                                                        <  8:320201/149‑1 (MQ=255)                                                                                      gCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGaa                                                        <  2:75225/149‑1 (MQ=255)                                                                                      gCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGACCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGaa                                                        <  1:165336/149‑1 (MQ=255)                                                                                              cGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAggg                                                     >  4:299439/1‑144 (MQ=255)                                                                                              cGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAggg                                                     <  3:299439/144‑1 (MQ=255)                                                                                                ggTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCttcttc                                                                      >  8:302261/1‑125 (MQ=255)                                                                                                ggTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCttcttc                                                                      <  7:302261/125‑1 (MQ=255)                                                                                                    tttCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTc                                                                                                 <  5:192931/94‑1 (MQ=255)                                                                                                    tttCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTc                                                                                                 >  6:192931/1‑94 (MQ=255)                                                                                                    tttCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTa                                          >  4:86578/1‑149 (MQ=255)                                                                                                      tCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTaca                                        >  3:58966/1‑149 (MQ=255)                                                                                                      tCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTAGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACgtgt                                           <  4:14615/146‑1 (MQ=255)                                                                                                       cAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTAc                                               >  1:167433/1‑141 (MQ=255)                                                                                                       cAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTAc                                               <  2:167433/141‑1 (MQ=255)                                                                                                        aaGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACgt                                             >  3:190839/1‑142 (MQ=255)                                                                                                        aaGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACgt                                             <  4:190839/142‑1 (MQ=255)                                                                                                          ggTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTAcaccca                                    >  4:115621/1‑149 (MQ=255)                                                                                                                 aGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGt                             >  6:186904/1‑149 (MQ=255)                                                                                                                 aGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCATTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGt                             >  1:89099/1‑149 (MQ=255)                                                                                                                  ggccggccTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGtt                            >  3:63625/1‑149 (MQ=255)                                                                                                                     cccccctCTGCAGCAGCAGCAGCGCACCGCCGCTGTTGCTGCCCTCCGTAGGGATCGCAACGGTTGCGCCATCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGACGGGCTCTACGTGTACACCCACCACGGCTTcc                         <  4:63625/143‑1 (MQ=255)                                                                                                                         cTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAggt                     <  3:244512/149‑1 (MQ=255)                                                                                                                            cTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAggtgg                   <  5:85049/148‑1 (MQ=255)                                                                                                                             tGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGAc                                                                             >  2:6941/1‑89 (MQ=255)                                                                                                                             tGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGAc                                                                             <  1:6941/89‑1 (MQ=255)                                                                                                                                 gcagcagcagcGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGccc             >  1:305468/1‑149 (MQ=255)                                                                                                                                   agcaacagcGGCCGTGCGCTGCTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCtt           >  1:49410/1‑149 (MQ=255)                                                                                                                                     cagcagcGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGc         >  6:356987/1‑149 (MQ=255)                                                                                                                                         agcGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTc     >  2:279735/1‑149 (MQ=255)                                                                                                                                          gcgcACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTc     <  5:292478/148‑1 (MQ=255)                                                                                                                                          gcgcACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTc     >  6:292478/1‑148 (MQ=255)                                                                                                                                          gcgcACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTCg    >  2:351249/1‑149 (MQ=255)                                                                                                                                            gcACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGGCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGAGGGTGCCCTTGCCTTCGtt  >  1:251803/1‑149 (MQ=255)                                                                                                                                            gcACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTCGtt  <  5:356987/149‑1 (MQ=255)                                                                                                                                             cACGGCCTCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTAc                                               <  5:95233/103‑1 (MQ=255)                                                                                                                                             cACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTAc                                               >  6:95233/1‑103 (MQ=255)                                                                                                                                              |                                                                                                                                                   CAGCACACGGGGCAGCATGGCCGACTCAAGCTCGCGGAACTTCAGCTTCTTCGGGTTCTCGGCGATATCCTTGGGGGTGGCCAGGGCATTCTTCGGGTCTTTCAAGGTAATCAGGCCGGCCTTCTGCAGCAGCAGCAGCGCACGGCCGCTGTTGCTGCCCTCGTTAGGGATGGCAACGGTTGCGCCTTCTTTCAGCTCGCTCAGGGCCTTGACCTTCTTCGAGTAGCCACCGAAGGGCTCTACGTGTACACCCACCACGGTTTCCAGGTGGGTGCCCTTGCCTTCGTT  >  NC_002947/5891069‑5891356

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