Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,164,586	A→T	12.1%	intergenic (+139/+62)	PP_4544 → / ← PP_4545	hypothetical protein/beta‑ketodecanoyl‑ACP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,164,586	0	A	T	12.1%	189.4 / 13.7	72	intergenic (+139/+62)	PP_4544/PP_4545	hypothetical protein/beta‑ketodecanoyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base A (35/25);  new base T (2/7);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 7.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

CCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCC  >  NC_002947/5164437‑5164735                                                                                                                                                       |                                                                                                                                                       ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGTAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGaa                                                                                                                                                        <  1:1497173/150‑1 (MQ=255) ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGaa                                                                                                                                                        <  1:4135542/150‑1 (MQ=255) ccAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGaa                                                                                                                                                        <  2:1454774/150‑1 (MQ=255)  cAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGTAAACGCCGCCCCCATCCGGGGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAAc                                                                                                                                                       <  1:1182344/150‑1 (MQ=255)      ttCTAAGCAACACCCTGCTCTGCTTCACCTGAAGGGGTAAACGCCGACCCCACCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACtgtg                                                                                                                                                   <  1:2743343/150‑1 (MQ=255)          aaGCAACACCCTGCTCTGCTTCACCTGAAGCGGTAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGt                                                                                                                                               <  1:2499833/150‑1 (MQ=255)                cacCCTGCTCTGCTTCACCTGAAGCGGTAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑tgt                                                                                                                                        <  1:178015/150‑1 (MQ=255)                 acCCTGCTGTGCTTCACCTGAAGCGGGAACCGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑tgtg                                                                                                                                       <  2:242534/150‑1 (MQ=255)                 acCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCATTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑tgtg                                                                                                                                       <  2:224708/150‑1 (MQ=255)                    ctgctctgctTCACCTGAAGGGTTAAAGGCCTCCCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGa                                                                                                                                    <  1:3695994/150‑1 (MQ=38)                      gctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGc                                                                                                                                  <  2:2048733/150‑1 (MQ=18)                      gctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGc                                                                                                                                  <  2:2048522/150‑1 (MQ=18)                      gctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGc                                                                                                                                  <  1:4913790/150‑1 (MQ=18)                        tctgctTCACCTGACGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGCgg                                                                                                                                <  1:4599766/150‑1 (MQ=17)                        tctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCgg                                                                                                                                <  1:2675459/150‑1 (MQ=255)                        tctgctACACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGCgg                                                                                                                                <  1:4612367/150‑1 (MQ=17)                             ttCACCTGAAGCGGGAACCGCCGACCCCCTCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTGGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAg                                                                                                                           <  2:1088606/150‑1 (MQ=255)                               cACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGcc                                                                                                                         <  2:3324716/150‑1 (MQ=255)                                  cTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGCGGGCAAGCCcgc                                                                                                                      <  2:2954478/150‑1 (MQ=17)                                         gggAAACGCCGACCCCAGCCGGCGTGTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCTTGCCcgctcccaca                                                                                                               <  2:2955875/150‑8 (MQ=37)                                                   gACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCGCCCACACT‑ATAACCACAGAACGTGAGGACGGGGGTGTCCC‑TGTGGGAGCGGGCCAGCCCGCGACGAGGGCGTAACAgg                                                                                                     >  2:3049097/1‑150 (MQ=255)                                                         aTCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGAACTTGAGAACAGTTCTGTAAC‑AGTTGGAGCGGGCAAGCCCCCGAAGAGGCCGTAACAGGCGCCtg                                                                                               >  2:4881667/1‑150 (MQ=39)                                                           cgggCGTTTCCCGCTTATTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGGTTGCCCGCTCCCACTGT‑ATCACCACAGACCTTGAGAAGTGTGCGTTACC‑TGTGGGAGGGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                                             <  2:2534438/148‑1 (MQ=37)                                                           ccGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCCCGCCCCCAGG‑AAAAACAAAGAACGTTATAACAGTTCTGGTCC‑TGTGGGGGCGGGCAAGCCCGCGAAGCGGCCGTGAACGGCGccggca                                                                                             >  2:767789/1‑146 (MQ=255)                                                           ccGGCGTTTCCCGCTTCTTTTGGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCtgat                                                                                             <  2:820848/150‑1 (MQ=255)                                                                gTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCCCTCCCACACG‑ATAACCCCACACCATGAGAACAGTTCTGTAAC‑AGTGGGTGCCGGCAACCCCCCCAAGAGGCCCGGACACGCGCCGGATGaaaa                                                                                        >  2:282831/1‑148 (MQ=255)                                                                     cctcttcttTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCACCCACACT‑ATGACCACAAACCTTGAGAACAGTGCTGTACC‑GGTGGGCGCGGGCAAGCCCCCGAACAGGCCGGAAAAGGCGCCTGAGGAATAAgcgc                                                                                   >  2:337657/4‑150 (MQ=255)                                                                     ccGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTTCCCGCGCCCACACT‑ATCACCACACACCATGAGAACTGTGCTGTAAC‑GGTGGGAGCGGGCCAGCCCGCGAAGAGGCCGGAACAGGCGCCGGATGGAACAgcgc                                                                                   >  2:5275017/1‑150 (MQ=255)                                                                     ccGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑AGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGGAACACGCGCCGGATGaaaaaccgc                                                                                   >  2:336919/1‑143 (MQ=37)                                                                       gCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCGTTCCCCCGCCCACAGT‑ATCACCACAGACCGTGAGAACTGTGCTGTACC‑AGTGGGAGCCGGCAAGCCCGCGCAGAGGCCGGAACAGGCGCCAGATGAATAAccgcgt                                                                                 >  2:394798/1‑148 (MQ=255)                                                                        ctacttTTTGACTGTGTCAAATGCGCCGCCCCCTTCGCGGGATTGCCCCCACCCACACT‑AACACCACACACCTTTAGAACTGTGCTTTACC‑TGTGGCAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTg                                                                                <  2:74309/147‑1 (MQ=37)                                                                             ttttGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTACCACAGG‑GTAACCCCAGAAAATGTGAACTGTGCGGTTCC‑GGGGGGAGCGGGCCAGCCCCCGGAGATGCCGGAAAAGGCGCCAGATGAAACAGCGCGGGCGCAg                                                                           >  2:2661309/1‑150 (MQ=255)                                                                              tttGCCTGTGTTAAAGGCGCCGGCCCCTTCGGGGGCTTGCCCGCTCCCCCAGT‑ACCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGCCGCGGCCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGa                                                                          <  1:1609363/150‑1 (MQ=255)                                                                                tGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCCCCCACCTT‑ATCACCACAGACCTTGTGAAGTGTGCTGTACC‑TGTGGGAGCGGGAAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAAt                                                                        <  2:1581539/150‑1 (MQ=39)                                                                                    tgtgtgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCCCCCACAGT‑ATCACCACCTCCCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATcaca                                                                    <  1:872019/150‑1 (MQ=255)                                                                                       gtgAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTa                                                                 <  2:4944011/150‑1 (MQ=255)                                                                                            aTGCGCCGGCCCCTTCGCGGGCGTGCCCGCTCCCACAGT‑ATCACCCCCCTCCGTGCGCCCGGTGCTGTGCC‑GGTGGGCGCGGGCAACCCCCCGAAGAGGCCGTCAACGGCGCCGGATGACTCAGCGCGGGCGCCGAAGCACCCGAcccct                                                            >  2:3581729/1‑147 (MQ=255)                                                                                             tGCGCCGGCCCCTTCGCGGGCTGCCCCGCTCCCACAGT‑AACACCACACACCTTGAGAACTGTGCGCTCCC‑GGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATc                                                           <  1:4881667/150‑1 (MQ=37)                                                                                                   ggCCCCGGCGCGGGGTGCCCCGC‑CCCACAATAACCACAACCTACTTTGAGACCTGTGC‑GTCCCGTGTGGCGGGGGCCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATAc                                                     <  1:2329215/150‑1 (MQ=255)                                                                                                       ccTTGGGGGGGTGGCCCGCCCCCACAGT‑AACACCACTTACTTTTAAATTTTGGCTTTCCC‑TGTGGGGGGGGGCACCCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCg                                                 <  1:3658250/150‑1 (MQ=255)                                                                                                        cTTCGCGGGCTTGCCCCCTCCCACAGT‑ATAAACACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCCCGAAGAGGCCGTAACAGGCGCCAGATTAAACACCGCTTTCGCAGAATCAAACAAACAATAGAATAACCCGc                                                >  2:4737297/1‑150 (MQ=25)                                                                                                           gggggggTTGCCCGCACCCCCATT‑AACACCACTGACATTTATAAGTGTGCTGTCCC‑GGTGGGCGGGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCAcc                                             <  1:1603034/147‑1 (MQ=255)                                                                                                            gcgGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGTTAC‑TGTGGGAGCGGGCAAG‑CCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAAAACCCCGCACCaa                                           >  1:4520055/1‑150 (MQ=17)                                                                                                               ggCTTGCCGGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TTTGGGTGGGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAAc                                         <  1:822363/150‑1 (MQ=255)                                                                                                                gCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACg                                        >  2:711754/1‑150 (MQ=255)                                                                                                                  ttGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAg                                      >  2:1741924/1‑150 (MQ=255)                                                                                                                  gtGCCGGCTCCCACTGT‑ATCACCACAGACGTTGTGAAGTGTGGTGTACG‑TGTGGGTGCGGGAATGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAg                                      <  2:1442797/149‑1 (MQ=255)                                                                                                                    gCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGCGGGCAAGCCCGCGAAGGGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCt                                    >  1:3791594/1‑150 (MQ=17)                                                                                                                     cccGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTc                                   >  2:3161660/1‑150 (MQ=255)                                                                                                                      ccGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                                  >  1:962518/1‑150 (MQ=255)                                                                                                                      ccGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCa                                  >  1:5105678/1‑150 (MQ=255)                                                                                                                        gCTCCCTCAGT‑ATCACCACTGACCTTGAGAAGTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAac                                <  2:4912793/150‑1 (MQ=255)                                                                                                                            ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                            >  2:798358/1‑150 (MQ=255)                                                                                                                            ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                            >  1:3558237/1‑150 (MQ=255)                                                                                                                            ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                            >  2:2170438/1‑150 (MQ=255)                                                                                                                            ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                            >  2:3611681/1‑150 (MQ=255)                                                                                                                             cacaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc                           >  1:2072868/1‑150 (MQ=255)                                                                                                                             cacaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc                           >  1:902629/1‑150 (MQ=255)                                                                                                                              acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                          >  1:4331470/1‑150 (MQ=255)                                                                                                                              acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                          >  2:3486274/1‑150 (MQ=255)                                                                                                                              acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                          >  1:2451145/1‑150 (MQ=255)                                                                                                                                     tCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAg                    >  1:2647475/1‑150 (MQ=255)                                                                                                                                        ccacAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCt                 >  1:3416680/1‑150 (MQ=255)                                                                                                                                         cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAAc                       >  2:1571198/1‑143 (MQ=255)                                                                                                                                         cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAAc                       <  1:1571198/143‑1 (MQ=255)                                                                                                                                         cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCtt                >  1:5260068/1‑150 (MQ=255)                                                                                                                                         cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCtt                >  1:4366316/1‑150 (MQ=255)                                                                                                                                         cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCtt                >  1:5186162/1‑150 (MQ=255)                                                                                                                                           attacgttGATAAGTGTGCTGTACG‑TGTGGGCGGGGGAACCCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTg              <  1:964525/144‑1 (MQ=255)                                                                                                                                             gcccTTGCGACGTGTGCTGTCCC‑TGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGc            <  1:2087593/148‑1 (MQ=255)                                                                                                                                               ccTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCa          >  2:3334486/1‑150 (MQ=255)                                                                                                                                                cTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAg         >  2:430170/1‑150 (MQ=255)                                                                                                                                                cTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAg         >  2:2215774/1‑150 (MQ=255)                                                                                                                                                 ttGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAgg        >  1:631091/1‑150 (MQ=255)                                                                                                                                                    agaACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCa     >  2:3997361/1‑150 (MQ=255)                                                                                                                                                      aaCTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATc   >  1:992744/1‑150 (MQ=255)                                                                                                                                                       aCTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATcc  <  1:394798/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CCAAGTTCTAAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCC  >  NC_002947/5164437‑5164735

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