Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,440,630	A→C	6.0%	K125Q (AAG→CAG)	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,630	0	A	C	6.0%	192.4 / 2.0	62	K125Q (AAG→CAG)	PP_3060	tail sheath protein
Reads supporting (aligned to +/- strand):  ref base A (25/33);  new base C (2/2);  total (27/35)
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

GGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCC  >  NC_002947/3440485‑3440771                                                                                                                                                  |                                                                                                                                                ggCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAg                                                                                                                                               >  3:397158/1‑148 (MQ=255)          cAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAACTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTg                                                                                                                                     >  2:277122/1‑149 (MQ=255)           aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGc                                                                                                                                    >  3:389982/1‑149 (MQ=255)            aCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                   <  5:72304/149‑1 (MQ=255)               ccTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCTACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTc                                                                                                                                <  4:397158/149‑1 (MQ=255)                cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGGCAATGCCGTGCTTCAGCTCGGCCTTGTCgccg                                                                                                                               <  1:388729/149‑4 (MQ=38)                     aCGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                          <  5:98529/149‑1 (MQ=255)                               ccGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAACTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 <  1:277122/148‑1 (MQ=255)                                cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                >  1:313762/1‑148 (MQ=255)                                   ggCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGt                                                                                                            >  6:307322/1‑149 (MQ=255)                                    gCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 >  6:148660/1‑143 (MQ=255)                                    gCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                 <  5:148660/143‑1 (MQ=255)                                      tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                >  7:199207/1‑142 (MQ=255)                                      tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCg                                                                                                                <  8:199207/142‑1 (MQ=255)                                       cGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                        <  8:287458/149‑1 (MQ=255)                                        gTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggcg                                                                                                       <  3:281027/149‑1 (MQ=255)                                        gTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCAGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggcg                                                                                                       >  1:322794/1‑149 (MQ=255)                                            aaGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGa                                                                                                   <  5:373542/149‑1 (MQ=255)                                                 ccAGAAGGGCGAAACCCTTGCGGTCGTCGTCACTCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCa                                                                                              >  7:53383/1‑149 (MQ=255)                                                        ggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                       <  1:181763/149‑1 (MQ=255)                                                        ggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                       <  5:307322/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:315978/149‑1 (MQ=255)                                                                                              ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCa                                               >  4:52450/1‑149 (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                                             <  1:303807/124‑1 (MQ=255)                                                                                                                         gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCAcc                                             >  2:303807/1‑124 (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                              >  5:32410/1‑126 (MQ=255)                                                                                                                                      aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGt                              <  6:32410/126‑1 (MQ=255)                                                                                                                                      aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt        >  6:157062/1‑148 (MQ=255)                                                                                                                                         gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc    >  8:237549/1‑149 (MQ=255)                                                                                                                                         gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc    >  3:3218/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)                                                                                                                                                  |                                                                                                                                                GGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGC‑ATGAAGC‑GTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCC  >  NC_002947/3440485‑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: ‑