Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,158,429	Δ1 bp	8.8%	intergenic (+95/‑12)	PP_5401 → / → PP_5739	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,158,429	0	A	.	8.8%	204.6 / 14.4	68	intergenic (+95/‑12)	PP_5401/PP_5739	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base . (3/3);  total (34/34)
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.34e-01

GCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCA  >  NC_002947/6158280‑6158565                                                                                                                                                       |                                                                                                                                         gCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCACGGTGTACTGATTGATCAGTAAGCGCGACGCTTCACAATCCTTGGAAGAACCCTGCCG‑GCCgcagca                                                                                                                                          >  2:187881/1‑150 (MQ=255)        cGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCAcggcagg                                                                                                                                   >  1:698742‑M2/1‑143 (MQ=255)        cGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCAcggcagg                                                                                                                                   >  2:148032‑M2/1‑143 (MQ=255)              ttGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGAt                                                                                                                             >  1:491577/1‑150 (MQ=255)              ttGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGAt                                                                                                                             <  2:1693240/150‑1 (MQ=255)                    ggAGAGATGATTAATAATGGCCCACCTTTTTACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACt                                                                                                                       <  1:267114/150‑1 (MQ=255)                    ggAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACt                                                                                                                       >  1:300092/1‑150 (MQ=255)                      agagaTGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTAt                                                                                                                     <  1:1803283/150‑1 (MQ=255)                         gatgatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATggg                                                                                                                  >  1:816675/1‑150 (MQ=255)                          atgatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGt                                                                                                                 >  1:606136/1‑150 (MQ=255)                          atgatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGt                                                                                                                 >  2:1285163/1‑150 (MQ=255)                            gatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCg                                                                                                               >  1:1926813/1‑150 (MQ=255)                            gatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCg                                                                                                               >  1:1172236/1‑150 (MQ=255)                             atTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGAt                                                                                                                             <  1:787816/135‑1 (MQ=255)                             atTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGAt                                                                                                                             >  2:787816/1‑135 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCTGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAgg                                                                                                            <  1:1191544/150‑1 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTACGAGCGCCGCCGAACAATAAATCGCACAACTCGGACA‑GCTCCTGCGCGGCAGGGTCGATGCAACTATGGGTCGAgg                                                                                                            <  1:529704/150‑1 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttcttccaaggattgtgcagc                                                                                                               >  1:1199754‑M2/1‑121 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttcttccaaggattgtgcagc                                                                                                               <  2:1199754‑M2/149‑29 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAgg                                                                                                            <  2:1926813/150‑1 (MQ=255)                               taataaTGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAgg                                                                                                            <  2:1101030/150‑1 (MQ=255)                                  taatGGCCCACCTTTTGACTGAGTAGCTACCAGGCGATTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGa                                                                                                         >  2:1202560/1‑150 (MQ=255)                                      ggCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCaa                                                                                                     <  2:179055/150‑1 (MQ=255)                                           aCCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCata                                                                                                >  1:751797/1‑150 (MQ=255)                                             cTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGa                                                                                                         >  2:33218/1‑139 (MQ=255)                                             cTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGa                                                                                                         <  1:33218/139‑1 (MQ=255)                                                             tACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAAc                                                                              <  2:1172236/150‑1 (MQ=255)                                                             tACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAAc                                                                              <  2:1061780/150‑1 (MQ=255)                                                              aCCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACa                                                                             >  1:1649333/1‑150 (MQ=255)                                                               ccAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACaa                                                                            >  2:578378/1‑150 (MQ=255)                                                                cAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAg                                                                           >  1:389837/1‑150 (MQ=255)                                                                  ggCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCg                                                                         >  2:1156626/1‑150 (MQ=255)                                                                   gcgcTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGa                                                                        <  2:606136/150‑1 (MQ=255)                                                                       tttCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttctt                                                                                                                               <  2:806176‑M2/93‑13 (MQ=255)                                                                       tttCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttctt                                                                                                                               >  1:806176‑M2/1‑81 (MQ=255)                                                                             gTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttcttccaaggatt                                                                                                                      <  2:1878836‑M2/96‑22 (MQ=255)                                                                             gTTCTCAAGGTGTACAGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCGTGCTGCGGCggcagggttcttccaaggatt                                                                                                                      >  1:1878836‑M2/1‑75 (MQ=255)                                                                              ttCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGaa                                                             <  2:491577/150‑1 (MQ=255)                                                                                  cAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGc                                                                                                   >  1:608162/1‑108 (MQ=255)                                                                                  cAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGc                                                                                                   <  2:608162/108‑1 (MQ=255)                                                                                  cAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAacaca                                                         <  1:1370980/150‑1 (MQ=255)                                                                                  cAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACAAGGCCGAACAAGCGACCTACGCTGAacaca                                                         <  1:1508044/150‑1 (MQ=255)                                                                                     ggTGTACTGATTGATCAGTAAGCGCGACGCTTCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTa                                                      <  1:187881/150‑1 (MQ=255)                                                                                        gtACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagc                                                    >  1:55863/1‑149 (MQ=255)                                                                                        gtACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagc                                                    <  2:55863/149‑1 (MQ=255)                                                                                         tACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagcga                                                  >  2:1017806/1‑150 (MQ=255)                                                                                          aCTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagcgag                                                 <  1:1017806/150‑1 (MQ=255)                                                                                            tgattgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagcgagcg                                               >  1:1510129/1‑150 (MQ=255)                                                                                              attgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGaa                                                             >  2:711561/1‑134 (MQ=255)                                                                                              attgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGaa                                                             <  1:711561/134‑1 (MQ=255)                                                                                              attgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGca                                             <  2:389837/150‑1 (MQ=255)                                                                                              attgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGca                                             >  1:1613158/1‑150 (MQ=255)                                                                                               ttgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGcac                                            >  2:558449/1‑150 (MQ=255)                                                                                                  atCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGAACTACGCTGAACACACTAGCGAGCGcacaca                                         >  2:1416433/1‑150 (MQ=255)                                                                                                   tCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACaa                                        <  2:1649333/150‑1 (MQ=255)                                                                                                                 gCTGCACACTCCTTGGACGACCCCCGACA‑TCCCCCTCACTCCATAGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATAc                          <  2:1667104/150‑1 (MQ=255)                                                                                                                  cTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACg                         <  1:1156626/150‑1 (MQ=255)                                                                                                                   tGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACgg                        <  1:1285163/150‑1 (MQ=255)                                                                                                                         aTCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGc                  <  2:300092/150‑1 (MQ=255)                                                                                                                             ttGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATaa              >  1:45052/1‑150 (MQ=255)                                                                                                                             ttGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATaa              >  1:871796/1‑150 (MQ=255)                                                                                                                                   gaaCCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTcc        >  1:1899018/1‑150 (MQ=255)                                                                                                                                    aaCCCTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcctggtagctactcagtcaaaaggtgggccattattaatcatctctcccagcaattt       <  2:1690321‑M2/150‑132 (MQ=255)                                                                                                                                     aCCCTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcctggtagctactcagtcaaaaggtgggccattattaa                          <  2:1449496‑M2/130‑113 (MQ=255)                                                                                                                                     aCCCTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcctggtagctactcagtcaaaaggtgggccattattaa                          >  1:1449496‑M2/1‑18 (MQ=255)                                                                                                                                       ccTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcct                                                             >  2:642188‑M2/1‑16 (MQ=255)                                                                                                                                       ccTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcct                                                             <  1:642188‑M2/93‑78 (MQ=255)                                                                                                                                        cTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTc   >  1:561106/1‑150 (MQ=255)                                                                                                                                         tGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGCGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCa  <  1:558449/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                         GCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCA  >  NC_002947/6158280‑6158565

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