Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,133,981	T→A	9.2%	Q350L (CAG→CTG)	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,981	0	T	A	9.2%	249.9 / 12.7	87	Q350L (CAG→CTG)	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 = 7.88e-01

CTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGT  >  NC_002947/4133838‑4134128                                                                                                                                                  |                                                                                                                                                    cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                                <  4:156962/149‑1 (MQ=255) cTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGt                                                                                                                                                <  4:151385/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                                                                                                                                              >  4:2570/1‑146 (MQ=255)      gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTg                                                                                                                                              <  3:2570/146‑1 (MQ=255)       ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                                          >  7:350261/1‑149 (MQ=255)       ccAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                                          >  2:96139/1‑149 (MQ=255)          aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                             >  4:79704/1‑143 (MQ=255)          aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGc                                                                                                                                             <  3:79704/143‑1 (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                                                                                                                            >  6:340406/1‑146 (MQ=255)                        cTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGTAGTTCGGCGCGCAATTg                                                                                                                            <  5:340406/146‑1 (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                                                                                                                                             >  6:37234/1‑124 (MQ=255)                             cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCAGGCAGGCGTTGc                                                                                                                                             <  5:37234/124‑1 (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                                                                                                                 <  5:356482/117‑1 (MQ=255)                                                                tgctgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTgg                                                                                                                 >  6:356482/1‑117 (MQ=255)                                                                   tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                                        <  7:167011/107‑1 (MQ=255)                                                                   tgcCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAcc                                                                                                                        >  8:167011/1‑107 (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                                                                        <  8:130626/149‑1 (MQ=255)                                                                         ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                        <  6:232049/149‑1 (MQ=255)                                                                         ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                        <  8:125234/149‑1 (MQ=255)                                                                          caccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCAAGTGCAACGCATGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCAATcggcg                                                                         >  1:293698/1‑149 (MQ=255)                                                                           accCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGc                                                                      <  8:350261/149‑1 (MQ=255)                                                                            cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                           <  1:47309/111‑1 (MQ=255)                                                                            cccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                           >  2:47309/1‑111 (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                                                          >  8:42597/1‑132 (MQ=255)                                                                                                        cACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGcc                                                          <  7:42597/132‑1 (MQ=255)                                                                                                         aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                  >  5:8210/1‑123 (MQ=255)                                                                                                         aCGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAAc                                                                  <  6:8210/123‑1 (MQ=255)                                                                                                               cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                                      <  5:283258/129‑1 (MQ=255)                                                                                                               cTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAt                                                      >  6:283258/1‑129 (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:352840/1‑147 (MQ=255)                                                                                                                                     tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGTGTGACGATCGGCACGAAGACCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc                >  1:347169/1‑147 (MQ=255)                                                                                                                                      gtaaCGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                         >  8:298883/3‑89 (MQ=255)                                                                                                                                      gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                             >  2:160829/1‑133 (MQ=255)                                                                                                                                      gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                             <  1:160829/133‑1 (MQ=255)                                                                                                                                      gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg             >  5:357925/1‑149 (MQ=255)                                                                                                                                      gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg             >  8:185992/1‑149 (MQ=255)                                                                                                                                       cAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGCTAAGCCTTGCCGCCACTGGCCTGCaa                   >  6:107969/1‑142 (MQ=255)                                                                                                                                           gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc                >  7:121647/1‑141 (MQ=255)                                                                                                                                           gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTc                <  8:121647/141‑1 (MQ=255)                                                                                                                                            cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                                   <  6:77756/73‑1 (MQ=255)                                                                                                                                            cctgcctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTccc                                                                                   >  5:77756/1‑73 (MQ=255)                                                                                                                                                 ctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGt  >  1:183012/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCA‑‑TGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGT  >  NC_002947/4133838‑4134128

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