Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,158,426	A→G	9.2%	intergenic (+92/‑15)	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,426	0	A	G	9.2%	181.7 / 7.3	68	intergenic (+92/‑15)	PP_5401/PP_5739	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base G (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 = 5.06e-01

CGGGCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCA  >  NC_002947/6158277‑6158565                                                                                                                                                       |                                                                                                                                            cGGGCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCCCCCACTGGCTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCgca                                                                                                                                             <  1:118185/150‑1 (MQ=255)  gggCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCgcag                                                                                                                                            >  1:179055/1‑150 (MQ=255)    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:1172236/1‑150 (MQ=255)                               gatTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCg                                                                                                               >  1:1926813/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                                                                                                         <  1:33218/139‑1 (MQ=255)                                                cTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGa                                                                                                         >  2:33218/1‑139 (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                                                                                                   <  2:608162/108‑1 (MQ=255)                                                                                     cAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGc                                                                                                   >  1:608162/1‑108 (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                                                    <  2:55863/149‑1 (MQ=255)                                                                                           gtACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTagc                                                    >  1:55863/1‑149 (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                                                             <  1:711561/134‑1 (MQ=255)                                                                                                 attgatCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGaa                                                             >  2:711561/1‑134 (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                          >  1:1449496‑M2/1‑18 (MQ=255)                                                                                                                                        aCCCTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcctggtagctactcagtcaaaaggtgggccattattaa                          <  2:1449496‑M2/130‑113 (MQ=255)                                                                                                                                          ccTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcct                                                             <  1:642188‑M2/93‑78 (MQ=255)                                                                                                                                          ccTGCCG‑GCCGCAGCAcggcagggttcttccaaggattgtgcagcgtcgcgcttactgatcaatcagtacaccttgagaacccgaaagcgcct                                                             >  2:642188‑M2/1‑16 (MQ=255)                                                                                                                                           cTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTc   >  1:561106/1‑150 (MQ=255)                                                                                                                                            tGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGCGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCa  <  1:558449/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                            CGGGCACGTTCGCAAATTGCTGGGAGAGATGATTAATAATGGCCCACCTTTTGACTGAGTAGCTACCAGGCGCTTTCGGGTTCTCAAGGTGTACTGATTGATCAGTAAGCGCGACGCTGCACAATCCTTGGAAGAACCCTGCCG‑GCCGCAGCACGGCAGGGTCGATGCAACTATGGGTCGAGGCGATCAAGCATATCCAGCACGCCGCTAAACAAGCGACCTACGCTGAACACACTAGCGAGCGCACACAAGCCTCGCCAGATACGGCAAAGCATAAGCTTCCAGCTCA  >  NC_002947/6158277‑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: ‑