Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,491,771	G→C	8.8%	L231V (CTG→GTG)	PP_3096 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,491,771	0	G	C	8.8%	206.0 / 9.4	68	L231V (CTG→GTG)	PP_3096	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base C (3/3);  total (31/37)
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.03e-01

TGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGT  >  NC_002947/3491626‑3491919                                                                                                                                                  |                                                                                                                                                     tGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGc                                                                                                                                                  <  1:577228/150‑1 (MQ=255)  ggAAGCGCTGGTGGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTCGTACCACGCCCGGTTGCCTGCGATGGCTCTGTTGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGca                                                                                                                                                 <  2:215435/150‑1 (MQ=255)  ggAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGca                                                                                                                                                 <  2:145916/150‑1 (MQ=255)      gcgcTGGTGGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGcagcca                                                                                                                                             <  1:1575910/150‑1 (MQ=255)      gcgcTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGcagcca                                                                                                                                             <  1:26821/150‑1 (MQ=255)      gcgcTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGcagcca                                                                                                                                             <  1:716071/150‑1 (MQ=255)        gcTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGcagccagc                                                                                                                                           >  1:818177/1‑150 (MQ=255)         cTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGGTGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCACCTGcagcccgaa                                                                                                                                          <  1:736846/150‑5 (MQ=255)            gTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGTGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGccc                                                                                                                                       >  2:699294/1‑150 (MQ=255)                   gCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAAc                                                                                                                                <  1:54766/150‑1 (MQ=255)                    ccAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACg                                                                                                                               >  2:1034133/1‑150 (MQ=255)                    ccAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACg                                                                                                                               >  1:708335/1‑150 (MQ=255)                       ggCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCt                                                                                                                            <  1:437450/150‑1 (MQ=255)                       ggCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCt                                                                                                                            <  1:764979/150‑1 (MQ=255)                        gCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTg                                                                                                                           >  2:1681995/1‑150 (MQ=255)                           gAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGTGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTAc                                                                                                                        <  1:1121062/150‑1 (MQ=255)                                          gATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCga                                                                                                         >  1:129749/1‑150 (MQ=255)                                              cGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAg                                                                                                     <  1:1189042/150‑1 (MQ=255)                                                gATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACTGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGTGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTa                                                                                                   <  1:548046/150‑1 (MQ=255)                                                gATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTa                                                                                                   >  1:1807212/1‑150 (MQ=255)                                                     aCGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGc                                                                                              >  1:186671/1‑150 (MQ=255)                                                     aCGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGc                                                                                              >  2:1059087/1‑150 (MQ=255)                                                       gACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCTGCTGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCg                                                                                            <  1:20109/150‑1 (MQ=255)                                                       gACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTGCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCg                                                                                            >  2:366043/1‑150 (MQ=255)                                                        aCTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGa                                                                                           <  1:1034133/150‑1 (MQ=255)                                                         cTGGTAGTCCCAGACCCGGTTGCCGGCGACAGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAg                                                                                          >  2:1114828/1‑150 (MQ=255)                                                          tGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCGGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGTATATTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGc                                                                                         <  1:425697/150‑1 (MQ=255)                                                            gTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGCTCGGCAGCGAGAAGTAGAAGCCGAGCAt                                                                                       <  1:1894811/150‑1 (MQ=255)                                                            gTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTa                                                                                                   >  1:1505475/1‑138 (MQ=255)                                                            gTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTa                                                                                                   <  2:1505475/138‑1 (MQ=255)                                                            gTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCAt                                                                                       <  1:22600/150‑1 (MQ=255)                                                            gTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCAt                                                                                       >  1:1605490/1‑150 (MQ=255)                                                              aGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGc                                                                                     <  1:1684558/150‑1 (MQ=255)                                                                  ccAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTc                                                                                 >  1:691873/1‑150 (MQ=255)                                                                     gACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAg                                                                              <  2:691873/150‑1 (MQ=255)                                                                     gACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCAGCTCAACCTCAAAATTTTTTTCTGTAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAg                                                                              <  1:36990/150‑1 (MQ=255)                                                                        ccGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTAc                                                                                                                        >  1:556777/1‑105 (MQ=255)                                                                        ccGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTAc                                                                                                                        <  2:556777/105‑1 (MQ=255)                                                                        ccGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGtt                                                                           >  2:1840706/1‑150 (MQ=255)                                                                        ccGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGtt                                                                           >  2:1253029/1‑150 (MQ=255)                                                                          ggTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTAGCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCa                                                                         >  1:294290/1‑150 (MQ=255)                                                                             tGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGTGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAgcg                                                                      <  1:699294/150‑1 (MQ=255)                                                                               ccGGCGACGGCGCTGTGGCCGGGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCa                                                                    <  2:36994/150‑1 (MQ=255)                                                                                 ggCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCAtt                                                                  >  2:1705860/1‑150 (MQ=255)                                                                                    gACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTgcg                                                               >  2:172353/1‑150 (MQ=255)                                                                                     aCGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTgcgc                                                              <  2:803403/150‑1 (MQ=255)                                                                                       ggCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCgg                                                            <  2:545815/150‑1 (MQ=255)                                                                                        gcgcTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGc                                                           <  1:620356/150‑1 (MQ=255)                                                                                               ggCCGAGCTGGGCATTGCGGCGGCCGCGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGc                                                    <  1:1840706/150‑1 (MQ=255)                                                                                                  cGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGa                                                 <  2:129749/150‑1 (MQ=255)                                                                                                     gCTNGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCa                                              >  2:1203976/1‑150 (MQ=255)                                                                                                       tGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAgg                                            >  1:1152501/1‑150 (MQ=255)                                                                                                         ggCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGcc                                          <  1:1114828/150‑1 (MQ=255)                                                                                                          gCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCg                                         <  2:898584/150‑1 (MQ=255)                                                                                                            aTTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGc                                       >  2:346348/1‑150 (MQ=255)                                                                                                                 ggcggcCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGAGACGAGCTGTACCTTGATCGGCATCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAAt                                  >  1:706409/1‑150 (MQ=255)                                                                                                                                tGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGTTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCac                   <  2:1068426/150‑1 (MQ=255)                                                                                                                                tGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCAAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCac                   >  1:895430/1‑150 (MQ=255)                                                                                                                                  cATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCacga                 <  2:1807212/150‑1 (MQ=255)                                                                                                                                          tggCTGCACCTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTgcgc                                                              <  2:1217838/96‑1 (MQ=255)                                                                                                                                          tggCTGCACCTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTgcgc                                                              >  1:1217838/2‑97 (MQ=255)                                                                                                                                          tggCTGCACCTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGc         >  2:725002/2‑150 (MQ=255)                                                                                                                                            gCTGCACCTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCt                              >  2:1396908/1‑127 (MQ=255)                                                                                                                                            gCTGCACCTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCt                              <  1:1396908/127‑1 (MQ=255)                                                                                                                                             cTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCACGGGTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGa      >  1:1067014/1‑150 (MQ=255)                                                                                                                                               gCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAg                                                                              <  2:1803321/76‑1 (MQ=255)                                                                                                                                               gCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAg                                                                              >  1:1803321/1‑76 (MQ=255)                                                                                                                                                 aGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGt  >  2:1108944/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGGAAGCGCTGGTAGTCGGCCAGGCTGAGCGGGCCCAGCTCGATGCGGATGCACGACTGGTAGTCCCAGACCCGGTTGCCGGCGACGGCGCTGTGGCCGAGCTGGGCATTGCGGCGGCCGAGCTGGGTGCATTGTTCGGGCTGCAGGTGCAGCCAGCGCCCGGCGAACGGCTGTACCTTGATCGGCAGCGAGAAGTAGAAGCCGAGCATGCGTTCAAGGTTCAGCGCATTGCGCGGCTTCTGGCTGAGCAGGCCGGCAAAATAGCTGAACAGCTCACGACGCAATGCTGACGGT  >  NC_002947/3491626‑3491919

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