Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,784	A→T	6.4%	I398L (ATA→TTA)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,784	0	A	T	6.4%	192.1 / 2.2	62	I398L (ATA→TTA)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base T (2/2);  total (26/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.80e-01

GCTGATCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTTGATGATGCT  >  NC_002947/3686636‑3686931                                                                                                                                                     |                                                                                                                                                    gCTGATCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGata                                                                                                                                                     >  2:46331/1‑149 (MQ=255)   tGATCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGAtata                                                                                                                                                   >  8:208500/1‑149 (MQ=255)   tGATCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGAtata                                                                                                                                                   >  6:35857/1‑149 (MQ=255)     aTCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATAc                                                                                                                                                  <  7:316485/148‑1 (MQ=255)     aTCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATAc                                                                                                                                                  <  7:292871/148‑1 (MQ=255)      tCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATAccc                                                                                                                                                <  5:54482/149‑1 (MQ=255)      tCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATAccc                                                                                                                                                <  5:60797/149‑1 (MQ=255)          cGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATg                                                                                                                                            <  8:246944/149‑1 (MQ=255)              ccGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGc                                                                                                                                        <  3:240284/149‑1 (MQ=255)                 aCGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCtt                                                                                                                                     <  3:384123/149‑1 (MQ=255)                   ggcggcGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTAGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTg                                                                                                                                    >  5:131632/1‑148 (MQ=255)                    gcggcgATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGaa                                                                                                                                  >  1:62173/1‑149 (MQ=255)                      ggcgATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAata                                                                                                                                <  6:264554/149‑1 (MQ=255)                       gcgATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAatat                                                                                                                               <  8:52879/149‑1 (MQ=255)                         gATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATg                                                                                                                              >  3:119863/1‑148 (MQ=255)                                 gAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCGGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAg                                                                                                                     >  6:413905/1‑149 (MQ=255)                                    ttttCAACGCCTATGAAAAACAGCTATGTCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGcc                                                                                                                  <  6:277229/149‑1 (MQ=255)                                       tCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCtt                                                                                                                >  6:308095/1‑148 (MQ=255)                                         aaCGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGa                                                                                                              <  6:56423/148‑1 (MQ=255)                                         aaCGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGa                                                                                                              <  7:208500/148‑1 (MQ=255)                                          aCGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATc                                                                                                            <  4:78725/149‑1 (MQ=255)                                           cGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATcc                                                                                                           <  7:226374/149‑1 (MQ=255)                                              cTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCa                                                                                                          >  4:46766/1‑147 (MQ=255)                                                   aaaaaCAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATg                                                                                                   <  7:255281/149‑1 (MQ=255)                                                    aaaaCAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATTTACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATgg                                                                                                  >  1:211830/1‑149 (MQ=255)                                                    aaaaCAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATTTACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATgg                                                                                                  >  1:205335/1‑149 (MQ=255)                                                     aaaCAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGc                                                                                                 <  2:212886/149‑1 (MQ=255)                                                      aaCATCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGcc                                                                                                <  8:161024/149‑1 (MQ=255)                                                         aGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGc                                                                                             <  8:114659/149‑1 (MQ=255)                                                                 cATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACcgcca                                                                                     >  4:16496/1‑149 (MQ=255)                                                                 cATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACcgcca                                                                                     >  4:268633/1‑149 (MQ=255)                                                                      cTGGATCGCCGACGATTGGCAGTCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACcgccacgcca                                                                                <  4:254865/149‑1 (MQ=255)                                                                      cTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACcgccacgcca                                                                                >  6:409389/1‑149 (MQ=255)                                                                       tGGATCGCCGACGATTGGCAGCCGTCCTTCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCgccacgccac                                                                               >  8:242657/1‑149 (MQ=255)                                                                        ggATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCACCCGCCACGCCAcc                                                                              >  2:249626/1‑149 (MQ=255)                                                                         gATCGCCGACGATTGGCAGCGGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCgccacgccac                                                                               <  3:268633/147‑1 (MQ=255)                                                                         gATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCAccc                                                                             <  5:35857/149‑1 (MQ=255)                                                                          aTCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCAcccc                                                                            <  2:21891/149‑1 (MQ=255)                                                                             gCCGACGATTGGCGGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGa                                                                         <  5:413905/149‑1 (MQ=255)                                                                              ccGACGATTTGCAGCCGTCCATCCGCCGCGTATCACATGTCAGATGGGGGCTGCCTGAGGAGCATGAAGATATACCATATGACGCCTTGAATATGGACCCCCTGACCTTGATCCAGAGGATGGCAGGCAACCGCCACGCACCCCCAGag                                                                        >  8:174983/1‑148 (MQ=255)                                                                              ccGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGAc                                                                        >  6:52051/1‑149 (MQ=255)                                                                                   gaTTGGCAGCCGTCCTTCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGAtt                                                                   <  5:365936/149‑1 (MQ=255)                                                                                   gaTTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCCCGCCACCCCAGACGGAtt                                                                   >  8:280827/1‑149 (MQ=255)                                                                                     ttGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATtgt                                                                 <  3:46766/149‑1 (MQ=255)                                                                                                 cATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCAcgcgcg                                                     >  3:140678/1‑149 (MQ=255)                                                                                                      gccgcGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTAc                                                <  6:250860/149‑1 (MQ=255)                                                                                                       ccgcGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACa                                               >  7:240753/1‑149 (MQ=255)                                                                                                        cgcgTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACAt                                              <  1:46331/149‑1 (MQ=255)                                                                                                           gTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATgcg                                           >  6:239158/1‑149 (MQ=255)                                                                                                              tCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCACCCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGcgcg                                         >  6:278324/1‑148 (MQ=255)                                                                                                                aCCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTa                                       <  5:52051/148‑1 (MQ=255)                                                                                                                     gCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCggg                                 <  2:62173/149‑1 (MQ=255)                                                                                                                            ggggCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTgcgc                           <  3:16496/148‑1 (MQ=255)                                                                                                                                     tGAGGAGCCTGAAGATATACCCTATGACGCATTGAATATGGACCCCAAGGCCTTGATCCAGACGATGGCCGGCAACCGACACTCCACCCCAGACCGATTGTGTGCCACGCGCGCCAACATGCGCGTAACCGGGCTGCGCCAGGAGGAAc                 >  2:101445/1‑149 (MQ=255)                                                                                                                                      gaggagCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAAcc                >  6:427832/1‑149 (MQ=255)                                                                                                                                         gagCCTGAAGATTTACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACcgcg             <  2:211830/149‑1 (MQ=255)                                                                                                                                         gagCCTGAAGATTTACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACcgcg             <  2:205335/149‑1 (MQ=255)                                                                                                                                           gCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGtt           <  4:140678/149‑1 (MQ=255)                                                                                                                                            ccTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGtt           >  7:154226/1‑148 (MQ=255)                                                                                                                                            ccTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTtg          >  3:282689/1‑149 (MQ=255)                                                                                                                                             cTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTtga         <  3:171133/149‑1 (MQ=255)                                                                                                                                             cTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTtga         <  5:308095/149‑1 (MQ=255)                                                                                                                                                aaGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTtgatga      >  1:346213/1‑149 (MQ=255)                                                                                                                                                  gATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATtgtg                                                                >  4:139807/1‑89 (MQ=255)                                                                                                                                                  gATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATtgtg                                                                <  3:139807/89‑1 (MQ=255)                                                                                                                                                   atatACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCAACCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTTGATGATGc   >  8:124734/1‑149 (MQ=255)                                                                                                                                                    tataCCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTTGATGATGct  <  2:346213/149‑1 (MQ=255)                                                                                                                                                    tataCCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTTGATGATGct  <  5:427832/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCTGATCGACGGGCCGACGGCGGCGATGTTCGGAGTTTTCAACGCCTATGAAAAACAGCTATGGCATGACTGGATCGCCGACGATTGGCAGCCGTCCATCCGCCGCGTATCACCTGGCAGTTGGGGGCTGCCTGAGGAGCCTGAAGATATACCCTATGACGCCTTGAATATGGACCCCAAGGCCTTGATCCAGAGGATGGCCGGCAACCGCCACGCCACCCCAGACGGATTGTGTGCCACGCGCGCCTACATGCGCGTAACCGGGCTGCGCCAGGAGGAACCGCGTTGATGATGCT  >  NC_002947/3686636‑3686931

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