Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,169,785	A→G	6.0%	H452R (CAC→CGC)	PP_4547 →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,169,785	0	A	G	6.0%	196.0 / 2.3	67	H452R (CAC→CGC)	PP_4547	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

GACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGATGTGG  >  NC_002947/5169639‑5169934                                                                                                                                                   |                                                                                                                                                       gACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCa                                                                                                                                                     <  2:398930/150‑1 (MQ=255)    cGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgtc                                                                                                                                                  <  2:1296002/150‑3 (MQ=255)    cGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                                  >  1:377407/1‑150 (MQ=255)        cTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCCGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGa                                                                                                                                              >  2:15494/1‑150 (MQ=255)                 cTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTaaa                                                                                                                                     >  1:263352/1‑150 (MQ=255)                       cATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACaa                                                                                                                              <  2:931720/150‑1 (MQ=255)                       cATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCCGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACaa                                                                                                                              <  1:15494/150‑1 (MQ=255)                         ttCCGATTGGGCCCGTGAAGGATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAgg                                                                                                                            >  1:440118/1‑150 (MQ=255)                         ttCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAgg                                                                                                                            <  2:624747/150‑1 (MQ=255)                          tacGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCGGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGa                                                                                                                           >  1:289615/3‑150 (MQ=255)                              aTTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAacac                                                                                                                       <  2:263352/150‑1 (MQ=255)                                     ccGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGaa                                                                                                                >  1:237745/1‑150 (MQ=255)                                        tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                               >  1:10135/1‑148 (MQ=255)                                        tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                               <  2:10135/148‑1 (MQ=255)                                        tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACgc                                                                                                             >  2:878975/1‑150 (MQ=255)                                              aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                                    <  2:122000/122‑1 (MQ=255)                                              aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                                    >  1:122000/1‑122 (MQ=255)                                              aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCAtt                                                                                                       <  2:237745/150‑1 (MQ=255)                                                tGGGCTGAGACTGGCTGGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGa                                                                                                     <  2:658757/150‑1 (MQ=255)                                                  ggCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGaac                                                                                                   >  1:942088/1‑150 (MQ=255)                                                     gAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGaacaac                                                                                                <  1:1039064/150‑1 (MQ=255)                                                              cTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                                  <  1:1110437/92‑1 (MQ=255)                                                              cTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                                  >  2:1110437/1‑92 (MQ=255)                                                                    gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgt                                                                                                                                                   >  2:884188/1‑84 (MQ=25)                                                                    gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgt                                                                                                                                                   <  1:884188/85‑2 (MQ=37)                                                                    gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAgcgc                                                                                 >  2:1337845/1‑150 (MQ=255)                                                                     tgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAgcgcg                                                                                >  2:1178279/1‑150 (MQ=255)                                                                          cTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGcc                                                                                                                                                 >  1:529483/1‑81 (MQ=255)                                                                          cTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGcc                                                                                                                                                 <  2:529483/81‑1 (MQ=255)                                                                           tGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGt                                                                          <  2:440118/150‑1 (MQ=255)                                                                           tGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCAGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAAACTACG‑ACATCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGt                                                                          >  1:1156085/1‑150 (MQ=255)                                                                            gAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTc                                                                         <  1:1178279/150‑1 (MQ=255)                                                                                gTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCctac                                                                     <  1:416991/150‑1 (MQ=255)                                                                                    aGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                               >  2:884152/1‑104 (MQ=255)                                                                                    aGCGCGCCGAGTACCGCCAGTTCATGGCAGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                               <  1:884152/104‑1 (MQ=255)                                                                                      cgcgcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAgc                                                               <  2:377407/150‑1 (MQ=255)                                                                                         gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGcc                                                            >  1:1197156/1‑150 (MQ=255)                                                                                         gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGcc                                                            >  1:1197863/1‑150 (MQ=255)                                                                                                  cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                                   >  1:606344/1‑150 (MQ=255)                                                                                                  cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                                   <  1:981082/150‑1 (MQ=255)                                                                                                  cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                                   >  2:876118/1‑150 (MQ=255)                                                                                                   gCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACg                                                  >  2:887143/1‑150 (MQ=255)                                                                                                   gCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACg                                                  >  1:1075286/1‑150 (MQ=255)                                                                                                    ccAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGa                                                 >  1:782670/1‑150 (MQ=255)                                                                                                    ccAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCAACAGGCCTGAACCCCTAAAT‑AAGGAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGa                                                 <  1:878975/150‑1 (MQ=255)                                                                                                     cAGCTCATGGCCGAAGGCAGCGAACAGCACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACCCCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGAc                                                <  2:145404/150‑1 (MQ=255)                                                                                                       gTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCAGGCGCAGCGCGCGTCGTCCTACTACAGCGCCACCCTCCACGACca                                              <  1:1176574/150‑1 (MQ=255)                                                                                                             ggCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGc                                        <  2:782670/150‑1 (MQ=255)                                                                                                             ggCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGc                                        <  2:956658/150‑1 (MQ=255)                                                                                                                cGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTa                                     <  2:289615/150‑1 (MQ=255)                                                                                                                        gCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCa                                                     <  2:505230/126‑1 (MQ=255)                                                                                                                        gCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCa                                                     >  1:505230/1‑126 (MQ=255)                                                                                                                         cGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACAGCGTACCCGCAACt                            >  2:369229/1‑150 (MQ=255)                                                                                                                          gAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTc                           >  2:1326321/1‑150 (MQ=255)                                                                                                                              aGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCAACCTCAACGACCCCCCTGCGTACCCGCAACTCAAgg                       <  2:1156085/150‑1 (MQ=255)                                                                                                                                gccTGGTGCAGGTACCAGCGCCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGc                     >  1:1353054/3‑150 (MQ=255)                                                                                                                                         tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg            <  1:1326321/150‑1 (MQ=255)                                                                                                                                         tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg            <  1:887143/150‑1 (MQ=255)                                                                                                                                         tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg            >  1:64770/1‑150 (MQ=255)                                                                                                                                           gTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGACGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGt          >  1:1017649/1‑150 (MQ=255)                                                                                                                                             aCCTGCAACAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTgg        >  2:1421429/1‑150 (MQ=255)                                                                                                                                              ccTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                                   >  2:345530/1‑106 (MQ=255)                                                                                                                                              ccTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                                   <  1:345530/106‑1 (MQ=255)                                                                                                                                               cTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg            <  2:746008/144‑1 (MQ=255)                                                                                                                                               cTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg            >  1:746008/1‑144 (MQ=255)                                                                                                                                               cTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGAt      >  2:505954/1‑150 (MQ=255)                                                                                                                                               cTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGAt      >  1:1254869/1‑150 (MQ=255)                                                                                                                                                   aCCAGGCCTGATCCCCTAAAT‑AACAAGGAACACCCATCAACGCAGCATTGAACAACGGCCCGGCGCATCGCGCGCCGTCCTACTACAGCGCCACCCTCAATGTCCACACTTCGTACCCTCAACTCAAGGGCACGGTGCAGGTGGATGTgg  <  1:574459/150‑1 (MQ=255)                                                                                                                                                   aCCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAgg           <  2:419461/141‑1 (MQ=255)                                                                                                                                                   aCCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCCCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAgg           >  1:419461/1‑141 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAGGTGGATGTGG  >  NC_002947/5169639‑5169934

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