Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,184,280	T→A	6.4%	W20R (TGG→AGG)	PP_3681 →	putative Helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,184,280	0	T	A	6.4%	188.4 / 2.6	62	W20R (TGG→AGG)	PP_3681	putative Helicase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/2);  total (32/30)
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.46e-01

ACTGACACGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGGTGTCGTTCAA  >  NC_002947/4184137‑4184428                                                                                                                                                 |                                                                                                                                                          aCTGACACGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGa                                                                                                                                                        <  5:143028/147‑1 (MQ=255)     acacGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCgtt                                                                                                                                                   <  4:175910/148‑1 (MQ=255)     acacGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCgtt                                                                                                                                                   <  3:255266/148‑1 (MQ=255)     acacGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCgtt                                                                                                                                                   <  4:188423/148‑1 (MQ=255)          gAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTc‑gc                                                                                                                                            >  8:284687/1‑149 (MQ=255)           aGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTc‑gc                                                                                                                                            >  5:90856/1‑148 (MQ=255)              tgctgGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATgg                                                                                                                                                         >  3:330085/1‑133 (MQ=255)              tgctgGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATgg                                                                                                                                                         <  4:330085/133‑1 (MQ=255)                ctgGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCtt                                                                                                                                      >  4:169064/1‑149 (MQ=255)                   gCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTgat                                                                                                                                   <  7:284687/149‑1 (MQ=255)                          tACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCaa                                                                                                                            >  3:146049/1‑149 (MQ=255)                          tACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCAAGGAGCGTTGTTC‑GCGAGCTTGATCGATCaa                                                                                                                            >  2:3533/1‑149 (MQ=255)                             caccaTCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATc                                                                                                                          <  5:10070/148‑1 (MQ=255)                             caccaTCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCt                                                                                                                         >  8:159794/1‑149 (MQ=255)                             caccaTCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCCATCAATCt                                                                                                                         >  6:171986/1‑149 (MQ=255)                             caccaTCCCTCAAACGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATc                                                                                                                          >  6:10070/1‑148 (MQ=255)                               ccaTCCCTCAATCGTCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTTATTTCTCATGATCA‑ATTGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTg                                                                                                                       <  3:136381/149‑1 (MQ=255)                               ccaTCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTg                                                                                                                       <  5:121078/149‑1 (MQ=255)                                caTCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGc                                                                                                                      >  2:42133/1‑149 (MQ=255)                                        aaTCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACcgcg                                                                                                              >  8:304061/1‑149 (MQ=255)                                          tCGGCTAGCAGACTGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCggg                                                                                                            >  4:79828/1‑149 (MQ=255)                                               gatcAGCCCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAg                                                                                                       <  4:239871/146‑1 (MQ=255)                                                 gCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCAGTTGGGCTCCCAGCGTCGGATTACGCAGGGAGAAAAATGCATGGACCGTTTTTG‑GCCTCCTTGAGCGAAAAAAGTTGCTCACCGGGGGGCCAgg                                                                                                      >  1:277662/1‑149 (MQ=255)                                                 gCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTgatcgat                                                                                                                               <  4:226212/123‑1 (MQ=255)                                                 gCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTgatcgat                                                                                                                               >  3:226212/1‑123 (MQ=255)                                                 gCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCAAGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCATGt                                                                                                     <  1:3533/149‑1 (MQ=255)                                                  cAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCAAGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTg                                                                                                    >  6:98585/1‑149 (MQ=255)                                                     aCCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGcc                                                                                                         <  7:145214/141‑1 (MQ=255)                                                     aCCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGcc                                                                                                         >  8:145214/1‑141 (MQ=255)                                                       cGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGc                                                                                            <  1:42133/149‑1 (MQ=255)                                                           tCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGa                                                                                        <  4:146049/149‑1 (MQ=255)                                                             aaaGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGAcc                                                                                      <  3:169064/149‑1 (MQ=255)                                                                gCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATTATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAg                                                                                   <  7:159794/149‑1 (MQ=255)                                                                gCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAg                                                                                   <  2:277662/149‑1 (MQ=255)                                                                gCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAg                                                                                   >  2:234330/1‑149 (MQ=255)                                                                          aGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCATCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGccc                                                                         <  6:90856/149‑1 (MQ=255)                                                                          aGGAATCTCCCGAGTGATTTCTCATGAGCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGATC‑GCGAGCTTGATCGATCAATCTTGCGCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGcc                                                                          <  5:359879/148‑1 (MQ=255)                                                                            gAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGTGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCaa                                                                       >  2:45332/1‑149 (MQ=255)                                                                             aaTCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAAc                                                                      <  5:96398/149‑1 (MQ=255)                                                                                 tcCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGc                                                                  >  6:173844/1‑149 (MQ=255)                                                                                   ccGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGTGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGcc                                                                 <  1:45332/148‑1 (MQ=255)                                                                                   ccGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGcc                                                                 >  5:365219/1‑148 (MQ=255)                                                                                    cGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAg                                                                                   <  4:99798/129‑1 (MQ=255)                                                                                    cGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAg                                                                                   >  3:99798/1‑129 (MQ=255)                                                                                    cGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCAAGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGa                                                               <  5:98585/149‑1 (MQ=255)                                                                                                       a‑aTGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑aa                                            <  6:365219/148‑1 (MQ=255)                                                                                                         aTGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGc                                                                                                                      >  2:231904/1‑77 (MQ=255)                                                                                                         aTGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGc                                                                                                                      <  1:231904/77‑1 (MQ=255)                                                                                                         aTGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCCCCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGCGGTACTGGCCCAACGTGCCGACTACCTGCTGCGCAccca                                             >  5:340361/1‑145 (MQ=255)                                                                                                                 ccTCGCTCGGAGTACGCCTCGAGAAAATTGCATGGAGCGTTGTAC‑CCGAGCATGAACCAACAATCAAGCTCACCGCGGGGCCAGGTCCC‑‑‑CGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGccc                                   >  1:21991/1‑147 (MQ=255)                                                                                                                   tcgctcgGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGaa                                                 <  7:83386/133‑1 (MQ=255)                                                                                                                   tcgctcgGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGaa                                                 >  8:83386/1‑133 (MQ=255)                                                                                                                     gctcgGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAAcc                                               <  7:174929/133‑1 (MQ=255)                                                                                                                     gctcgGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAAcc                                               >  8:174929/1‑133 (MQ=255)                                                                                                                            gTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAgc                        <  3:79828/148‑1 (MQ=255)                                                                                                                                  cTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCGTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCt                 >  6:169142/1‑149 (MQ=255)                                                                                                                                    ggAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTg               >  2:187920/1‑149 (MQ=255)                                                                                                                                     gagaAAAATGCATGGAGCGTTGTTCGGCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTg               >  1:5178/1‑149 (MQ=255)                                                                                                                                     gagaAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCCTTAGGCAAGACCGAGGTACTGGCCAAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCt                 >  1:111229/1‑149 (MQ=255)                                                                                                                                     gagaAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGc              >  7:101328/1‑149 (MQ=255)                                                                                                                                      agaAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCg             >  8:191479/1‑149 (MQ=255)                                                                                                                                         aaaaTGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACCAAACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGgt           <  1:234330/149‑1 (MQ=255)                                                                                                                                                 tGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGGTGTCGTTCaa  <  5:173844/149‑1 (MQ=255)                                                                                                                                                 tGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GACAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGgtg          >  2:110264/1‑141 (MQ=255)                                                                                                                                                 tGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GACAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGgtg          <  1:110264/141‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                          ACTGACACGGAGCTGCTGGCCTCACTACCACCATCCCTCAATCGGCTAGCAGACCGGGTCAAAGCCAAGCTTGAGGAATCTCCCGAGTGATTTCTCATGATCA‑ATGGGCTCCCTCGCTCGGAGTACGCCTGGAGAAAAATGCATGGAGCGTTGTTC‑GCGAGCTTGATCGATCAATCTTGCTCACCGCGGGGCCAGGTGCC‑‑‑GGCAAGACCGAGGTACTGGCCCAACGTGCCGACTACCTTCTGCGAACC‑AACAGCTGCCGCTTCCCCAAGCGCATCCTTGCGGTGTCGTTCAA  >  NC_002947/4184137‑4184428

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