Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,375	T→G	5.5%	A496A (GCA→GCC)	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,375	0	T	G	5.5%	234.3 / 4.0	73	A496A (GCA→GCC)	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base G (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

TTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCAC  >  NC_002947/5798234‑5798518                                                                                                                                              |                                                                                                                                                ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                          >  3:302880/1‑149 (MQ=255) ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                          >  4:47125/1‑149 (MQ=255) ttAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCATGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                                                                                                                          >  5:170437/1‑149 (MQ=255)        cTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                   <  4:172200/149‑1 (MQ=255)         tGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                                                                                                                   >  7:172446/1‑148 (MQ=255)           ggcggcGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCGACGCAGCGGTTgcgg                                                                                                                                  <  8:217143/147‑4 (MQ=255)            gcggcgGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGaaa                                                                                                                               >  2:196802/1‑149 (MQ=255)               gcggGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTg                                                                                                                            >  4:176108/1‑149 (MQ=255)                         gccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCa                                                                                                                  >  8:118853/1‑149 (MQ=255)                          ccgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                                 >  1:140856/1‑149 (MQ=255)                           cgccCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                                                                                                                 >  1:26630/1‑148 (MQ=255)                                ggAAGGTATCAGCGGGCCGTGGCCAGCTCTTGGTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTa                                                                                                           <  6:323988/149‑1 (MQ=255)                                       aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTgg                                                                                                                                    <  6:127234/117‑1 (MQ=255)                                       aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTgg                                                                                                                                    >  5:127234/1‑117 (MQ=255)                                        tCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCt                                                                                                   >  7:100244/1‑149 (MQ=255)                                          aGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg                                                                                                 <  2:26630/149‑1 (MQ=255)                                           gCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGc                                                                                                <  6:107923/149‑1 (MQ=255)                                                ccGTGGCCAGCTCTTGTTGGACGCTGCCTACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGc                                                                                           <  3:176108/149‑1 (MQ=255)                                                ccGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTg                                                                                            >  2:264269/1‑148 (MQ=255)                                                      ccAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGttt                                                                                     >  7:232726/1‑149 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCtt                                                                                <  8:100244/147‑1 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                               >  6:29713/1‑148 (MQ=255)                                                             ttgttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                               >  5:74494/1‑148 (MQ=255)                                                               gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                             <  7:169144/148‑1 (MQ=255)                                                               gttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGaa                                                                             <  4:302880/148‑1 (MQ=255)                                                                ttgGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAg                                                                            >  4:173795/1‑148 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                         <  1:196802/148‑1 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                         <  1:264269/148‑1 (MQ=255)                                                                   gACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAg                                                                         <  7:118853/148‑1 (MQ=255)                                                                      gCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGAt                                                                     >  6:62020/1‑149 (MQ=255)                                                                          ccGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTc                                                                 >  4:250377/1‑149 (MQ=255)                                                                           cGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTct                                                                <  4:174688/149‑2 (MQ=255)                                                                           cGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCg                                                                <  4:174689/149‑1 (MQ=255)                                                                                    gTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTgg                                                       >  8:278964/1‑149 (MQ=255)                                                                                       gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                               >  5:61785/1‑122 (MQ=255)                                                                                       gCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTg                                                                               <  6:61785/122‑1 (MQ=255)                                                                                         cAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                  >  2:104641/1‑149 (MQ=255)                                                                                         cAAGGTTTACCAGCCTGGACATTAGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                  >  2:221973/1‑149 (MQ=255)                                                                                             gTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGcc                                                              >  7:61970/1‑133 (MQ=255)                                                                                             gTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGcc                                                              <  8:61970/133‑1 (MQ=255)                                                                                                 aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTaa                                           <  3:47125/148‑1 (MQ=255)                                                                                                 aCCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTaa                                           <  4:118015/148‑1 (MQ=255)                                                                                                   cAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAg                                          <  6:40998/147‑1 (MQ=255)                                                                                                   cAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAg                                          >  5:40998/1‑147 (MQ=255)                                                                                                     gCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTg                                      <  6:74494/149‑1 (MQ=255)                                                                                                     gCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTg                                      >  8:291076/1‑149 (MQ=255)                                                                                                      ccTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTgg                                     <  7:31967/149‑1 (MQ=255)                                                                                                       cTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                      <  7:246000/131‑1 (MQ=255)                                                                                                       cTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTggg                                                      >  8:246000/1‑131 (MQ=255)                                                                                                        tGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACTCTGGGCCTTGACGTAAGCCTGGGc                                   >  8:314457/1‑149 (MQ=255)                                                                                                                cGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCt                           >  4:233757/1‑149 (MQ=255)                                                                                                                 gCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTc                          <  5:62020/149‑1 (MQ=255)                                                                                                                      cAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGc                     >  6:80401/1‑149 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGCGCt                    >  8:301537/1‑149 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                    <  6:170437/149‑1 (MQ=255)                                                                                                                       aaGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCt                    <  8:172446/149‑1 (MQ=255)                                                                                                                        aGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTg                   <  3:233757/149‑1 (MQ=255)                                                                                                                          ttCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTgg                  >  8:154249/1‑148 (MQ=255)                                                                                                                             ttGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCg                         >  6:52493/1‑138 (MQ=255)                                                                                                                             ttGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCg                         <  5:52493/138‑1 (MQ=255)                                                                                                                               gATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTTGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTAcc              <  4:232747/147‑1 (MQ=255)                                                                                                                                aTTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCaa            <  5:80401/148‑1 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCa             >  2:82276/1‑146 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCa             <  1:82276/146‑1 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg           >  4:205453/1‑148 (MQ=255)                                                                                                                                 ttGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAg           >  3:93282/1‑148 (MQ=255)                                                                                                                                   gCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGa         <  1:104641/148‑1 (MQ=255)                                                                                                                                    ccGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGAt        <  3:250377/148‑1 (MQ=255)                                                                                                                                     gccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCt      >  7:290051/3‑149 (MQ=255)                                                                                                                                      ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                  <  8:75426/103‑1 (MQ=255)                                                                                                                                      ccAACCGCGACGCAGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCtt                                                  >  7:75426/2‑104 (MQ=255)                                                                                                                                        aaCCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGca   >  1:135577/1‑149 (MQ=255)                                                                                                                                         aCCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGcac  <  7:278964/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                TTAGCTACTGGGCGGCGGGGTGATGCCGCCCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTGCCCTGCAGGTTTTTCTTGAAGTAGCGATAGTCGCCACCCTGGGCCTTGACGTAAGCCTGGGCCGACGCCTCGGAGCTGGTACCAAGGATCTGCAC  >  NC_002947/5798234‑5798518

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