Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,744,531	C→T	7.9%	P24S (CCG→TCG)	PP_2399 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,744,531	0	C	T	7.9%	180.9 / 4.4	63	P24S (CCG→TCG)	PP_2399	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (3/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

CCCTTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTACGTGATT  >  NC_002947/2744392‑2744672                                                                                                                                             |                                                                                                                                               cccTTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCg                                                                                                                                       <  1:314585/149‑1 (MQ=255)   cTTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGc                                                                                                                                     >  3:109189/1‑149 (MQ=255)    ttGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCg                                                                                                                                    >  3:143331/1‑149 (MQ=255)         cGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAg                                                                                                                               <  4:143331/149‑1 (MQ=255)         cGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAg                                                                                                                               >  3:263936/1‑149 (MQ=255)                ttGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGgc                                                                                                                        >  5:49404/1‑149 (MQ=255)                    ccGAGACGGGGCGTATGCGTGGCGTGCAATGCGCTTGCGGTGGGCAATCAGTTGAATGAGACCCTGCGTGCGCAGGAGTATGCCTTGGCCCGTCCCCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGCCGGGCCTCGACGccg                                                                                                                    <  8:50281/149‑2 (MQ=255)                       agaCGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGaa                                                                                                                 <  6:49404/149‑1 (MQ=255)                        gaCGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAg                                                                                                                <  5:230956/149‑1 (MQ=255)                        gaCGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAg                                                                                                                <  8:107648/149‑1 (MQ=255)                        gaCGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAg                                                                                                                <  7:265290/149‑1 (MQ=255)                                      gtgccgtgcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCg                                                                                                                                 <  8:117967/118‑1 (MQ=255)                                      gtgacgtgcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCg                                                                                                                                 >  7:117967/1‑118 (MQ=255)                                         ccgtgcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAgg                                                                                                               <  1:285791/133‑1 (MQ=255)                                         ccgtgcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAgttgtt                                                                                               <  3:291303/149‑1 (MQ=255)                                         cagtgcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTACGCGACACACACCGCC‑ATCCCCCGCCGCCAGGCGTCGAGGCCCCGCCCCTGaacg                                                                                                               >  2:285791/3‑131 (MQ=255)                                             gcAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTAt                                                                                           >  2:240144/1‑149 (MQ=255)                                              cAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATg                                                                                          <  2:240358/149‑1 (MQ=255)                                                aTGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGa                                                                                                     >  4:52721/1‑136 (MQ=255)                                                aTGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGa                                                                                                     <  3:52721/136‑1 (MQ=255)                                                   cgcTGCCGGAGGGCAAGCAGATGAATGACACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGgcgc                                                                                      <  4:241783/148‑1 (MQ=255)                                                        ccGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCg                                                                                 >  3:230262/1‑148 (MQ=255)                                                           gAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAgg                                                                             <  1:240144/149‑1 (MQ=255)                                                            aGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGc                                                                            >  5:141532/1‑149 (MQ=255)                                                              ggCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCtt                                                                          >  5:71005/1‑149 (MQ=255)                                                                 aaGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTgct                                                                       <  6:257935/149‑1 (MQ=255)                                                                     aGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTgctggct                                                                   <  4:230262/149‑1 (MQ=255)                                                                         gaatgaGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAgttgtt                                                                                               >  1:309041/1‑117 (MQ=255)                                                                         gaatgaGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAgttgtt                                                                                               <  2:309041/117‑1 (MQ=255)                                                                            tgaGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCt                                                            <  4:111335/149‑1 (MQ=255)                                                                               gaCCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTgctgg                                                                     <  4:296576/137‑1 (MQ=255)                                                                                aCCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTccc                                                        >  6:133797/1‑149 (MQ=255)                                                                                  ccTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCgg                                                      >  5:41587/1‑149 (MQ=255)                                                                                    tgcgttcgCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGACGAGGCCCCGCGCCAGAAGGGTTACCGTGAGGAGTTCTATGCCGCCATCGAAGGCTTGCAGGCTGGCAGATTCCGGGt                                                     >  8:289930/1‑148 (MQ=255)                                                                                     gcgtgcgCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGa                                                   >  1:290709/1‑149 (MQ=255)                                                                                       gtgcgCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCACCGCCGCACGGAGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCAGGCTGGCAGCTGCCCGGTGATg                                                 >  1:292694/1‑149 (MQ=255)                                                                                            cagcagTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCa                                            <  4:109189/149‑1 (MQ=255)                                                                                                    tGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTGTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTgg                                    <  5:160299/149‑1 (MQ=255)                                                                                                      ccTTGGCCCGTCACCTGCGCGACCCACAGCGCCAATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTggg                                   >  1:205352/1‑149 (MQ=255)                                                                                                      ccTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATc                                                                                  <  4:130522/101‑1 (MQ=255)                                                                                                      ccTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTTATGCTCCAGGCGTTTGGc                                  >  4:72611/1‑149 (MQ=255)                                                                                                             ccGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGa                                >  7:54121/1‑144 (MQ=255)                                                                                                             ccGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGa                                <  8:54121/144‑1 (MQ=255)                                                                                                                tCACGTTCGCGACCCTCAGCGAG‑ATGCCCCGCCGCCCGGCGTCGAGGGCCTGCGCCTGAATGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTTCTGGCTGGCAGCTTCCCGGTTATTCGCCAGGCGTTGGGCGAGCAGCG‑Ct                       <  2:292694/149‑1 (MQ=255)                                                                                                                  aCCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATg                                                 <  6:203770/122‑1 (MQ=255)                                                                                                                  aCCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATg                                                 >  5:203770/1‑122 (MQ=255)                                                                                                                    cTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCACCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCa                   >  1:113417/1‑149 (MQ=255)                                                                                                                      gcgcAACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGc                    >  1:199920/1‑146 (MQ=255)                                                                                                                      gcgcAACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGc                    <  2:199920/146‑1 (MQ=255)                                                                                                                         cgACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATgc                                                <  3:2335/116‑1 (MQ=255)                                                                                                                         cgACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATgc                                                >  4:2335/1‑116 (MQ=255)                                                                                                                              cacaGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTa         <  6:41587/149‑1 (MQ=255)                                                                                                                              cacaGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCA‑GCCTTTGGCCAGCAACGCCTGGCTGGCCCGAgaa         >  2:221336/1‑147 (MQ=255)                                                                                                                               acaGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTAc        >  5:57072/1‑149 (MQ=255)                                                                                                                                  gcgcC‑GGGCCTCGACGCCGGGCGGCGGGGCATGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCAGTGATGCGCCAGGCGTTg                                     <  7:218689/118‑1 (MQ=25)                                                                                                                                  gcgcC‑GGGCCTCGACGCCGGGCGGCGGGGCATGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCAGTGATGCGCCAGGCGTTg                                     >  8:218689/1‑118 (MQ=25)                                                                                                                                     gccgggCCTCGACGCCGGGCGGCGGGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGag                               <  6:265837/117‑1 (MQ=255)                                                                                                                                     gccgggCCTCGACGCCGGGCGGCGGGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGag                               >  5:265837/6‑122 (MQ=255)                                                                                                                                     gccgggCCTCGACGCCGGGCGGCGGGGCATGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTACGTGAt   >  8:164799/6‑149 (MQ=38)                                                                                                                                      c‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGa            <  4:178477/138‑1 (MQ=255)                                                                                                                                      c‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGa            >  3:178477/1‑138 (MQ=255)                                                                                                                                      c‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTACGTGAtt  <  7:36854/148‑1 (MQ=255)                                                                                                                                      c‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTACGTGAtt  >  8:36854/1‑148 (MQ=255)                                                                                                                                             ccgccgCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGc                    <  6:173390/124‑1 (MQ=255)                                                                                                                                             |                                                                                                                                               CCCTTGGCCGAACTGTTGGCCGAGACGGCGCGTATGCGTGCCGTGCAATGCGCTGCCGGAGGGCAAGCAGATGAATGAGACCCTGCGTGCGCAGCAGTATGCCTTGGCCCGTCACCTGCGCGACCCACAGCGCC‑ATGCCCCGCCGCCCGGCGTCGAGGCCCGGCGCCTGAAGGTTTACCGTGAGTTGTTCTATGGCGCCATCGAAGGCTTGCTGGCTGGCAGCTTCCCGGTGATGCGCCAGGCGTTGGGCGAGCAGCG‑CTGGCATGCCCGAGTACGTGATT  >  NC_002947/2744392‑2744672

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