Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,164,592	C→G	12.7%	intergenic (+145/+56)	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,592	0	C	G	12.7%	189.9 / 17.7	77	intergenic (+145/+56)	PP_4544/PP_4545	hypothetical protein/beta‑ketodecanoyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base C (39/28);  new base G (3/7);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

AAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTGGTAC  >  NC_002947/5164446‑5164741                                                                                                                                                    |                                                                                                                                                       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                                                                                                                                        <  1:4913790/150‑1 (MQ=18)             gctctgctTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGG‑TACAGCACAGTTCTCAAGGTCTGTGGTGATAC‑TGTGGGAGc                                                                                                                                        <  2:2048522/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:3611681/1‑150 (MQ=255)                                                                                                                   ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                                  >  1:3558237/1‑150 (MQ=255)                                                                                                                   ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                                  >  2:798358/1‑150 (MQ=255)                                                                                                                   ccACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGc                                  >  2:2170438/1‑150 (MQ=255)                                                                                                                    cacaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc                                 >  1:902629/1‑150 (MQ=255)                                                                                                                    cacaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGcc                                 >  1:2072868/1‑150 (MQ=255)                                                                                                                     acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                                >  1:4331470/1‑150 (MQ=255)                                                                                                                     acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                                >  1:2451145/1‑150 (MQ=255)                                                                                                                     acaGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGccc                                >  2:3486274/1‑150 (MQ=255)                                                                                                                            tCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAg                          >  1:2647475/1‑150 (MQ=255)                                                                                                                               ccacAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCt                       >  1:3416680/1‑150 (MQ=255)                                                                                                                                cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAAc                             <  1:1571198/143‑1 (MQ=255)                                                                                                                                cacaGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAAc                             >  2:1571198/1‑143 (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)                                                                                                                                                tgtgCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTg      >  1:1152440/1‑150 (MQ=255)                                                                                                                                                tgtgCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTg      >  1:507646/1‑150 (MQ=255)                                                                                                                                                ggggCTGCACG‑TGTGGGAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTg      <  1:3804578/147‑1 (MQ=255)                                                                                                                                                 gggCTGTACC‑TGTGGGAGCGGGCACGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTgg     <  1:2650841/148‑1 (MQ=255)                                                                                                                                                  tgCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAACTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTGGt    >  1:4147248/1‑150 (MQ=255)                                                                                                                                                    cTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTGGTAc  >  1:1350051/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       AAGCAACACCCTGCTCTGCTTCACCTGAAGCGGGAAACGCCGACCCCATCCGGCGTTTCCCGCTTCTTTTTGTCTGTGTGAAATGCGCCGGCCCCTTCGCGGGCTTGCCCGCTCCCACAGT‑ATCACCACAGACCTTGAGAACTGTGCTGTACC‑TGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGTAACAGGCGCCTGATGAATCAGCGCTTGCGCAGAATCACACTACCAATCGAATACCCCGCACCAAACGAGCTCAACACGCCCAACGAGCCTTTGGCCAGGTCATCCTGGTAC  >  NC_002947/5164446‑5164741

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