Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,440,625	A→T	7.3%	E123V (GAG→GTG)	PP_3060 →	tail sheath protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,440,625	0	A	T	7.3%	182.9 / 3.1	63	E123V (GAG→GTG)	PP_3060	tail sheath protein
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base T (2/3);  total (27/36)
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.23e-01

GGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCC  >  NC_002947/3440480‑3440771                                                                                                                                                   |                                                                                                                                                     ggCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTATTAAGAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAg                                                                                                                                                     <  1:339842/148‑1 (MQ=255)      ggCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAg                                                                                                                                               >  3:397158/1‑148 (MQ=255)               cAACGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAACTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTg                                                                                                                                     >  2:277122/1‑149 (MQ=255)                aaCGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGc                                                                                                                                    >  3:389982/1‑149 (MQ=255)                 aCGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                   <  5:72304/149‑1 (MQ=255)                    ccTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCTACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTc                                                                                                                                <  4:397158/149‑1 (MQ=255)                     cTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGGCAATGCCGTGCTTCAGCTCGGCCTTGTCgccg                                                                                                                               <  1:388729/149‑4 (MQ=38)                          aCGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                          <  5:98529/149‑1 (MQ=255)                                    ccGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAACTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 <  1:277122/148‑1 (MQ=255)                                     cGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                >  1:313762/1‑148 (MQ=255)                                        ggCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGt                                                                                                            >  6:307322/1‑149 (MQ=255)                                         gCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 <  5:148660/143‑1 (MQ=255)                                         gCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 >  6:148660/1‑143 (MQ=255)                                           tCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                <  8:199207/142‑1 (MQ=255)                                           tCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                >  7:199207/1‑142 (MQ=255)                                            cGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                        <  8:287458/149‑1 (MQ=255)                                             gT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggcg                                                                                                       <  3:281027/149‑1 (MQ=255)                                             gT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCAGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggcg                                                                                                       >  1:322794/1‑149 (MQ=255)                                                  aaGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGa                                                                                                   <  5:373542/149‑1 (MQ=255)                                                       ccAGAAGGGCGAAACCCTTGCGGTCGTCGTCACTCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCa                                                                                              >  7:53383/1‑149 (MQ=255)                                                              ggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                       <  5:307322/149‑1 (MQ=255)                                                              ggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                       <  1:181763/149‑1 (MQ=255)                                                                cGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACgaga                                                                                     <  3:300095/149‑1 (MQ=255)                                                                 gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAt                                                                                    <  6:202570/149‑1 (MQ=255)                                                                 gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAt                                                                                    >  6:265300/1‑149 (MQ=255)                                                                       cTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTTATCTACGAGATTGAccc                                                                              >  6:114521/1‑149 (MQ=255)                                                                             gtcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCgg                                                                        <  2:81334/149‑1 (MQ=255)                                                                                gtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑At                                                                    <  4:328616/149‑1 (MQ=255)                                                                                   gtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTc                                                                                               <  5:281841/120‑1 (MQ=255)                                                                                   gtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTc                                                                                               >  6:281841/1‑120 (MQ=255)                                                                                       cccTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑gt                                                            <  5:265300/149‑1 (MQ=255)                                                                                               ggCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCaa                                                    <  8:53383/149‑1 (MQ=255)                                                                                                    ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCa                                               >  4:52450/1‑149 (MQ=255)                                                                                                    ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCa                                               <  4:315978/149‑1 (MQ=255)                                                                                                      tGAAAGAGCTGGACATGGGCGACTGTAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcg                                                                                                          <  7:268454/90‑1 (MQ=255)                                                                                                      tGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcg                                                                                                          >  8:268454/1‑90 (MQ=255)                                                                                                      tGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCAcc                                             >  3:275957/1‑149 (MQ=255)                                                                                                         aaGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGa                                           <  2:313762/148‑1 (MQ=255)                                                                                                         aaGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGAc                                          <  6:54418/149‑1 (MQ=255)                                                                                                                      gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAAc                                                   >  3:14532/1‑127 (MQ=255)                                                                                                                      gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAAc                                                   <  4:14532/127‑1 (MQ=255)                                                                                                                       ggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTTATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTc                             <  5:114521/148‑1 (MQ=255)                                                                                                                        gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAg                           <  4:275957/149‑1 (MQ=255)                                                                                                                        gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAg                           >  1:401727/1‑149 (MQ=255)                                                                                                                               gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCAcc                                             >  2:303807/1‑124 (MQ=255)                                                                                                                               gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCAcc                                             <  1:303807/124‑1 (MQ=255)                                                                                                                               gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                     <  3:52450/148‑1 (MQ=255)                                                                                                                               gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCg                    <  4:389982/149‑1 (MQ=255)                                                                                                                                     cGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCg              >  2:315489/1‑149 (MQ=255)                                                                                                                                         ggcAATGCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc            >  3:45981/3‑147 (MQ=255)                                                                                                                                         ggcAATGCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc            <  4:45981/145‑1 (MQ=255)                                                                                                                                          aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg         >  4:376302/1‑149 (MQ=255)                                                                                                                                           cAAGGCCGTGCTGAAGCACGGCATTGCCGTCACCTAAAACAACCGCGAAGTCGGCGGCGTGGTAATCTACGCTATTGCCCGGGTCGCCTATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg         <  1:115825/149‑1 (MQ=255)                                                                                                                                           cAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg         >  4:397481/1‑148 (MQ=255)                                                                                                                                            aatgCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt        >  6:251914/4‑148 (MQ=255)                                                                                                                                            aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGt                              <  6:32410/126‑1 (MQ=255)                                                                                                                                            aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGt                              >  5:32410/1‑126 (MQ=255)                                                                                                                                            aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt        >  6:157062/1‑148 (MQ=255)                                                                                                                                               gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc    >  3:3218/1‑149 (MQ=255)                                                                                                                                               gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc    >  8:237549/1‑149 (MQ=255)                                                                                                                                               gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc    <  3:397481/149‑1 (MQ=255)                                                                                                                                                ccGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTAATCTACGAGATTGACCCGGTCGGC‑ATGAAGCTATGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAc     <  1:315489/148‑1 (MQ=255)                                                                                                                                                 cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcccc  <  4:3218/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGT‑TCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCC  >  NC_002947/3440480‑3440771

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