Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,411,648	G→T	5.8%	A10A (GCC→GCA)	PP_3895 ←	regulatory protein Cro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,411,648	0	G	T	5.8%	209.1 / 3.0	70	A10A (GCC→GCA)	PP_3895	regulatory protein Cro
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base T (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GAGCCTCCTTCTTGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTAAAGACCGAATCAAA  >  NC_002947/4411500‑4411788                                                                                                                                                     |                                                                                                                                             gagCCTCCTTCATGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAg                                                                                                                                              <  6:294197/149‑1 (MQ=255)  agCCTCCTTCTTGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAgg                                                                                                                                             <  1:3195/149‑1 (MQ=255)    cctcctTCTTGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGc                                                                                                                                            >  2:71010/1‑148 (MQ=255)     ctcctTCTTGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGcga                                                                                                                                          >  8:310446/1‑149 (MQ=255)        cttcttGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGcgacga                                                                                                                                       >  5:22034/1‑149 (MQ=255)        cttcttGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGcgacga                                                                                                                                       <  6:204736/149‑1 (MQ=255)         ttcttGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAg                                                                                                                                      >  5:340675/1‑149 (MQ=255)         ttcttGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAg                                                                                                                                      >  6:176859/1‑149 (MQ=255)           cttGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGtt                                                                                                                                    >  6:136019/1‑149 (MQ=255)               aaTTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCttt                                                                                                                                <  5:335721/149‑1 (MQ=255)               aaTTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCttt                                                                                                                                >  6:3400/1‑149 (MQ=255)                 tttGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATATGCGACGAGTTCTTTg                                                                                                                               >  8:64407/1‑148 (MQ=255)                 tttGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGt                                                                                                                              <  5:3400/149‑1 (MQ=255)                 tttGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCAGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATATGCGACGAGTTCTTTg                                                                                                                               >  6:96222/1‑148 (MQ=255)                     cccTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGa                                                                                                                           <  5:43394/148‑1 (MQ=255)                      ccTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGa                                                                                                                           >  7:116704/1‑147 (MQ=255)                      ccTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATa                                                                                                                         <  5:125521/149‑1 (MQ=255)                       cTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAg                                                                                                                        >  5:102109/1‑149 (MQ=255)                          agtcagGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTg                                                                                                                      <  2:7515/148‑1 (MQ=255)                          agtcagGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGc                                                                                                                     <  5:239517/149‑1 (MQ=255)                              agGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCAt                                                                                                                 >  8:227571/1‑149 (MQ=255)                                gCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATgg                                                                                                               >  5:357325/1‑149 (MQ=255)                                gCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAATGCTCATgg                                                                                                               <  7:310446/149‑1 (MQ=255)                                        cTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCACGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCcctc                                                                                                       >  5:12301/1‑149 (MQ=255)                                               aGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATg                                                                                                 <  8:134795/148‑1 (MQ=255)                                                gCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATATGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATggg                                                                                               <  7:64407/149‑1 (MQ=255)                                                gCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATgg                                                                                                <  8:116704/148‑1 (MQ=255)                                                gCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCAGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATATGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATggg                                                                                               <  5:96222/149‑1 (MQ=255)                                                   accacAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAAt                                                                                            >  2:48292/1‑149 (MQ=255)                                                    ccacAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAAtt                                                                                           >  5:282509/1‑149 (MQ=255)                                                     cacaGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATtc                                                                                          >  1:185800/1‑149 (MQ=255)                                                      acaGGAGACATCCCGTGCTTACCGCTCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATg                                                                                                 <  1:48007/141‑1 (MQ=255)                                                      acaGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATg                                                                                                 >  2:48007/1‑141 (MQ=255)                                                      acaGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATtct                                                                                         <  4:388707/149‑1 (MQ=255)                                                       caGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATtctc                                                                                        <  2:105308/149‑1 (MQ=255)                                                       caGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATtctc                                                                                        <  6:357325/149‑1 (MQ=255)                                                                tCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTAt                                                                               <  5:352078/149‑1 (MQ=255)                                                                tCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGCTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTAt                                                                               <  4:185331/149‑1 (MQ=255)                                                                   cGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATaa                                                                                    >  7:144437/1‑141 (MQ=255)                                                                   cGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATaa                                                                                    <  8:144437/141‑1 (MQ=255)                                                                      gCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATg                                                                         >  3:233533/1‑149 (MQ=255)                                                                      gCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGAAGATAGTGCTCAt                                                                                                                 <  4:153135/109‑1 (MQ=255)                                                                      gCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGAAGATAGTGCTCAt                                                                                                                 >  3:153135/1‑109 (MQ=255)                                                                       cTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCATCCTTATGGGAATTCTCATAAGGTATCCCATgg                                                                        >  4:67463/1‑149 (MQ=255)                                                                             gcgcACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAacac                                                                  <  1:4289/149‑1 (MQ=255)                                                                             acgcACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAacac                                                                  >  8:219384/2‑149 (MQ=255)                                                                                 aCCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCAt                                                              <  1:71010/149‑1 (MQ=255)                                                                                  cccAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCAtt                                                             <  5:176859/149‑1 (MQ=255)                                                                                   ccAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCAttt                                                            <  7:38361/149‑1 (MQ=255)                                                                                             aCTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAatat                                                  <  2:185800/149‑1 (MQ=255)                                                                                               tgtgCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATgg                                                >  1:210668/1‑149 (MQ=255)                                                                                                gtgCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATgg                                                >  1:350099/1‑148 (MQ=255)                                                                                                 agCTCTGGTCTACATTGAGTCGCTCAGCGGTAACCTCCTGCCTCCCCATAAAGGCCACGTGTTCTTTGTAGATTGTGCTCATGGCACCCCTCCTTATGTGAATTCTCATATGTTATCCCATGGGAACACTCATTTGCAAGAATATGGGa                                              <  2:152981/148‑1 (MQ=255)                                                                                                   ctctGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAgc                                            <  7:227571/149‑1 (MQ=255)                                                                                                       ggTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGGTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACCCTCATTTGCAAGAATATGGGAGCGccc                                        >  4:355901/1‑149 (MQ=255)                                                                                                        gTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGcccc                                       >  5:14073/1‑149 (MQ=255)                                                                                                         tCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATCGGAATTCTCATAAGGTATCCCAAGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCt                                      <  6:282509/149‑1 (MQ=255)                                                                                                           aaCATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGa                                    >  4:213858/1‑149 (MQ=255)                                                                                                           aaCATTCAGTCGCTCAGCGCTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTATGTAGATAGTGCTCATGGCACCCCTCCTTATCGGAATTCTCATAAGGTATCCCATGTGAACACTCATTTGCAAGAATATGGGAGCGCCACTg                                     >  6:178621/1‑148 (MQ=255)                                                                                                            aCATTGAGTCGCTCAGCGGTAACCTCCTGCTTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGAt                                   <  6:22034/149‑1 (MQ=255)                                                                                                             cATTGAGTCGCTCAGCGGTAACCTCCTGCGTCACGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATa                                  >  1:245327/1‑149 (MQ=255)                                                                                                              aTTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATAc                                 >  2:28692/1‑149 (MQ=255)                                                                                                                  aGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTc                              <  6:340675/148‑1 (MQ=255)                                                                                                                  aGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCt                             <  6:14073/149‑1 (MQ=255)                                                                                                                   gTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCt                             >  6:325040/1‑148 (MQ=255)                                                                                                                    tCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTg                            <  1:48292/148‑1 (MQ=255)                                                                                                                       cTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGatat                        <  5:136019/149‑1 (MQ=255)                                                                                                                                 aCCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTaaa              <  6:102109/149‑1 (MQ=255)                                                                                                                                 aCCTCCTGCGTACCGAATTTCGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTaaa              <  5:362040/149‑1 (MQ=255)                                                                                                                                    tcctGCTTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATTAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGa                                    <  3:317890/124‑1 (MQ=255)                                                                                                                                         cGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTAAAGACCGAAt      <  2:350099/149‑1 (MQ=255)                                                                                                                                             ccGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCTTGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTAAAGACCGAATCaaa  <  4:233533/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GAGCCTCCTTCTTGAATTTGCCCTCAGTCAGGCGTTCAGCTCGTTTAGCCACCACAGGAGACATCCCGTGCTTACCGCGCACCCAGCCGGAGACTGTGCTCTGGTCAACATTGAGTCGCTCAGCGGTAACCTCCTGCGTCCCGAAATAGGCGACGAGTTCTTTGTAGATAGTGCTCATGGCACCCCTCCTTATGGGAATTCTCATAAGGTATCCCATGGGAACACTCATTTGCAAGAATATGGGAGCGCCCCTGATACTTCTGATATGGAACTTAAAGACCGAATCAAA  >  NC_002947/4411500‑4411788

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