Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	704,281	A→C	6.2%	Q125P (CAG→CCG)	PP_5745 →	paraquat‑inducible protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	704,281	0	A	C	6.2%	199.8 / 4.2	65	Q125P (CAG→CCG)	PP_5745	paraquat‑inducible protein A
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base C (2/2);  total (33/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 = 9.06e-01

CCCGGACACCATCATGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAGTGTCGAAGCCAACCTGGGC  >  NC_002947/704135‑704426                                                                                                                                                   |                                                                                                                                                    accGGACATGATGGTGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGNCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGc                                                                                                                                                   >  1:109653/2‑149 (MQ=255) accGGACATGATGGTGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGc                                                                                                                                                   >  1:105775/2‑149 (MQ=255)        acCATCATGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCTGCCGCCGCtgaa                                                                                                                                            <  5:227849/149‑4 (MQ=255)         ccATCATGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCg                                                                                                                                           <  3:150677/149‑1 (MQ=255)          catcatGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCgg                                                                                                                                          >  1:210911/1‑149 (MQ=255)           atcatGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGc                                                                                                                                         >  4:12383/1‑149 (MQ=255)           atcatGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGc                                                                                                                                         >  5:67282/1‑149 (MQ=255)            tcatGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGcc                                                                                                                                        >  4:252476/1‑149 (MQ=255)                          atcaCCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATAtt                                                                                                                          >  7:34698/1‑149 (MQ=255)                             aCCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGAt                                                                                                                       <  3:232575/149‑1 (MQ=255)                             aCCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGAt                                                                                                                       <  5:117722/149‑1 (MQ=255)                               ccTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGt                                                                                                                     <  3:12383/149‑1 (MQ=255)                                cTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTa                                                                                                                    <  6:67282/149‑1 (MQ=255)                                   cTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCttt                                                                                                                                             >  5:4995/1‑121 (MQ=255)                                   cTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCttt                                                                                                                                             <  6:4995/121‑1 (MQ=255)                                    tCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGACATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGt                                                                                                                     >  5:303419/1‑144 (MQ=255)                                    tCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGACATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGt                                                                                                                     <  6:303419/144‑1 (MQ=255)                                          aCGGCATGTTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCGACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTGCGTTCAGCGGCGGCAGCCGCTTTCGTCCCGGCAGCGGATATTGATGTACCGCTTCATc                                                                                                          <  1:338320/149‑1 (MQ=255)                                             gCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGaa                                                                                                       <  8:66965/149‑1 (MQ=255)                                                    gCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGc                                                                                                                >  3:272125/1‑133 (MQ=255)                                                    gCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGc                                                                                                                <  4:272125/133‑1 (MQ=255)                                                        aTTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGgcg                                                                                            >  5:272099/1‑149 (MQ=255)                                                        aTTGCGGCTGTGGTGTTAATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGgcg                                                                                            >  1:79158/1‑149 (MQ=255)                                                                    gtgtTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGc                                                                                 <  2:210911/148‑1 (MQ=255)                                                                     tgtTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTg                                                                               >  5:247483/1‑149 (MQ=255)                                                                      gtTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTg                                                                               >  7:216016/1‑148 (MQ=255)                                                                       ttCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTg                                                                                                >  7:313966/1‑130 (MQ=255)                                                                       ttCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTg                                                                                                <  8:313966/130‑1 (MQ=255)                                                                       ttCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑a                                                                           <  2:109653/149‑1 (MQ=255)                                                                       ttCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑a                                                                           <  2:105775/149‑1 (MQ=255)                                                                              ccAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTc                                                                    <  6:247483/149‑1 (MQ=255)                                                                                  tATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCa                                                                <  3:252476/149‑1 (MQ=255)                                                                                     ccTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTg                                                             >  8:248335/1‑149 (MQ=255)                                                                                         gTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGcc                                                           <  8:216016/147‑1 (MQ=255)                                                                                            cccAGGTTCAAGCTGGTGTGCTTCGGGTTGGTGTTGTAGTCCGTTCTGCGGCGGCAGCCGCTTTCGGCCCGGGAGCGGATATTGATGTAGCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCAtt                                                              <  2:152094/141‑1 (MQ=255)                                                                                            cccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGGTTCATCGAATTCATTGGGCGCTGGTCCATGCTCG‑‑ACATCTTCGTCAtt                                                              >  1:152094/1‑141 (MQ=255)                                                                                            cccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCAtt                                                              <  2:146729/141‑1 (MQ=255)                                                                                            cccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCAtt                                                              >  1:146729/1‑141 (MQ=255)                                                                                            cccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCt                                                      <  2:330595/149‑1 (MQ=255)                                                                                            cccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCt                                                      >  1:330595/1‑149 (MQ=255)                                                                                             ccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTCCCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtg                                                     >  8:316112/1‑149 (MQ=255)                                                                                             ccACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtg                                                     >  1:135502/1‑149 (MQ=255)                                                                                             ccACGTACAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGTCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTCTTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtg                                                     >  6:338612/1‑149 (MQ=255)                                                                                              cACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtgg                                                    >  5:76367/1‑149 (MQ=255)                                                                                              cACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtgg                                                    >  7:284935/1‑149 (MQ=255)                                                                                               aCGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtggt                                                   <  7:316112/149‑1 (MQ=255)                                                                                                     aaGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAg                                              <  6:76367/148‑1 (MQ=255)                                                                                                      aGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑Ac                                                                          <  4:306183/119‑1 (MQ=255)                                                                                                      aGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑Ac                                                                          >  3:306183/1‑119 (MQ=255)                                                                                                               gCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtggtg                                                  <  5:70618/134‑1 (MQ=255)                                                                                                               gCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtggtg                                                  >  6:70618/1‑134 (MQ=255)                                                                                                                cATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTg                                  >  2:318118/1‑149 (MQ=255)                                                                                                                    ggCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATTCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGcc                               <  8:34698/148‑1 (MQ=255)                                                                                                                          cTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGc                        <  7:248335/149‑1 (MQ=255)                                                                                                                           tgttgtACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCtggt                                                   >  2:190698/1‑121 (MQ=255)                                                                                                                               gtACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGAAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCCCAGGAATAAGCTGGATAGATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAg                     >  1:112133/1‑149 (MQ=255)                                                                                                                                tACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAgtg                   <  8:284935/148‑1 (MQ=255)                                                                                                                                tACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAgtg                   <  6:272099/148‑1 (MQ=255)                                                                                                                                   tCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAGTGTCg                >  6:31429/1‑148 (MQ=255)                                                                                                                                    gccgaaaGCGGCTGCCGCCGCTTTCGGCGCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAGTGTCg                >  1:196530/7‑147 (MQ=255)                                                                                                                                       ttCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAAt                           <  6:5105/133‑1 (MQ=255)                                                                                                                                       ttCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAAt                           >  5:5105/1‑133 (MQ=255)                                                                                                                                        aaaGCGGCTGCCGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAgtggt                                          >  6:66166/3‑117 (MQ=255)                                                                                                                                        aaaGCGGCTGCCGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAgtggt                                          <  5:66166/115‑1 (MQ=255)                                                                                                                                                ggcAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAGTGTCGAAGCCAACCTGGgc  <  5:31429/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCCGGACACCATCATGTCCGGTGTCATCACCCTGCTCAAGCACGGCATGGTGCCCATTGCGGCTGTGGTGTTCATTGCCAGTATCCTGGTCCCCACGTTCAAGCTGGTGGGCATCGGCTTGCTGTTGTACTCCGTTCAGCGGCGGCAGCCGCTTTCGGCCCGGCAGCGGATATTGATGTACCGCTTCATCGAATTCATTGGGCGCTGGTCCATGCTGG‑‑ACATCTTCGTCATTGCCATCCTGGTGGCAGTGGTGAATTTTGGCCGAATAGCCAGTGTCGAAGCCAACCTGGGC  >  NC_002947/704135‑704426

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