Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,105,663	T→C	8.6%	S817G (AGC→GGC)	ydiJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,105,663	0	T	C	8.6%	197.2 / 7.3	70	S817G (AGC→GGC)	ydiJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base C (2/4);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCC  >  NC_002947/5105519‑5105811                                                                                                                                                 |                                                                                                                                                     gATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTAc                                                                                                                                                  >  4:229707/1‑149 (MQ=255)    acacCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgc                                                                                                                                               <  7:340553/149‑1 (MQ=255)      acCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcgc                                                                                                                                             >  3:243100/1‑149 (MQ=255)      acCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGACGTACAgcgc                                                                                                                                             >  3:161183/1‑149 (MQ=255)        cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcgc                                                                                                                                             >  2:284079/1‑147 (MQ=255)         aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                          <  3:57579/149‑2 (MQ=255)          gCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaac                                                                                                                                         <  8:261727/149‑3 (MQ=255)          gCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaac                                                                                                                                         <  4:161183/149‑3 (MQ=255)              cATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGCTGGCGCCGTACAGCGCCaacggca                                                                                                                                     <  7:332303/149‑7 (MQ=255)                   ccGGGTCCTGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCgg                                                                                                                                >  6:52991/1‑149 (MQ=255)                   ccGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCgg                                                                                                                                >  7:1354/1‑149 (MQ=255)                   ccGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACCCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCgg                                                                                                                                >  1:267832/1‑149 (MQ=255)                    cGGGTCCCGCCCCACCTGCGCCACACCGCAAGCGGCCAGGGCCTTGAGAGGCGTGGGGTTGCATAAGGCTGCATTGGCAAAGCGGCCAAGGAAACCCTGTACGTCCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGt                                                                                                                               <  1:76646/149‑1 (MQ=255)                        tCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGcc                                                                                                                           >  4:127776/1‑149 (MQ=255)                          cAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGccc                                                                                                                          <  5:188554/148‑1 (MQ=255)                             ccccACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTg                                                                                                                              >  3:9172/1‑141 (MQ=255)                             ccccACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTg                                                                                                                              <  4:9172/141‑1 (MQ=255)                             ccccACCAGCGGCACAACGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGACCAGGc                                                                                                                      >  2:322997/1‑149 (MQ=255)                               ccaccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCgg                                                                                                                    >  6:124245/1‑149 (MQ=255)                                 accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTg                                                                                                                  >  2:346772/1‑149 (MQ=255)                                        gCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAgg                                                                                                           <  5:124245/149‑1 (MQ=255)                                         cacaCCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGt                                                                                                          <  8:255772/149‑1 (MQ=255)                                            acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGa                                                                                                       <  3:229707/149‑1 (MQ=255)                                              cGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                      >  8:220180/1‑148 (MQ=255)                                              cGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                      >  6:246406/1‑148 (MQ=255)                                                      ggctgggCCTTGAGCTGCGTGGAGTTGCGGATCTCTTCGTTGGCAAACGCGCCAAGGAAAGGCTGTACTTGCAACGGCTTGCCGTTGGCGATGTACGGCGCCAGTAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCTATGAAGGCCGa                                                                                             <  8:273393/145‑1 (MQ=255)                                                             ccTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAc                                                                                            >  6:87101/1‑143 (MQ=255)                                                             ccTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAc                                                                                            <  5:87101/143‑1 (MQ=255)                                                               ttGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAgcgg                                                                                    >  5:359065/1‑149 (MQ=255)                                                                 gAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAgcggc                                                                                   <  5:219616/148‑1 (MQ=255)                                                                  aGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGt                                                                                 >  2:206215/1‑149 (MQ=255)                                                                             gTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGc                                                                      >  6:75090/1‑149 (MQ=255)                                                                                gCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCggg                                                                   >  8:112146/1‑149 (MQ=255)                                                                                     tCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                              <  2:229168/149‑1 (MQ=255)                                                                                     tCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                              <  5:246406/149‑1 (MQ=255)                                                                                     tCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                              >  1:184840/1‑149 (MQ=255)                                                                                       gCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGc                                                            <  3:127776/149‑1 (MQ=255)                                                                                        cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    <  3:80198/140‑1 (MQ=255)                                                                                        cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    >  4:80198/1‑140 (MQ=255)                                                                                        cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCAGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCg                                                           >  4:232193/1‑149 (MQ=255)                                                                                           ccTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTcc                                                        <  2:355138/149‑1 (MQ=255)                                                                                            cTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                  <  7:324226/106‑1 (MQ=255)                                                                                            cTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                  >  8:324226/1‑106 (MQ=255)                                                                                            cTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCt                                                       <  3:232193/149‑1 (MQ=255)                                                                                               ggCAAACGCGCCAAGGAATGCCTGTACGTGCAACGGCTTGGCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGc                                                     <  6:359065/148‑1 (MQ=255)                                                                                               ggCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCa                                                    >  3:261282/1‑149 (MQ=255)                                                                                                gCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCggg                                                                   >  5:94946/1‑133 (MQ=255)                                                                                                gCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCggg                                                                   <  6:94946/133‑1 (MQ=255)                                                                                                 cAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCAcc                                                  <  8:1354/149‑1 (MQ=255)                                                                                                  aaaCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCa                                                 >  4:361044/1‑149 (MQ=255)                                                                                                    aCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGc                                               <  5:52991/149‑1 (MQ=255)                                                                                                         ccAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGAt                                          >  3:15510/1‑149 (MQ=255)                                                                                                          cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCt                                                       >  5:82271/1‑135 (MQ=255)                                                                                                          cAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCt                                                       <  6:82271/135‑1 (MQ=255)                                                                                                                 aCCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCt                                  >  7:120399/1‑149 (MQ=255)                                                                                                                   ccTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcg                                >  5:72517/1‑149 (MQ=255)                                                                                                                    cTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcgc                               >  1:24778/1‑149 (MQ=255)                                                                                                                    cTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcgc                               <  1:284079/149‑1 (MQ=255)                                                                                                                        aCGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTg                           >  8:157406/1‑149 (MQ=255)                                                                                                                            gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                              >  3:290172/1‑110 (MQ=255)                                                                                                                            gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                              <  4:290172/110‑1 (MQ=255)                                                                                                                              aaCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCgg                     >  4:146002/1‑149 (MQ=255)                                                                                                                               aCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                    <  6:238860/149‑1 (MQ=255)                                                                                                                                cGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCgggg                   >  8:3021/1‑149 (MQ=255)                                                                                                                                    ttGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGt                  >  7:153651/1‑146 (MQ=255)                                                                                                                                    ttGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGt                  <  8:153651/146‑1 (MQ=255)                                                                                                                                       ttcctGGCGCCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTc            >  3:245673/5‑149 (MQ=255)                                                                                                                                             gcgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAg      <  7:157406/149‑1 (MQ=255)                                                                                                                                              cgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAgg     <  1:346772/149‑1 (MQ=255)                                                                                                                                               gcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAggg    >  7:204061/1‑149 (MQ=255)                                                                                                                                                 tGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGcc  >  4:336965/1‑149 (MQ=255)                                                                                                                                                 tGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGcc  >  5:127525/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCC  >  NC_002947/5105519‑5105811

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