Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,706,378	Δ1 bp	10.0%	coding (401/924 nt)	PP_1500 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,706,378	0	T	.	10.0%	230.9 / 21.6	80	coding (401/924 nt)	PP_1500	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base . (3/5);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

GGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGCC  >  NC_002947/1706236‑1706525                                                                                                                                               |                                                                                                                                                     ggTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑Tgg                                                                                                                                              <  1:337023/150‑1 (MQ=255)  gTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGc                                                                                                                                             <  2:332064/150‑1 (MQ=255)   tAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCt                                                                                                                                            >  1:219669/1‑150 (MQ=255)     aacTGCCAGGGCACCTGTGGGGCAGCTACTGGCTTCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCA‑GGCCATGGCTgc                                                                                                                                          <  2:792511/148‑1 (MQ=255)      ccTGCCAGGGCTCATTTTGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGCTCAAGTACCGCGGCTACGGCTTTTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCA‑TGGCTgca                                                                                                                                         <  1:1169565/150‑1 (MQ=255)      ccTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCA‑GGCCATGGCTgca                                                                                                                                         <  1:456855/150‑1 (MQ=255)       cTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCA‑GGCCATGGCTgcag                                                                                                                                        >  2:775910/1‑150 (MQ=255)        tGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGAACGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGGACCAGGACATTGCTGCAGCCATGGCC‑TGGCTgcagc                                                                                                                                       >  1:382652/1‑150 (MQ=255)            aGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGgc                                                                                                                                   <  2:735224/150‑1 (MQ=255)               gCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCA‑GGCCATGGCTgcagcaatgtcc                                                                                                                                <  1:775910/150‑7 (MQ=255)                 aCCTGCGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATCGCC‑TGGCTGCAGCAAGCGCCCg                                                                                                                              >  1:1263357/1‑150 (MQ=255)                  ccTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGTCATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGa                                                                                                                             <  1:1309201/150‑1 (MQ=255)                  ccTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGa                                                                                                                             <  1:498384/150‑1 (MQ=255)                   cTGGGGGGTAGCTACTGGGTGGCGGAGCAGGGCTACCAGGTGCTGATGATCGAGGGCCGCGGCTATGGCTTGTCTCAGGGCCGGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            <  1:917949/150‑1 (MQ=255)                   cTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            <  2:1324994/150‑1 (MQ=255)                   cTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            <  1:755193/150‑1 (MQ=255)                   cTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            >  2:269377/1‑150 (MQ=255)                       ggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTgca                                                                                                                                         >  1:59919/1‑133 (MQ=255)                       ggggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTgca                                                                                                                                         <  2:59919/133‑1 (MQ=255)                        gggTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCaa                                                                                                                       <  1:818110/150‑1 (MQ=255)                          gTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAgg                                                                                                                     <  2:869800/150‑1 (MQ=255)                          gTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAAGTGTACCAGGACATTGCTGCAGCCAGGGCC‑TGGCTGCCGCAGGTGCCCGACGCCAAgg                                                                                                                     >  1:675206/1‑150 (MQ=255)                              cTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCaa                                                                                                                 >  2:643576/1‑150 (MQ=255)                                   ggCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGc                                                                                                            <  1:643576/150‑1 (MQ=255)                                    gCTGCCGGAGCAGGGCTACCAGGTGCTGATGCTCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCAGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCt                                                                                                           >  1:1182548/1‑150 (MQ=255)                                         cGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTAc                                                                                                      >  1:735047/1‑150 (MQ=255)                                               aGGGCTACCAGGTGCTGATGATCGATCACCGCGGTTCTGTCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                <  2:1182548/150‑1 (MQ=255)                                                gggCTACCAGGTGCTTATGATCGATTCCCGGCGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTggg                                                                                               <  2:675206/150‑1 (MQ=255)                                                gggCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTggg                                                                                               >  1:840867/1‑150 (MQ=255)                                                gggCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTggg                                                                                               >  2:1211104/1‑150 (MQ=255)                                                 ggCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgggg                                                                                              >  2:908477/1‑150 (MQ=255)                                                   cTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTActgctgagcca                                                                                            <  2:1039127/150‑6 (MQ=255)                                                   cTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCa                                                                                            >  1:657084/1‑150 (MQ=255)                                                      ccAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCAGGCCA‑TGGCTgcagc                                                                                                                                       >  1:1164312/1‑104 (MQ=255)                                                      ccAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCA‑GGCCATGGCTgcagc                                                                                                                                       <  2:1164312/104‑1 (MQ=39)                                                        aGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGAAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGcc                                                                                       >  2:699330/1‑150 (MQ=255)                                                           tGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTgg                                                                                    >  1:895402/1‑150 (MQ=255)                                                             cTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGc                                                                                  <  1:908477/150‑1 (MQ=255)                                                              tgatgaTCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                >  1:44838/1‑135 (MQ=255)                                                              tgatgaTCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                <  2:44838/135‑1 (MQ=255)                                                               gatgatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                <  1:602028/134‑1 (MQ=255)                                                               gatgatCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                >  2:602028/1‑134 (MQ=255)                                                                  gatCGATTACCGCAAGTATAACTTGTCTCACGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGc                                                                             <  2:682872/150‑1 (MQ=255)                                                                     cGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            <  1:1112740/100‑1 (MQ=255)                                                                     cGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAg                                                                                                                            >  2:1112740/1‑100 (MQ=255)                                                                      gATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATg                                                                         <  2:895402/150‑1 (MQ=255)                                                                       aTTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATgg                                                                        <  2:657084/150‑1 (MQ=255)                                                                        ttACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGc                                                                       >  1:430658/1‑150 (MQ=255)                                                                                    tGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcat                                                                    >  2:857016/1‑141 (MQ=255)                                                                                    tGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcat                                                                    <  1:857016/141‑1 (MQ=255)                                                                                       cTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGCGCCTGGGCGGTGCCATGGCCATTCATTACCTggc                                                        >  2:220120/1‑150 (MQ=255)                                                                                       cTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTggc                                                        <  1:1147863/150‑1 (MQ=255)                                                                                       cTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTggc                                                        <  2:840867/150‑1 (MQ=255)                                                                                        ttGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcatt                                                                   <  1:1010791/138‑1 (MQ=255)                                                                                        ttGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCcatt                                                                   >  2:1010791/1‑138 (MQ=255)                                                                                             tcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCa                                                  >  1:215837/1‑150 (MQ=255)                                                                                             tcAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCa                                                  >  2:1113457/1‑150 (MQ=255)                                                                                                  gCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCg                                             >  2:309075/1‑150 (MQ=255)                                                                                                    cAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGag                                           <  1:344745/150‑1 (MQ=255)                                                                                                            ccTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCa                                   <  1:120124/150‑1 (MQ=255)                                                                                                             cTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAg                                  >  2:31545/1‑150 (MQ=255)                                                                                                                ccGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGGGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAgcgc                               <  2:1263357/150‑1 (MQ=255)                                                                                                                        gtACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAgcgc                       >  2:795927/1‑150 (MQ=255)                                                                                                                          aCCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGccc                     >  1:1179773/1‑150 (MQ=255)                                                                                                                           ccAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCt                    <  1:269377/150‑1 (MQ=255)                                                                                                                           ccAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAACAGGCGCCCGAGGTCACGGGCAAGCCGCTGGTACTGCTGGGGTAGAGCCTGTGCGGTGCCATTGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCt                    >  2:520090/1‑150 (MQ=255)                                                                                                                            cAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTg                   <  1:504036/150‑1 (MQ=255)                                                                                                                               gACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                <  2:957707/70‑1 (MQ=255)                                                                                                                               gACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTgg                                                                                                >  1:957707/1‑70 (MQ=255)                                                                                                                                 gcctgctgcAGCCA‑GGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGgt                 >  2:233681/4‑147 (MQ=255)                                                                                                                                 gcctgctgcAGCCA‑GGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGgt                 <  1:233681/144‑1 (MQ=255)                                                                                                                                      ctgcAGCCA‑GGCCATGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATg         >  2:889152/1‑150 (MQ=255)                                                                                                                                          aGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGgtgt     >  2:724104/1‑150 (MQ=255)                                                                                                                                          aGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGgtgt     <  1:31545/150‑1 (MQ=255)                                                                                                                                            ccATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGagcagc                                       >  2:1239374/1‑114 (MQ=255)                                                                                                                                            ccATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGagcagc                                       <  1:1239374/114‑1 (MQ=255)                                                                                                                                            ccATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGc   >  2:693442/1‑150 (MQ=255)                                                                                                                                             cATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGcc  >  1:1368924/1‑150 (MQ=255)                                                                                                                                              aTGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGtt              <  2:352765/137‑1 (MQ=255)                                                                                                                                              aTGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGtt              >  1:352765/1‑137 (MQ=255)                                                                                                                                               |                                                                                                                                                     GGTAACCTGCCAGGGCACCTGGGGGGTAGCTACTGGCTGCCGGAGCAGGGCTACCAGGTGCTGATGATCGATTACCGCGGCTATGGCTTGTCTCAGGGCCAGCCGAGCCTGCCGGAGGTGTACCAGGACATTGCTGCAGCCATGGCC‑TGGCTGCAGCAGGCGCCCGAGGTCAAGGGCAAGCCGCTGGTACTGCTGGGGCAGAGCCTGGGCGGTGCCATGGCCATTCATTACCTGGCGGCCCACCCCGAGCAGCGTCAGCGCTTCAGCGCCCTGGTGTTCGATGGTGTGCC  >  NC_002947/1706236‑1706525

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