Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,864,753	G→A	5.7%	A29A (GCG→GCA)	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,753	0	G	A	5.7%	262.6 / 4.2	87	A29A (GCG→GCA)	hda	DnaA regulatory inactivator Hda
Reads supporting (aligned to +/- strand):  ref base G (44/38);  new base A (2/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTG  >  NC_002947/1864606‑1864898                                                                                                                                                    |                                                                                                                                                  gTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCgg                                                                                                                                                  <  1:125522/149‑1 (MQ=255) gTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCAGGCGCCAATGCTGCgg                                                                                                                                                  >  7:277131/1‑149 (MQ=255)      gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                             <  7:11815/149‑1 (MQ=255)      gACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTg                                                                                                                                             <  4:186664/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                                                                                            <  7:307077/149‑1 (MQ=255)                                                       cgcgAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGcc                                                                                            <  1:290841/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:243558/1‑149 (MQ=255)                                                                      accGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTggg                                                                             >  2:106798/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                                                                                                                                                   <  1:98811/65‑1 (MQ=255)                                                                                    cccTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCg                                                                                                                                                   >  2:98811/1‑65 (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                                                             <  1:195174/149‑1 (MQ=255)                                                                                      cTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTgg                                                             >  2:6596/1‑149 (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                                      <  8:219928/149‑1 (MQ=255)                                                                                                             cGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGgccg                                      <  5:307381/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:157968/1‑149 (MQ=255)                                                                                                                        tcaACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGc                           >  5:325036/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)                                                                                                                                                   cGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTg  <  5:172437/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GTATAGACTTTGCAGATTGTTCACCAAAGGCCCCCAGCGGTCCCGTCGACCGCCCGCGAGCATGAAACCACCGATCCAGTTGCCCCTGGGTGTGCGCCTGCGCGATGACGCCACCTTCATCAACTACTATCCGGGCGCCAATGCTGCGGCATTGGGTTACGTCGAGCGGCTATGCGAAGCCGACGCCGGCTGGACCGAGAGCCTTATCTACCTGTGGGGCAAGCAGGGTGTTGGCCGTTCCCACCTGCTGCAGGCCGCCACCCATCGCTTCCAGCAACGTGGCGAGCCCGCTG  >  NC_002947/1864606‑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: ‑