Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,622,907	C→G	5.7%	G150G (GGG→GGC)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,622,907	0	C	G	5.7%	220.0 / 4.2	70	G150G (GGG→GGC)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base G (2/2);  total (35/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 = 4.66e-01

GGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGA  >  NC_002947/1622761‑1623052                                                                                                                                                   |                                                                                                                                                  ggAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGccc                                                                                                                                                 >  2:194366/1‑149 (MQ=255)  gAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGTTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGACATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGccca                                                                                                                                                >  4:57260/1‑149 (MQ=255)         ccAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCCTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAgg                                                                                                                                         >  8:16538/1‑149 (MQ=255)           aGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTg                                                                                                                                       <  2:231841/149‑1 (MQ=255)               cTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCt                                                                                                                                   >  6:362311/1‑149 (MQ=255)                  gAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGt                                                                                                                                >  8:53869/1‑149 (MQ=255)                  gAGGATCATGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTg                                                                                                                                  <  5:331663/147‑1 (MQ=255)                           gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaa                                                                                                                             <  5:358202/143‑3 (MQ=255)                           gTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGGCCACCTTGGGTGGGCTGGTcaa                                                                                                                             >  6:358202/1‑141 (MQ=255)                            tCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTga                                                                                                                       <  1:194366/148‑1 (MQ=255)                              ccccAGCACCGGCTGGAGTTGTTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGACATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                    <  3:57260/149‑1 (MQ=255)                              ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                           >  7:288854/1‑142 (MQ=255)                              ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                           <  8:288854/142‑1 (MQ=255)                              ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                     <  6:214708/148‑1 (MQ=255)                              ccccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                    <  2:187278/149‑1 (MQ=255)                               cccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAgg                                                                                                                   <  8:49726/149‑1 (MQ=255)                               cccAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAgg                                                                                                                   >  7:49726/1‑149 (MQ=255)                                   gCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGgtgt                                                                                                                >  6:104792/1‑148 (MQ=255)                                     aCCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAgg                                                                                                                                         >  7:153483/1‑121 (MQ=255)                                     aCCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAgg                                                                                                                                         <  8:153483/121‑1 (MQ=255)                                         gCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGa                                                                                                         >  2:334929/1‑149 (MQ=255)                                           tGGAGTTGCTCGACCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAg                                                                                                                                          >  8:159471/1‑114 (MQ=255)                                            ggAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACgg                                                                                                      >  5:261415/1‑149 (MQ=255)                                                ttGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                  <  4:149094/149‑1 (MQ=255)                                                ttGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGa                                                                                                  <  5:362311/149‑1 (MQ=255)                                                  gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCa                                                                                                    <  6:192763/145‑1 (MQ=255)                                                  gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCa                                                                                                    >  5:192763/1‑145 (MQ=255)                                                  gCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTCACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAATCGACGGCAGAtt                                                                                                >  7:59443/1‑149 (MQ=255)                                                    tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                     >  8:189063/1‑126 (MQ=255)                                                    tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTgaga                                                                                                                     <  7:189063/126‑1 (MQ=255)                                                    tCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTc                                                                                              >  6:10558/1‑149 (MQ=255)                                                      gATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCg                                                                                             <  5:350724/148‑1 (MQ=255)                                                       aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGTGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                           <  7:16538/149‑1 (MQ=255)                                                       aTCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTa                                                                                           <  5:104792/149‑1 (MQ=255)                                                        tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                           <  6:1945/116‑1 (MQ=255)                                                        tCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGt                                                                                                                           >  5:1945/1‑116 (MQ=255)                                                          acaTAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCt                                                                                        <  5:10558/149‑1 (MQ=255)                                                                     gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                             >  5:351193/1‑149 (MQ=255)                                                                     gCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCa                                                                             >  4:191337/1‑149 (MQ=255)                                                                       aaaCCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACt                                                                           >  8:205523/1‑149 (MQ=255)                                                                            gggTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAg                                                                      <  6:351193/149‑1 (MQ=255)                                                                                   acaacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                               >  4:146077/1‑149 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                     <  1:102514/125‑1 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                     >  2:98523/1‑125 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                     >  2:102514/1‑125 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGAt                                                                                     <  1:98523/125‑1 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                             >  8:196437/1‑149 (MQ=255)                                                                                     aacaGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGc                                                             >  4:286549/1‑149 (MQ=255)                                                                                       caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                    <  8:272923/92‑1 (MQ=255)                                                                                       caGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGagag                                                                                                                    >  7:272923/1‑92 (MQ=255)                                                                                        aGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCa                                                          <  1:334929/149‑1 (MQ=255)                                                                                         gTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                            >  1:174590/1‑146 (MQ=255)                                                                                          ttCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGcc                                                            >  6:115934/1‑145 (MQ=255)                                                                                                   ctctTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCa                                               <  3:191337/149‑1 (MQ=255)                                                                                                      ttCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGtt                                            >  5:130761/1‑149 (MQ=255)                                                                                                        cGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTCACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGttt                                           <  8:59443/148‑1 (MQ=255)                                                                                                             tAGATGAGTAGCGTCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTc                                     <  6:261415/149‑1 (MQ=255)                                                                                                                     tACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                               >  7:56868/1‑115 (MQ=255)                                                                                                                     tACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTg                                                               <  8:56868/115‑1 (MQ=255)                                                                                                                        cGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGc                           >  4:77354/1‑148 (MQ=255)                                                                                                                          gCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTg                                        >  7:266200/1‑133 (MQ=255)                                                                                                                          gCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTg                                        <  8:266200/133‑1 (MQ=255)                                                                                                                             gcgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCgg                     <  3:77354/149‑1 (MQ=255)                                                                                                                              cgcgcCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCggg                    <  7:205523/149‑1 (MQ=255)                                                                                                                                  ccAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATg                <  6:130761/149‑1 (MQ=255)                                                                                                                                      tCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaa             >  2:46850/1‑148 (MQ=255)                                                                                                                                       cTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACTGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGaaaa           >  1:123515/1‑149 (MQ=255)                                                                                                                                          cccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                 >  4:151837/4‑140 (MQ=255)                                                                                                                                          cccgACCAGGCCACCTTGGGTGGGCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCAt                 <  3:151837/137‑1 (MQ=255)                                                                                                                                             gACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc      <  7:196437/148‑1 (MQ=255)                                                                                                                                             gACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGc      <  3:286549/148‑1 (MQ=255)                                                                                                                                                cAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGa  >  5:113706/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GGAAACTCCCAGGCCTTGAGGATCTTGTCCCCCAGCACCGGCTGGAGTTGCTCGATCACATAATGCAGGCAAACCGGGTCGGACAACAGTTCGTTGTGCTCTTCGGCATAGATGAGTACCGGCAGCGCGCCAATCTGGTTGACCAGCCCACCCAAGGTGGCCTGGTCGGGTTTGAGAGGTGTGAAGCGACGGCAGATTTCGTAGCTGATGCCCGCCACTTCCAGGCTACTGGCCCAGATATCGCGCAGTTTCTGCTCCACTGCCGGCGAGCGGGCATGGAAAATCTGCTCGA  >  NC_002947/1622761‑1623052

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