Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	648,291	G→T	8.3%	intergenic (+60/+31)	acoR → / ← accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	648,291	0	G	T	8.3%	219.8 / 8.7	72	intergenic (+60/+31)	acoR/accC	acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base G (26/40);  new base T (3/3);  total (29/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

GCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGAC  >  NC_002947/648148‑648438                                                                                                                                                 |                                                                                                                                                    gCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCgg                                                                                                                                                <  2:53730/149‑1 (MQ=255)  cACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCggg                                                                                                                                               <  4:193980/149‑1 (MQ=255)  cACGGGCTGTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGCGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCggg                                                                                                                                               <  4:297328/149‑1 (MQ=255)   aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCGCGGCGCGGGCCCCGCTTGTGGGCGGGAATAAAAAACAACCGCACA‑GGGCCgggg                                                                                                                                              >  3:319721/1‑149 (MQ=255)   aCGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCGGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCT‑GTGCCgggg                                                                                                                                              <  3:244655/149‑1 (MQ=255)       gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCT‑GTGCCGGGGtttt                                                                                                                                          <  2:56594/149‑3 (MQ=255)       gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCt                                                                                                                                          <  6:223552/149‑1 (MQ=255)       gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCt                                                                                                                                          <  1:39516/149‑1 (MQ=255)           tGGGCATTTCGCGGGCGACGCTTTATCGGCGTGTGCATGAGCACCGCATCGACATGCCGCGTATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGt                                                                                                                                       <  1:201398/148‑1 (MQ=255)           tGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTgg                                                                                                                                        <  5:119716/147‑1 (MQ=255)           ggggCATTTCGCGGGCGACGCTTAATCGGCTAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGt                                                                                                                                       <  7:265541/147‑1 (MQ=255)               cATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGGCCCGCCCCGGGCCCAGATTCTGGCACGCAAAAAAAAACCCCCGCACA‑GGGCCGGGGGTCTGGTTCATc                                                                                                                                  >  5:177509/1‑149 (MQ=255)               cATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATc                                                                                                                                  <  4:44438/149‑1 (MQ=255)                aTTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCTCATCGACATCCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATca                                                                                                                                 <  7:69422/149‑1 (MQ=255)                aTTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATca                                                                                                                                 <  2:307900/149‑1 (MQ=255)                   tCGCGGGCGACGCTTTATCGGCGAGTGCAGGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCGGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATcaca                                                                                                                               <  5:292109/148‑1 (MQ=255)                         gCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCCCGCCGCGGCCCCAGAGTCGGGGAGGCAAATAAAAACACCCGCACACGGGCCGGGG‑TATGGGacaacacacaacac                                                                                                                         >  6:406840/1‑148 (MQ=255)                          cGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACATGa                                                                                                                       <  1:42372/149‑1 (MQ=255)                           gACGCTTTAGCGCCGCGTGCATGAGCACCGCATCGACAGGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGAt                                                                                                                      <  2:437827/149‑1 (MQ=255)                                       gCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                          <  7:331379/149‑1 (MQ=255)                                             gCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGCCGCGGGCCCAGATTCGGGGCGGCAATAAAAAACCCCCGCACA‑GGGCCGGGGTTGTGGGTCATCACACTACACGAGAACTACAGGTTgggccgt                                                                                                    >  4:54184/1‑144 (MQ=255)                                                    cACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGcc                                                                                                       <  5:358037/139‑1 (MQ=255)                                                     aCCGCATCTACATGCCGCGGATGATGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTc                                                                                            <  2:184991/149‑1 (MQ=255)                                                       cGCATCGACATGCCGCGGATGACGGGGTCGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                          <  4:287606/149‑1 (MQ=255)                                                       cGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGa                                                                                          <  1:292562/149‑1 (MQ=255)                                                         cATCGAAATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAgg                                                                                        <  7:73537/149‑1 (MQ=255)                                                           tCGACATGCGGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTa                                                                                      <  7:124381/149‑1 (MQ=255)                                                             gccATGCCGCGGATGAAGGGGTAGCGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGt                                                                                    <  7:61793/147‑1 (MQ=255)                                                             gACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTgg                                                                                                             <  5:325903/124‑1 (MQ=255)                                                             gACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAg                                                                                     <  8:76459/148‑1 (MQ=255)                                                             gACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGt                                                                                    <  8:238698/149‑1 (MQ=255)                                                                   ccGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGt                                                                              <  8:46136/149‑1 (MQ=255)                                                                    cgcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                             <  1:133894/149‑1 (MQ=255)                                                                            aaGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAcgccg                                                                      <  1:391239/148‑1 (MQ=255)                                                                                 gTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                                <  4:319721/149‑1 (MQ=255)                                                                                       cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAg                                                          >  5:379632/1‑149 (MQ=255)                                                                                       cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTAAACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAg                                                          >  3:1323/1‑149 (MQ=255)                                                                                       cccTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGGTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAg                                                          >  3:127509/1‑149 (MQ=255)                                                                                              gCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTGGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcat                                                   >  7:156235/1‑149 (MQ=255)                                                                                                cgcgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTtcatca                                                 <  6:177509/149‑1 (MQ=255)                                                                                                  cgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACg                                               >  1:202596/1‑149 (MQ=255)                                                                                                  cgACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACg                                               >  4:16377/1‑149 (MQ=255)                                                                                                            ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcg                                     >  8:84254/1‑149 (MQ=255)                                                                                                            ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcg                                     >  5:132939/1‑149 (MQ=255)                                                                                                            ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACAAGGTcg                                     >  6:368249/1‑149 (MQ=255)                                                                                                              aGATTCTGGCAGGCAATAAAAAACCCCGGCAGA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgcg                                   <  3:44480/149‑1 (MQ=255)                                                                                                                 ttCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCt                                                       >  8:193993/1‑126 (MQ=255)                                                                                                                 ttCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCt                                                       <  7:193993/126‑1 (MQ=255)                                                                                                                 ttCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTg                                <  2:202596/149‑1 (MQ=255)                                                                                                                     ggcaggcaATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                            >  6:338914/1‑149 (MQ=255)                                                                                                                        aggcaATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                            <  3:415069/146‑1 (MQ=255)                                                                                                                        aggcaATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGc                            >  4:415069/1‑146 (MQ=255)                                                                                                                         ggcaATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGa                        <  3:16377/149‑1 (MQ=255)                                                                                                                          gcaATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGAt                       >  2:6524/1‑149 (MQ=255)                                                                                                                           caATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATg                      <  5:195698/149‑1 (MQ=255)                                                                                                                            aaTAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                     <  6:132939/149‑1 (MQ=255)                                                                                                                            aaTAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                     >  1:209901/1‑149 (MQ=255)                                                                                                                               aaaaaaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTgg                  >  8:114430/1‑149 (MQ=255)                                                                                                                                 aaaaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                     >  3:387493/1‑144 (MQ=255)                                                                                                                                 aaaaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGt                     <  4:387493/144‑1 (MQ=255)                                                                                                                                 aaaaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTc                >  3:14181/1‑149 (MQ=255)                                                                                                                                  gaaCCCCGGCCCT‑GTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt               >  8:409591/2‑149 (MQ=255)                                                                                                                                  gaaCCCCGGCCCT‑GTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt               >  8:403940/2‑149 (MQ=255)                                                                                                                                  aaaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGGTGGTCt               >  8:346041/1‑149 (MQ=255)                                                                                                                                   aaCCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCtt              <  5:338914/149‑1 (MQ=255)                                                                                                                                     ccccGGCCCT‑GTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt                 <  5:305837/144‑1 (MQ=255)                                                                                                                                     ccccGGCCCT‑GTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGt                 >  6:305837/1‑144 (MQ=255)                                                                                                                                     ccccGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa            >  3:365346/1‑149 (MQ=255)                                                                                                                                      cccGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGAt           >  4:288006/1‑149 (MQ=255)                                                                                                                                        cGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGc         >  6:248153/1‑149 (MQ=255)                                                                                                                                           caca‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGt      <  5:272847/149‑1 (MQ=255)                                                                                                                                              a‑GGGCCGGGGTTCTGGTACATCACACTACAGGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTcga   >  1:36018/1‑149 (MQ=255)                                                                                                                                                gggCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTcgac  <  4:365346/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACA‑GGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCCGTCGAC  >  NC_002947/648148‑648438

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