Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,699,460	C→T	6.0%	intergenic (+150/‑119)	PP_5002 → / → phaA	hypothetical protein/poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,699,460	0	C	T	6.0%	196.6 / 2.3	67	intergenic (+150/‑119)	PP_5002/phaA	hypothetical protein/poly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base T (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GGCCTGCGGGGTTTAGAGCGCATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACGAT  >  NC_002947/5699316‑5699599                                                                                                                                                 |                                                                                                                                            ggCCTGCGGGGTTTAGAGCGCATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAAc                                                                                                                                        <  2:1504644/150‑1 (MQ=255)  gCCTGCGGGGTTTAGAGCGCATTTTCTAAATTCATCTGTTTGAATGACCTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAAcc                                                                                                                                       <  1:1214597/150‑1 (MQ=255)    cTGCGGGGTTTAGAGCGCATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCCAGTCGGCAGAACCCa                                                                                                                                     >  2:617306/1‑150 (MQ=255)                   cgcATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGtt                                                                                                                      <  2:568868/150‑1 (MQ=255)                   cgcATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGtt                                                                                                                      >  2:182697/1‑150 (MQ=255)                     cATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTcc                                                                                                                    <  1:357185/150‑1 (MQ=255)                       ttttCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCtt                                                                                                                  >  1:632487/1‑150 (MQ=255)                       ttttCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGAGATGTCTCGCGCATTACACGAAAGTCGGGTAACTAATGGGTGTGGCAAGTTTCCTGCATGACTTTGCAGGCCGGCACAACCCACGAAGCACCGCTGTCCCtt                                                                                                                  <  2:28533/150‑1 (MQ=255)                          tCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATc                                                                                                               <  1:182697/150‑1 (MQ=255)                            tAAATTCATCTGTTTGAATGACTTGAAAACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCAc                                                                                                             >  1:863266/1‑150 (MQ=255)                                       gTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTGCCCTGCATGACTTTGCAGGTCGGCAGACCCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGa                                                                                                  <  2:145796/150‑1 (MQ=255)                                        tttGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGt                                                                                                         >  2:764275/1‑142 (MQ=255)                                        tttGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGt                                                                                                         <  1:764275/142‑1 (MQ=255)                                        tttGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAg                                                                                                 >  1:407373/1‑150 (MQ=255)                                          tGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTa                                                                                               <  2:632487/150‑1 (MQ=255)                                            aaTGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTAGACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTGCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGc                                                                                             <  2:1402812/150‑1 (MQ=255)                                                  ttGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCTGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTAcc                                                                                       <  1:617306/150‑1 (MQ=255)                                                      aGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGc                                                                                   >  2:694895/1‑150 (MQ=255)                                                       gACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAAGCCGCAGAGTACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCt                                                                                  >  1:918814/1‑150 (MQ=255)                                                       gACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGt                                                                                          <  1:955295/142‑1 (MQ=255)                                                       gACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGt                                                                                          >  2:955295/1‑142 (MQ=255)                                                       gACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCt                                                                                  >  1:615151/1‑150 (MQ=255)                                                            cccAGTGAAGGCCTGGCTTCCGCATTCCACGAAACTCGGCTAACCAATGGGTGTGGCCAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTAGCACTGGTCACCCGAGTAGCAGTACCGGGCTAAgga                                                                             <  1:744756/150‑2 (MQ=255)                                                                  gAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGc                                                                       <  2:615151/150‑1 (MQ=255)                                                                                       cacGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCa                                                                                                              >  1:702316/1‑90 (MQ=255)                                                                                       cacGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCa                                                                                                              <  2:702316/90‑1 (MQ=255)                                                                                        acGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCAccc                                                                                                    >  1:847184/1‑99 (MQ=255)                                                                                        acGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCAccc                                                                                                    <  2:847184/99‑1 (MQ=255)                                                                                        acGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAggg                                                                             >  2:1084447/1‑122 (MQ=255)                                                                                        acGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAggg                                                                             <  1:1084447/122‑1 (MQ=255)                                                                                        acGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGAATCAACCCCCTGAACAAATCACTCTACCCGCCACAGCAACCGGTACt                                                 >  2:15237/1‑150 (MQ=255)                                                                                                   gTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGAc                                      >  2:957804/1‑150 (MQ=255)                                                                                                       ccAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGTAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCtgt                                                                         <  1:663403/111‑1 (MQ=255)                                                                                                       ccAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGTAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCtgt                                                                         >  2:663403/1‑111 (MQ=255)                                                                                                          aTGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAg                               <  2:407373/150‑1 (MQ=255)                                                                                                              gtgtGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTAt                                                                                                                >  2:1326155/1‑65 (MQ=255)                                                                                                              gtgtGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTAt                                                                                                                <  1:1326155/65‑1 (MQ=255)                                                                                                                   gCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGa                      <  1:957804/150‑1 (MQ=255)                                                                                                                   gCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGGCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAAGACAACGGAGCGTCGTAGa                      >  2:19361/1‑150 (MQ=255)                                                                                                                    aaaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGAt                     <  2:411108/149‑1 (MQ=255)                                                                                                                     aaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAggg                                                                             >  1:50776/1‑93 (MQ=255)                                                                                                                     aaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAggg                                                                             <  2:50776/93‑1 (MQ=255)                                                                                                                     aaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTa                        >  1:797559/1‑146 (MQ=255)                                                                                                                     aaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTa                        <  2:797559/146‑1 (MQ=255)                                                                                                                     aaGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATg                    >  2:1595699/1‑150 (MQ=255)                                                                                                                        ttCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAg                               <  2:1482397/136‑1 (MQ=255)                                                                                                                        ttCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAg                               >  1:1482397/1‑136 (MQ=255)                                                                                                                        ttCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGt                 >  1:1535361/1‑150 (MQ=255)                                                                                                                        cgaccTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGt                                              >  1:1062223/4‑121 (MQ=255)                                                                                                                        cgaccTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGt                                              <  2:1062223/118‑1 (MQ=255)                                                                                                                           ccTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAAc                                                        >  2:1413764/1‑108 (MQ=255)                                                                                                                           ccTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAAc                                                        <  1:1413764/108‑1 (MQ=255)                                                                                                                               cATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTa                                                                                               >  2:1607286/1‑65 (MQ=255)                                                                                                                               cATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTa                                                                                               <  1:1607286/65‑1 (MQ=255)                                                                                                                                    cTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGaacaac     >  1:1079461/1‑150 (MQ=255)                                                                                                                                     tttGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACg    >  2:666539/1‑150 (MQ=255)                                                                                                                                     tttGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACg    >  2:444001/1‑150 (MQ=255)                                                                                                                                      ttGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTAACACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACGa   <  1:15237/150‑1 (MQ=255)                                                                                                                                       tGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACGAt  <  2:863266/150‑1 (MQ=255)                                                                                                                                        gCAGGTCGGTAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTa                        >  2:572863/1‑127 (MQ=255)                                                                                                                                        gCAGGTCGGTAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTa                        <  1:572863/127‑1 (MQ=255)                                                                                                                                          aGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAg                                                                               <  1:1473742/70‑1 (MQ=255)                                                                                                                                          aGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAg                                                                               >  2:1473742/1‑70 (MQ=255)                                                                                                                                          aGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGaa         <  1:1495753/140‑1 (MQ=255)                                                                                                                                          aGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGaa         >  2:1495753/1‑140 (MQ=255)                                                                                                                                              cGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAg                       <  1:299227/122‑1 (MQ=255)                                                                                                                                              cGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAg                       >  2:299227/1‑122 (MQ=255)                                                                                                                                                 |                                                                                                                                            GGCCTGCGGGGTTTAGAGCGCATTTTCTAAATTCATCTGTTTGAATGACTTGCAGACAGCCCAGTGAAGGGCTGTCTTGCGCATTACACGAAAGTCGGGTAACCAATGGGTGTGGCAAGTTCCCTGCATGACTTTGCAGGTCGGCAGAACCCACGCAGCACCGCTGTTCCTTATCACTGGTCACCCGAGTAGCAGTACCGGGCTCAGGGCTGTGCACCCGCCACAGCAACCGGTACTCGTCTCAGGACAACGGAGCGTCGTAGATGAGTAACAAGAACAACGAT  >  NC_002947/5699316‑5699599

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