Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,133,977	T→A	9.2%	A351A (GCA→GCT) ‡	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,133,977	0	T	A	9.2%	250.5 / 12.2	87	A351A (GCA→GCT) ‡	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (37/42);  new base A (6/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAG  >  NC_002947/4133830‑4134117                                                                                                                                                      |                                                                                                                                             ccACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACgcctg                                                                                                                                             >  5:231143/1‑149 (MQ=255) ccACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACgcctg                                                                                                                                             <  3:110492/149‑1 (MQ=255)         cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                     <  4:151385/149‑1 (MQ=255)         cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                     <  4:156962/149‑1 (MQ=255)            gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGtt                                                                                                                                    <  6:231143/147‑1 (MQ=255)            gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   <  8:102555/148‑1 (MQ=255)            gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  <  5:147325/149‑1 (MQ=255)              gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   <  3:2570/146‑1 (MQ=255)              gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                   >  4:2570/1‑146 (MQ=255)               ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                               >  2:96139/1‑149 (MQ=255)               ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                               >  7:350261/1‑149 (MQ=255)                  aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  <  3:79704/143‑1 (MQ=255)                  aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                  >  4:79704/1‑143 (MQ=255)                  aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCg                                                                                                                            <  8:334271/149‑1 (MQ=255)                         tCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGgcgcgc                                                                                                                      >  2:365433/1‑148 (MQ=255)                              aCCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCCCGCTCA‑‑TGCAACGCCTGCCTCGCGTTGCAGTTCGGCGCGCAAttcc                                                                                                                <  4:7324/149‑3 (MQ=255)                                cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTg                                                                                                                 <  5:340406/146‑1 (MQ=255)                                cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTg                                                                                                                 >  6:340406/1‑146 (MQ=255)                                 tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             >  6:69861/1‑149 (MQ=255)                                 tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAc                                                                                                              >  5:109162/1‑148 (MQ=255)                                   aCCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGACCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAg                                                                                                           >  4:45746/1‑149 (MQ=255)                                    ccTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          <  7:46084/149‑1 (MQ=255)                                     cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCAGGCAGGCGTTGc                                                                                                                                  <  5:37234/124‑1 (MQ=255)                                     cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCAGGCAGGCGTTGc                                                                                                                                  >  6:37234/1‑124 (MQ=255)                                        gTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      <  6:129816/149‑1 (MQ=255)                                              gCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  4:288537/1‑149 (MQ=255)                                                  cAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGa                                                                                             <  2:75594/148‑1 (MQ=255)                                                   aGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATc                                                                                           <  2:277092/149‑1 (MQ=255)                                                     gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCgg                                                                                         <  6:109162/149‑1 (MQ=255)                                                        gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          <  3:3606/129‑1 (MQ=255)                                                        gAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                          >  4:3606/1‑129 (MQ=255)                                                         aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCAc                                                                                      >  5:232049/1‑148 (MQ=255)                                                         aaCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCAc                                                                                      >  7:71766/1‑148 (MQ=255)                                                               gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                <  5:145944/132‑1 (MQ=255)                                                               gcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  6:145944/1‑132 (MQ=255)                                                                  aGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGAc                                                                            <  3:288537/149‑1 (MQ=255)                                                                        tgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      >  6:356482/1‑117 (MQ=255)                                                                        tgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                      <  5:356482/117‑1 (MQ=255)                                                                           tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             >  8:167011/1‑107 (MQ=255)                                                                           tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                             <  7:167011/107‑1 (MQ=255)                                                                             ccAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcg                                                                 >  1:310069/1‑149 (MQ=255)                                                                               aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggc                                                               >  2:249085/1‑149 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              <  1:88964/149‑1 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              <  2:264600/149‑1 (MQ=255)                                                                                accaccCTCACGCGGGTAGCCAAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGATGGCCGACTCCCGCCATcggcc                                                              >  6:85999/1‑148 (MQ=255)                                                                                 ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  6:232049/149‑1 (MQ=255)                                                                                 ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  8:125234/149‑1 (MQ=255)                                                                                 ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                             <  8:130626/149‑1 (MQ=255)                                                                                  caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCAAGTGCAACGCATGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCAATcggcg                                                              >  1:293698/1‑149 (MQ=255)                                                                                   accCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                           <  8:350261/149‑1 (MQ=255)                                                                                    cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                >  2:47309/1‑111 (MQ=255)                                                                                    cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                <  1:47309/111‑1 (MQ=255)                                                                                         aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAccc                                                     >  5:115306/1‑149 (MQ=255)                                                                                         aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAccc                                                     <  1:365433/149‑1 (MQ=255)                                                                                         aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       <  5:260801/147‑1 (MQ=255)                                                                                               gTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               >  4:247677/1‑149 (MQ=255)                                                                                                  gCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTa                                            >  2:155607/1‑149 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGc                                                                                                         <  5:367897/86‑1 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGc                                                                                                         >  6:367897/1‑86 (MQ=255)                                                                                                    ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                          >  3:169542/1‑149 (MQ=255)                                                                                                                cACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               <  7:42597/132‑1 (MQ=255)                                                                                                                cACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                               >  8:42597/1‑132 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       >  5:8210/1‑123 (MQ=255)                                                                                                                 aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                       <  6:8210/123‑1 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                           >  6:283258/1‑129 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                           <  5:283258/129‑1 (MQ=255)                                                                                                                       cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgcc                       <  2:293698/149‑1 (MQ=255)                                                                                                                          cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     >  1:372880/1‑148 (MQ=255)                                                                                                                          cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     <  4:169542/148‑1 (MQ=255)                                                                                                                          cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccg                      <  2:310069/147‑1 (MQ=255)                                                                                                                          cAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGAGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgccgc                     >  8:215329/1‑148 (MQ=255)                                                                                                                            gCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  <  1:155607/149‑1 (MQ=255)                                                                                                                                    gcgcTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGc          <  6:115306/149‑1 (MQ=255)                                                                                                                                       ccgaactGCAACGCCAGGCAGGCGTTGCA‑‑TGAGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                >  8:126061/7‑140 (MQ=37)                                                                                                                                       ccgaactGCAACGCCAGGCAGGCGTTGCA‑‑TGAGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                <  7:126061/134‑1 (MQ=37)                                                                                                                                        tCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGAGCAATTGCACCAGGCTGGGCGTGAAGATCGGCACGAAGGCCGACTCCAGCCATCGGCGGGCAAACCCTTCCCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTCGCCTGCAACt      >  6:326020/1‑149 (MQ=255)                                                                                                                                        ggatctGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCaa        <  5:107969/144‑1 (MQ=255)                                                                                                                                          a‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                           <  7:190400/94‑1 (MQ=255)                                                                                                                                          a‑‑TGCAACGCCTGCCTGGCGTTCCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGAACGAAGGCCGACTCCCGCCAACGGCggga                                                           >  8:190400/1‑93 (MQ=255)                                                                                                                                           actGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  2:42243/3‑149 (MQ=255)                                                                                                                                             tGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCATGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCa   <  3:45746/149‑1 (MQ=255)                                                                                                                                             tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGACCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  1:347169/1‑147 (MQ=255)                                                                                                                                             tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGACCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  1:352840/1‑147 (MQ=255)                                                                                                                                              gtaaCGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                              >  8:298883/3‑89 (MQ=255)                                                                                                                                              gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  <  1:160829/133‑1 (MQ=255)                                                                                                                                              gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                  >  2:160829/1‑133 (MQ=255)                                                                                                                                              gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg  >  5:357925/1‑149 (MQ=255)                                                                                                                                              gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg  >  8:185992/1‑149 (MQ=255)                                                                                                                                               cAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGCTAAGCCTTGCCGCCACTGGCCTGCaa        >  6:107969/1‑142 (MQ=255)                                                                                                                                                   gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     <  8:121647/141‑1 (MQ=255)                                                                                                                                                   gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc     >  7:121647/1‑141 (MQ=255)                                                                                                                                                    cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                        >  5:77756/1‑73 (MQ=255)                                                                                                                                                    cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                        <  6:77756/73‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             CCACACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAG  >  NC_002947/4133830‑4134117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑