Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,864,750	T→C	5.6%	A28A (GCT→GCC)	hda →	DnaA regulatory inactivator Hda

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,864,750	0	T	C	5.6%	261.8 / 4.0	89	A28A (GCT→GCC)	hda	DnaA regulatory inactivator Hda
Reads supporting (aligned to +/- strand):  ref base T (45/39);  new base C (2/3);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CGGGTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTG  >  NC_002947/1864603‑1864898                                                                                                                                                    |                                                                                                                                                     cGGGTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAAtgctg                                                                                                                                                     <  3:48990/149‑1 (MQ=255)    gTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCgg                                                                                                                                                  <  1:125522/149‑1 (MQ=255)    gTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCAGGCGCCAATGCTGCgg                                                                                                                                                  >  7:277131/1‑149 (MQ=255)         gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                             <  4:186664/149‑1 (MQ=255)         gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                             <  7:11815/149‑1 (MQ=255)          aCTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTgg                                                                                                                                            >  1:60823/1‑149 (MQ=255)           cTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggc                                                                                                                                           <  3:237730/149‑2 (MQ=255)                cAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTAc                                                                                                                                       <  1:153019/148‑1 (MQ=255)                cAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCCGCAGCATTggcgccc                                                                                                                                       <  4:365304/148‑6 (MQ=255)                     tGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGAATGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAg                                                                                                                                 >  5:152566/1‑149 (MQ=255)                       ttCACCATAGGCCCGCAGCGGTCCGGTCGACCGCCCGCGAGCATCAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGc                                                                                                                                <  8:140934/148‑1 (MQ=255)                       ttCAACAAAGGCCCCGAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCGTGTTGCCCCTGGGTGTGCGCCTGCGCGATGTCGCCACGTTGATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGtggtgc                                                                                                                                <  8:135865/148‑6 (MQ=255)                         cACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGc                                                                                                                             >  6:141990/1‑149 (MQ=255)                          aCCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGc                                                                                                                             >  1:110793/1‑148 (MQ=255)                          aCCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGACACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCt                                                                                                                            >  6:373501/1‑149 (MQ=255)                            cAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTAt                                                                                                                          >  1:197405/1‑149 (MQ=255)                               aGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGc                                                                                                                        >  4:366138/1‑148 (MQ=255)                               aGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGc                                                                                                                        >  2:195174/1‑148 (MQ=255)                                 gCCCCCAGCGGTCCCGTCGACCGCCAGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGaa                                                                                                                     <  8:342621/149‑1 (MQ=255)                                   ccccAGCGGTCCCGTCGACCGCCCGCGCGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGc                                                                                                                   >  8:269786/1‑149 (MQ=255)                                     ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCg                                                                                                                                   <  8:189440/131‑1 (MQ=255)                                     ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCg                                                                                                                                   >  7:189440/1‑131 (MQ=255)                                     ccAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGcc                                                                                                                  >  5:375246/1‑148 (MQ=255)                                      cAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGa                                                                                                                <  2:60823/149‑1 (MQ=255)                                            tCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCgg                                                                                                          >  6:307381/1‑149 (MQ=255)                                                   gACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGAcc                                                                                                   <  7:268948/149‑1 (MQ=255)                                                    aCCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGAcc                                                                                                   >  8:307077/1‑148 (MQ=255)                                                     ccgcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCga                                                                                                 <  5:141990/149‑1 (MQ=255)                                                     ccgcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTCCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCga                                                                                                 <  5:333763/149‑1 (MQ=255)                                                      cgcccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCgag                                                                                                <  5:298823/149‑1 (MQ=255)                                                         ccgcGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGc                                                                                             <  2:110793/149‑1 (MQ=255)                                                          cgcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                            <  1:290841/149‑1 (MQ=255)                                                          cgcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                            <  7:307077/149‑1 (MQ=255)                                                           gcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCt                                                                                           >  1:27252/1‑149 (MQ=255)                                                           gcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCt                                                                                           <  2:197405/149‑1 (MQ=255)                                                                 aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCg                                                                                                                 >  3:266298/1‑121 (MQ=255)                                                                 aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCg                                                                                                                 <  4:266298/121‑1 (MQ=255)                                                                 aTGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACAACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACC‑‑GAGCCTTATCTAcc                                                                                   >  7:244239/1‑149 (MQ=255)                                                                  tGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAc                                                                                    <  6:375246/149‑1 (MQ=255)                                                                   gAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAcc                                                                                   >  6:216613/1‑149 (MQ=255)                                                                   gAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTAcc                                                                                   >  1:359317/1‑149 (MQ=255)                                                                         accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTTCCTGTggg                                                                             >  3:271878/1‑149 (MQ=255)                                                                         accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                             <  7:269786/149‑1 (MQ=255)                                                                         accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                             >  2:106798/1‑149 (MQ=255)                                                                         accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                             >  2:243558/1‑149 (MQ=255)                                                                                  gTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAgg                                                                    <  2:27252/149‑1 (MQ=255)                                                                                     gCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGgtg                                                                 >  6:172437/1‑149 (MQ=255)                                                                                      ccccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGgtg                                                                 <  7:141367/148‑1 (MQ=255)                                                                                       cccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                   >  2:98811/1‑65 (MQ=255)                                                                                       cccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                   <  1:98811/65‑1 (MQ=255)                                                                                       cccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGtt                                                               >  5:256456/1‑149 (MQ=255)                                                                                        ccggggTGTGCGCCTGCGCGATGACGCCACCTGCACCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGCTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTAATCTACCTCTGGGGCAAGCAGGGTGTTg                                                              <  6:98004/146‑1 (MQ=255)                                                                                        ccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTg                                                              >  5:127855/1‑149 (MQ=255)                                                                                        ccTGGGGGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTg                                                              <  4:178335/149‑1 (MQ=255)                                                                                         cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                             >  2:6596/1‑149 (MQ=255)                                                                                         cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                             <  1:195174/149‑1 (MQ=255)                                                                                         cTGGGAATGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                             <  6:152566/149‑1 (MQ=255)                                                                                              tgtgCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGtt                                                        >  8:97945/1‑149 (MQ=255)                                                                                               gtgCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTc                                                       <  7:353345/149‑1 (MQ=255)                                                                                                tgCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTTGCCGTTcc                                                      <  5:373501/149‑1 (MQ=255)                                                                                                  cgcCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCa                                                    >  7:288815/1‑149 (MQ=255)                                                                                                   gcCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCAc                                                   >  2:119880/1‑149 (MQ=255)                                                                                                      tGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACctg                                                >  3:164490/1‑149 (MQ=255)                                                                                                           gATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGt                                                                                                                                     <  3:315814/59‑1 (MQ=255)                                                                                                           gATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGt                                                                                                                                     >  4:315814/1‑59 (MQ=255)                                                                                                           gATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCa                                           <  6:256456/149‑1 (MQ=255)                                                                                                            aTGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAg                                          <  1:6596/149‑1 (MQ=255)                                                                                                                cGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgccg                                      <  5:307381/149‑1 (MQ=255)                                                                                                                cGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgccg                                      <  8:219928/149‑1 (MQ=255)                                                                                                                      cTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCAcc                                 >  5:120780/1‑148 (MQ=255)                                                                                                                          atcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATc                             <  6:127855/148‑1 (MQ=255)                                                                                                                          atcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCg                            >  3:188682/1‑149 (MQ=255)                                                                                                                           tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGc                           >  5:325036/1‑149 (MQ=255)                                                                                                                           tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGc                           >  5:157968/1‑149 (MQ=255)                                                                                                                           tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGc                           <  1:119880/149‑1 (MQ=255)                                                                                                                              actactATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTc                        >  8:214045/1‑149 (MQ=255)                                                                                                                                tactaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTc                        >  6:120339/1‑147 (MQ=255)                                                                                                                                  ctaTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAg                     >  6:339802/1‑148 (MQ=255)                                                                                                                                   taTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGc                    >  6:49749/1‑148 (MQ=255)                                                                                                                                    aTCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCaa                  >  1:18935/1‑149 (MQ=255)                                                                                                                                     tCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAAc                 >  6:182765/1‑149 (MQ=255)                                                                                                                                        gggCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTg              <  3:366138/149‑1 (MQ=255)                                                                                                                                         ggCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGGCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTgg             >  1:83621/1‑149 (MQ=255)                                                                                                                                          gcgcCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACC‑‑GAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGa          <  8:244239/149‑1 (MQ=255)                                                                                                                                            cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGc            >  6:194975/2‑147 (MQ=255)                                                                                                                                            cccAATGCCGCAGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGATGCAGGCCGCCACCCATCGATTCAAGCAACGTGCCGa          >  6:360911/2‑149 (MQ=255)                                                                                                                                            cccAATGCCGCAGCATTGGGTTACGGCGAGCGGCTATGCGTATCCGTCGCCGGCTGGACCGAGATCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGc            <  5:360911/146‑1 (MQ=255)                                                                                                                                              cAATGCTGCGGCTTTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGc        >  2:125813/1‑149 (MQ=255)                                                                                                                                              cAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGc        <  1:243558/149‑1 (MQ=255)                                                                                                                                               aaTGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGCGTGTTCGCCGATCCCACCTCCTGCAGGCCGCAACCCATCGCTTCCAGCAACGTGGCGAgac       >  6:200709/1‑147 (MQ=255)                                                                                                                                                    tgcGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTg  >  7:308043/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGGGTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTG  >  NC_002947/1864603‑1864898

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