Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,240,571	C→G	11.8%	L45V (CTG→GTG)	PP_4618 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,240,571	0	C	G	11.8%	206.4 / 18.2	76	L45V (CTG→GTG)	PP_4618	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base G (3/6);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-02

ATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTGCA  >  NC_002947/5240432‑5240719                                                                                                                                               |                                                                                                                                                     aTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                             <  7:244693/149‑2 (MQ=255) aTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                             <  4:63433/149‑2 (MQ=255)    ccccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGTGCCTGTCGACCTGCTTGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCgag                                                                                                                                           <  4:219018/148‑1 (MQ=255)    acccATGCACCCACGGCAATGCTTCCTGGCCTTCCTGGCAAGCGATCCTGTCGTACTGCTTGTATCCGCACTTTTGCTTGCTCCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCCGCGCTGCAGGTCGTGCAgcgca                                                                                                                                             <  1:272049/148‑2 (MQ=255)     cccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCCGCGCTGCAGGTCGTGCAgcgcagcaa                                                                                                                                         <  5:93665/149‑6 (MQ=255)      ccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                         >  5:275256/1‑148 (MQ=255)       cATGAAGCCACGTCAAGCCTGCCTGGCCTCCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CGGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                         >  8:238614/1‑147 (MQ=255)         tGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                     <  3:61465/149‑1 (MQ=255)                aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAcc                                                                                                                              >  6:187202/1‑149 (MQ=255)                   tCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCCGCAGCGCGAGATAAGCGCCAACCTGc                                                                                                                           <  4:254817/149‑1 (MQ=255)                       gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGcctct                                                                                                                       >  6:234592/1‑146 (MQ=255)                       gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                       <  8:86540/149‑1 (MQ=255)                       gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                       >  6:6400/1‑149 (MQ=255)                       gcctgcatggcctgcatgGAACCCGAGCCTGTCGACCACCTCGAACCCGCAATGAGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGCCCCCGCTCCCAGCCATGCTGAGCTGCACCACCTGCAGCCCGAGAAAACCGCCAACCTGCCCAt                                                                                                                       >  6:332195/1‑149 (MQ=255)                            cctggcctacctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                                  <  6:226904/149‑1 (MQ=255)                                gcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTTGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACt                                                                                                              <  5:187202/149‑1 (MQ=255)                                  ctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTAt                                                                                                            >  2:260765/1‑149 (MQ=255)                                  ctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTAt                                                                                                            >  4:245487/1‑149 (MQ=255)                                     cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                          >  5:119556/1‑148 (MQ=255)                                     cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                          >  3:95786/1‑148 (MQ=255)                                     cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTAGAGCCCGCTGCCAGCCGTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                          <  5:234592/148‑1 (MQ=255)                                      ctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCg                                                                                                         >  1:47126/1‑148 (MQ=255)                                       tggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAACCGCCAACATGCCGATGCTTCAACTATGCGaa                                                                                                       >  7:2724/1‑149 (MQ=255)                                        ggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                             >  3:109713/1‑109 (MQ=255)                                        ggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                             <  4:109713/110‑2 (MQ=255)                                          aaCGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTg                                                                                                    <  4:95786/149‑1 (MQ=255)                                           aCGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTgg                                                                                                   >  5:140763/1‑149 (MQ=255)                                               gAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCACCCTGCCGATGCTTCCACTATGCGAACTGGcc                                                                                                 >  5:112101/1‑147 (MQ=255)                                                 gCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGc                                                                                             <  6:119556/149‑1 (MQ=255)                                                    tGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCAc                                                                                          >  6:93343/1‑149 (MQ=255)                                                        gACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgct                                                                                      >  7:30865/1‑149 (MQ=255)                                                           ctgctgGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcg                                                                                   >  2:214646/1‑149 (MQ=255)                                                              ctgGAAGTCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTCCCGATGCTTCCACTATGCGAACTGGCACAGCCACTGCTccgc                                                                                  >  8:178968/1‑147 (MQ=255)                                                               tgGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAg                                                                               >  7:195784/1‑149 (MQ=255)                                                                 gAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCt                                                                             >  6:187416/1‑149 (MQ=255)                                                                    gCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGaa                                                                          <  5:73229/149‑1 (MQ=255)                                                                     ccGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAAt                                                                         >  7:153093/1‑149 (MQ=255)                                                                       gCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGc                                                                                             <  7:280960/127‑1 (MQ=255)                                                                       gCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGc                                                                                             >  8:280960/1‑127 (MQ=255)                                                                       gCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATgc                                                                       >  2:277518/1‑149 (MQ=255)                                                                           tgtgGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTg                                                                   >  2:329972/1‑149 (MQ=255)                                                                            gtgGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                                  >  5:109655/1‑149 (MQ=255)                                                                             tgGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                                 <  7:178968/149‑1 (MQ=255)                                                                             tgGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTg                                                                   <  1:329972/147‑1 (MQ=255)                                                                                aTTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCtt                                                              <  6:275256/149‑1 (MQ=255)                                                                                  tg‑‑‑ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                             >  1:57272/1‑148 (MQ=255)                                                                                   g‑‑‑ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggttctgat                                                           >  8:124826/1‑142 (MQ=255)                                                                                       ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                          <  3:245487/149‑1 (MQ=255)                                                                                          ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagcag                                                       <  3:325617/149‑1 (MQ=255)                                                                                          ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                        <  8:153093/148‑1 (MQ=255)                                                                                          ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                        <  1:214646/148‑1 (MQ=255)                                                                                          ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                        <  5:6400/148‑1 (MQ=255)                                                                                           cGAGCACGACCCTGGTGTCGATCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAgg                                                      <  6:112101/149‑1 (MQ=255)                                                                                                aCGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGaa                                                 <  1:83454/149‑1 (MQ=255)                                                                                                     ccTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                            <  2:47126/149‑1 (MQ=255)                                                                                                         ttgtCGAGCCCGCTGCCAGCCTTGGTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAACTGTCGATCCTTCCACTATGCGAACTGGCCCATCGACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAcc                                          <  5:332195/146‑1 (MQ=255)                                                                                                            gcGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCtt                                      <  6:279306/147‑1 (MQ=255)                                                                                                              gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACtgctgc                                                                                    <  2:270468/100‑1 (MQ=255)                                                                                                              gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACtgctgc                                                                                    >  1:270468/1‑100 (MQ=255)                                                                                                              gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                            <  7:151860/140‑1 (MQ=255)                                                                                                              gAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                            >  8:151860/1‑140 (MQ=255)                                                                                                                 cccGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGccc                                >  1:326441/1‑149 (MQ=255)                                                                                                                          aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATgcggc                       >  5:52121/1‑149 (MQ=255)                                                                                                                                gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgatg                   <  1:277518/147‑1 (MQ=255)                                                                                                                                  tgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATgg               <  7:238614/149‑1 (MQ=255)                                                                                                                                     gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg                >  1:327498/1‑145 (MQ=255)                                                                                                                                     gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg                <  2:327498/145‑1 (MQ=255)                                                                                                                                     gctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGAc             >  4:250564/1‑148 (MQ=255)                                                                                                                                        gcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga                     <  1:172363/137‑1 (MQ=255)                                                                                                                                        gcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCga                     >  2:172363/1‑137 (MQ=255)                                                                                                                                         cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgc                                    >  4:47310/1‑121 (MQ=255)                                                                                                                                         cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgc                                    <  3:47310/121‑1 (MQ=255)                                                                                                                                         cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg         <  6:140763/148‑1 (MQ=255)                                                                                                                                         cACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAg          <  5:93343/147‑1 (MQ=255)                                                                                                                                            gACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACt     >  1:213032/1‑149 (MQ=255)                                                                                                                                             aCCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTg    >  2:221015/1‑149 (MQ=255)                                                                                                                                               cTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTGCa  <  8:30865/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                     ATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTG‑‑‑CTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTACTGCA  >  NC_002947/5240432‑5240719

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