Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,973,259	Δ1 bp	6.5%	coding (615/1713 nt)	PP_2600 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,973,259	0	A	.	6.5%	185.8 / 6.9	62	coding (615/1713 nt)	PP_2600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base . (2/2);  total (30/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 = 1.00e+00

GCGTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTC  >  NC_002947/2973111‑2973406                                                                                                                                                     |                                                                                                                                                     gCGTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGa                                                                                                                                                      >  2:63133/1‑149 (MQ=255)            cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGCCCATGG‑CGGGc                                                                                                                                          >  1:285330/1‑149 (MQ=255)            cACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGc                                                                                                                                          >  1:197282/1‑149 (MQ=255)             aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cggg                                                                                                                                           >  6:75034/1‑147 (MQ=255)             aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cggg                                                                                                                                           <  5:75034/147‑1 (MQ=255)             aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGcc                                                                                                                                         <  1:283613/149‑1 (MQ=255)             aCGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGc                                                                                                                                          >  2:191573/1‑148 (MQ=255)                 cccGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggt                                                                                                                                     >  1:308458/1‑149 (MQ=255)                   cgcgCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtgg                                                                                                                                   <  1:337128/149‑1 (MQ=255)                   cgcgCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtgg                                                                                                                                   >  2:155241/1‑149 (MQ=255)                    gcgcACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAggtgg                                                                                                                                   >  5:159429/1‑148 (MQ=255)                        aCCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAg                                                                                                                              >  6:359553/1‑149 (MQ=255)                        aCCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAg                                                                                                                              <  2:28473/149‑1 (MQ=255)                          cTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGt                                                                                                                            >  5:68179/1‑149 (MQ=255)                                  ccTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGcc                                                                                                                     >  2:144230/1‑148 (MQ=255)                                      gATCAGCAGGCGTTTGGCCGGCGACGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAg                                                                                                                 <  6:320465/148‑1 (MQ=255)                                       aTCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAgcg                                                                                                               <  6:66566/149‑1 (MQ=255)                                        tCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAgcgc                                                                                                              <  8:228027/149‑1 (MQ=255)                                         cagcagGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCt                                                                                                             >  4:135874/1‑149 (MQ=255)                                                  tttGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGt                                                                                                    <  6:190859/149‑1 (MQ=255)                                                    tGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTg                                                                                                   <  5:359553/148‑1 (MQ=255)                                                          gcgcCGACGCTTGCAGGTTCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAgg                                                                                            <  8:83955/149‑1 (MQ=255)                                                              cGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCtt                                                                                        <  1:155241/149‑1 (MQ=255)                                                               gACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTg                                                                                                           >  4:297516/1‑129 (MQ=255)                                                               gACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTg                                                                                                           <  3:297516/129‑1 (MQ=255)                                                                 cGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTgg                                                                                     <  1:191573/149‑1 (MQ=255)                                                                   cTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGc                                                                                   <  2:308458/149‑1 (MQ=255)                                                                   cTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGc                                                                                   <  1:63133/149‑1 (MQ=255)                                                                    ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cgg                                                                                                                                            <  4:249997/91‑1 (MQ=255)                                                                    ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑Cgg                                                                                                                                            >  3:249997/1‑91 (MQ=255)                                                                    ttGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTAAAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCg                                                                                  >  1:278303/1‑149 (MQ=255)                                                                         ggTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGGGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCgg                                                                             <  1:176199/149‑1 (MQ=255)                                                                           tGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGt                                                                           <  8:217285/149‑1 (MQ=255)                                                                              acacCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCAt                                                                        <  2:197282/149‑1 (MQ=255)                                                                               cacCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATg                                                                       <  6:159429/149‑1 (MQ=255)                                                                                acCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGt                                                                      >  4:29269/1‑149 (MQ=255)                                                                                acCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGt                                                                      >  5:77590/1‑149 (MQ=255)                                                                                       tCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGaa                                                               <  6:77590/149‑1 (MQ=255)                                                                                          tgtgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGtt                                                                     >  6:283486/1‑140 (MQ=255)                                                                                          tgtgCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGtt                                                                     <  5:283486/140‑1 (MQ=255)                                                                                              ccGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAACGCGCTTGCCCACGTGCAGCAGGGCTGTGGGCGAGCGGGTCAGGTTGAGGAAggcggc                                                         >  7:263853/1‑148 (MQ=255)                                                                                                 aTTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAgg                                                     >  3:320553/1‑149 (MQ=255)                                                                                                     aGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCg                                                  <  2:285330/148‑1 (MQ=255)                                                                                                      gCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCgg                                                 <  3:135874/148‑1 (MQ=255)                                                                                                             gcgcCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGc                                          >  7:252387/1‑148 (MQ=255)                                                                                                                          cacaCCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGc                            >  3:255371/1‑149 (MQ=255)                                                                                                                               cGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCAt                        >  1:232458/1‑148 (MQ=255)                                                                                                                               cGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCAt                        >  3:8027/1‑148 (MQ=255)                                                                                                                                 ttGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATc                                   >  2:257779/1‑135 (MQ=255)                                                                                                                                 ttGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATc                                   <  1:257779/135‑1 (MQ=255)                                                                                                                                     ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAg                 <  1:144230/149‑1 (MQ=255)                                                                                                                                     ccAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAg                 >  3:35039/1‑149 (MQ=255)                                                                                                                                      cAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAgg                >  5:152506/1‑149 (MQ=255)                                                                                                                                          tGCATGGCCCGCCCATGG‑CGGGCCAGGTGGTAGCGGTGGCGGGCTAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGGACCGGCCATGATCGGCATGCCCATGAACGAGGTTTc            <  4:256170/149‑1 (MQ=255)                                                                                                                                           acctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  <  6:87917/139‑1 (MQ=255)                                                                                                                                           acctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGa                  >  5:87917/4‑142 (MQ=255)                                                                                                                                           acctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGt               <  6:369745/142‑1 (MQ=255)                                                                                                                                           acctGGCCCG‑CCATGGTCGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGt               >  5:369745/4‑145 (MQ=255)                                                                                                                                             aTGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTcgcg         >  5:156673/1‑149 (MQ=255)                                                                                                                                                  ccGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGTAt                              >  2:265819/1‑123 (MQ=255)                                                                                                                                                  ccGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGTAt                              <  1:265819/123‑1 (MQ=255)                                                                                                                                                    gACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTc  <  6:68179/149‑1 (MQ=255)                                                                                                                                                    gACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTc  <  8:252387/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGTCGCGACCCACGCCCCGCGCACCTGGCCACCCTCGATCAGCAGGCGTTTGGCCGGCGCCGACGCTTGCAGGTGCACACCCAGGTCCTGTGCCGATTTAGCCAGGCGCGCCACCAGTGCCACACCGTTGACCAGGTGCATGGCCCGACCATGG‑CGGGCCAGGTGGTACAGGTGGCGGCCAAAGCGCTTGCACACGTGCAGCAGGGCTTTGGGCGAGCGGGTCATGTTGAGGAAGGCGGCAAGGTCGGCACCGGCCATGATCGGCATGCCCATGAACGAGGTTTCGCGCATGGTC  >  NC_002947/2973111‑2973406

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