Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,591,021	G→A	6.4%	A73T (GCC→ACC)	PP_3169 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,591,021	0	G	A	6.4%	184.8 / 2.1	62	A73T (GCC→ACC)	PP_3169	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

ACGCTCGCTGCTGTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTCTGGCCTTTGCACACCACCCC  >  NC_002947/3590876‑3591169                                                                                                                                                   |                                                                                                                                                      aCGCTCGCTGCTGTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCct                                                                                                                                                    <  2:80693/149‑1 (MQ=255)   gctcgctGCTGTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCctgg                                                                                                                                                  <  5:147713/149‑1 (MQ=255)       gctgctGTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCctggctgg                                                                                                                                              <  8:234066/149‑1 (MQ=255)        ctgctgTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGt                                                                                                                                             <  7:346000/149‑1 (MQ=255)           ctgTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCaa                                                                                                                                          >  6:271394/1‑149 (MQ=255)             gTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACa                                                                                                                                        >  2:363374/1‑149 (MQ=255)              tATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAg                                                                                                                                       <  7:26268/149‑1 (MQ=255)               aTGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGAccagccag                                                                                                                                               <  4:177026/140‑7 (MQ=255)               aTGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGAccagccag                                                                                                                                               >  3:177026/1‑134 (MQ=255)                tGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGgc                                                                                                                                     <  6:359119/149‑1 (MQ=255)                   tGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCt                                                                                                                                  >  4:290281/1‑149 (MQ=255)                   tGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCt                                                                                                                                  >  2:48720/1‑149 (MQ=255)                     cagcaCTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTgg                                                                                                                                >  3:99777/1‑149 (MQ=255)                             ggCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCctgg                                                                                                                                                  <  3:124630/123‑1 (MQ=255)                             ggCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGAGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCctgg                                                                                                                                                  >  4:124630/1‑123 (MQ=255)                              gCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAg                                                                                                                       <  5:84977/149‑1 (MQ=255)                               cACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAgg                                                                                                                      <  5:266745/149‑1 (MQ=255)                                aCGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAggg                                                                                                                     <  5:191373/149‑1 (MQ=255)                                aCGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAggg                                                                                                                     <  8:273959/149‑1 (MQ=255)                                aCGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAggg                                                                                                                     <  5:26708/149‑1 (MQ=255)                                aCGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAggg                                                                                                                     <  4:99777/149‑1 (MQ=255)                                  gCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGt                                                                                                                   <  1:48720/149‑1 (MQ=255)                                  gCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGt                                                                                                                   <  3:290281/149‑1 (MQ=255)                                           ccTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCACCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGTATACCCCTTGc                                                                                                          <  7:54167/149‑1 (MQ=255)                                            cTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGGTATTCTCCCACTCGTTGTACGTGCTAACCGGCCTCGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCACGGCAAGGGGTATACCTTGc                                                                                                          >  1:285911/1‑149 (MQ=255)                                            cTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGACCAGCCAGGCCAACAg                                                                                                                                       <  3:108538/119‑1 (MQ=255)                                            cTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGACCAGCCAGGCCAACAg                                                                                                                                       >  4:108538/1‑119 (MQ=255)                                             tGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGTCCCGGCAAGGGGTATACCTTGcc                                                                                                         >  7:117865/1‑148 (MQ=255)                                               gATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAg                                                                                                      >  6:183375/1‑149 (MQ=255)                                                aTGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAg‑g                                                                                                    <  5:271394/149‑1 (MQ=255)                                                   tgATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGccc                                                                                                                             <  2:15550/122‑1 (MQ=255)                                                   tgATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGccc                                                                                                                             >  1:15550/1‑122 (MQ=255)                                                     aTCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCtttt                                                                                               >  3:12627/1‑149 (MQ=255)                                                        aTTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCtttttt                                                                                             <  5:129175/148‑1 (MQ=255)                                                            gCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCg                                                                                        >  5:300036/1‑149 (MQ=255)                                                              tATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGt                                                                                       >  5:340250/1‑148 (MQ=255)                                                                tGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTc                                                                                     <  4:12627/148‑1 (MQ=255)                                                                 gCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTct                                                                                    >  4:229733/1‑148 (MQ=255)                                                                 gCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTct                                                                                    <  3:229733/148‑1 (MQ=255)                                                                   cGACCCGGTGTCGCAGATGACTTATCATGGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTg                                                                                           >  1:207315/1‑139 (MQ=255)                                                                   cGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTg                                                                                           <  2:207315/139‑1 (MQ=255)                                                                       ccGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTctggctgg                                                                              >  2:80208/1‑148 (MQ=255)                                                                        cGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTTGCTGGTa                                                                            >  6:328094/1‑149 (MQ=255)                                                                             tCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGc                                                                                                                          <  4:243403/99‑1 (MQ=255)                                                                             tCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGc                                                                                                                          >  3:243403/1‑99 (MQ=255)                                                                                     gACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCCGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTGCTGGCTGGTATTGGTAACCAATc                                                               >  6:208501/1‑149 (MQ=255)                                                                                       cTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCg                                                             <  5:183375/149‑1 (MQ=255)                                                                                        ttATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTctggct                                                                                >  7:189308/1‑129 (MQ=255)                                                                                        ttATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTctggct                                                                                <  8:189308/129‑1 (MQ=255)                                                                                           tcatcGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGc                                                         >  5:349783/1‑149 (MQ=255)                                                                                              tcGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCg                                                      <  1:363374/149‑1 (MQ=255)                                                                                                         tcCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCAc                                            >  2:309459/1‑149 (MQ=255)                                                                                                           ccACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACgg                                          >  1:180548/1‑149 (MQ=255)                                                                                                           ccACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACgg                                          >  1:174823/1‑149 (MQ=255)                                                                                                             aCTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGAc                                                    >  7:295104/1‑137 (MQ=255)                                                                                                             aCTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGAc                                                    <  8:295104/137‑1 (MQ=255)                                                                                                                            tCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAgttg                         >  4:307758/1‑149 (MQ=255)                                                                                                                              aCCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAgttgtt                       >  6:204639/1‑149 (MQ=255)                                                                                                                                  gcctggcctTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTCTgg                   <  6:349783/149‑1 (MQ=255)                                                                                                                                                ttGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAGAGCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTCTGGCCTTTGCAcacca      >  1:96084/1‑149 (MQ=255)                                                                                                                                                ttGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTc                      <  4:24962/132‑1 (MQ=255)                                                                                                                                                ttGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTc                      >  3:24962/1‑132 (MQ=255)                                                                                                                                                   gCCTGGCTGGTCAACAGTCGCTGGCCCGGCAAGTGGTATACCTTGCCCAG‑GCt                                                                                                  <  3:94908/53‑1 (MQ=255)                                                                                                                                                   gCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGACTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTCTGGCCTTTGCACACCAcccc  >  2:128471/1‑149 (MQ=255)                                                                                                                                                   gCCTGGCTGGTCAACAGGCGCTCGCCCGGCAAGGGGTATACCTTGCCCAG‑GCt                                                                                                  >  4:94908/1‑53 (MQ=255)                                                                                                                                                   |                                                                                                                                                      ACGCTCGCTGCTGTATGGTGCAGCACTGGGCACGCTGCCGGACCTGGATGTGATCATTCGCTATGCCGACCCGGTGTCGCAGATGACTTATCATCGCGG‑CTTCTCCCACTCGTTGTTCGTGCTCACCGGCCTGGCCTTGCTGTTGGCCTGGCTGGTCAACAGGCGCTGGCCCGGCAAGGGGTATACCTTGCCCAG‑GCTTTTTTTGGCGTTCTGGCTGGTATTGGTAACCCATCCGGTGCTCGACGCCTTCACGGTTTACGGCACGCAGTTGTTCTGGCCTTTGCACACCACCCC  >  NC_002947/3590876‑3591169

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