Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,426,692	A→G	9.2%	D109G (GAC→GGC)	PP_3920 →	putative Phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,426,692	0	A	G	9.2%	174.8 / 7.2	65	D109G (GAC→GGC)	PP_3920	putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (3/3);  total (33/32)
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.76e-01

AAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGCGT  >  NC_002947/4426544‑4426839                                                                                                                                                     |                                                                                                                                                    aaGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGa                                                                                                                                                     >  1:337253/1‑149 (MQ=255)        tGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGTCGATCGGCGccgag                                                                                                                                              <  1:39662/149‑4 (MQ=255)        tGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATc                                                                                                                                              >  5:282826/1‑149 (MQ=255)        tGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATc                                                                                                                                              <  1:850/149‑1 (MQ=255)               ggAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCgg                                                                                                                                       <  6:36366/149‑1 (MQ=255)                 aaCTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCAAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGcc                                                                                                                                     >  6:47861/1‑149 (MQ=255)                   cTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCAt                                                                                                                                   >  5:57354/1‑149 (MQ=255)                   cTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCAt                                                                                                                                   <  7:196307/149‑1 (MQ=255)                         aCGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATcgc                                                                                                                             <  2:337253/149‑1 (MQ=255)                                 gggAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTAcc                                                                                                                     >  2:142572/1‑149 (MQ=255)                                            acgaGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCgg                                                                                                          >  5:101121/1‑149 (MQ=255)                                              gaGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTc                                                                                                        >  3:262723/1‑149 (MQ=255)                                                  aCCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGc                                                                                                     >  2:75654/1‑148 (MQ=255)                                                    ctgctgAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGaa                                                                                                  <  2:244735/149‑1 (MQ=255)                                                    ctgctgAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGaa                                                                                                  <  8:223306/149‑1 (MQ=255)                                                         gAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGcc                                                                                                                                     >  2:258013/1‑109 (MQ=255)                                                         gAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGcc                                                                                                                                     <  1:258013/109‑1 (MQ=255)                                                          aaGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATcgcg                                                                                                                            <  4:51782/117‑1 (MQ=255)                                                          aaGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATcgcg                                                                                                                            >  3:51782/1‑117 (MQ=255)                                                          aaGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGAcc                                                                                              <  4:310316/147‑1 (MQ=255)                                                             aTCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGc                                                                                         <  6:327031/149‑1 (MQ=255)                                                                 tCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCaaa                                                                                      >  4:310770/1‑148 (MQ=255)                                                                  cccAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGAGCAATGTGTTCTACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCaaa                                                                                      >  6:270854/1‑147 (MQ=255)                                                                     aGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAgg                                                                                  >  7:291419/1‑148 (MQ=255)                                                                       gCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGa                                                                                >  8:246953/1‑148 (MQ=255)                                                                         aaGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAgg                                                                                           >  1:148036/1‑135 (MQ=255)                                                                         aaGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAAt                                                                             >  8:212930/1‑149 (MQ=255)                                                                          aGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATc                                                                            <  1:142572/149‑1 (MQ=255)                                                                          aGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATc                                                                            >  8:214386/1‑149 (MQ=255)                                                                               gTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGTCGCCGATCGACGCAATCACCCCGTCCATGATCGCGCAGTACAGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGTAAATCGCCCt                                                                       <  6:127613/149‑1 (MQ=255)                                                                                tCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTg                                                                      <  6:101121/149‑1 (MQ=255)                                                                                       ccAGAAGGACAACATCTACGAGCTCAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTccc                                                               <  5:101182/149‑1 (MQ=255)                                                                                         aGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCAc                                                             <  1:75654/149‑1 (MQ=255)                                                                                         aGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCAc                                                             <  6:282826/149‑1 (MQ=255)                                                                                              gACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTAtt                                                        >  8:238959/1‑149 (MQ=255)                                                                                                caacaTCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTAATc                                                       >  2:16167/1‑148 (MQ=255)                                                                                                 aacaTCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCaa                                                     >  2:186058/1‑149 (MQ=255)                                                                                                  acaTCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAac                                                    <  6:81434/149‑1 (MQ=255)                                                                                                        aCGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCg                                              >  3:214587/1‑149 (MQ=255)                                                                                                          gAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAg                                            >  5:306494/1‑149 (MQ=255)                                                                                                           aGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAgg                                           <  6:57354/149‑1 (MQ=255)                                                                                                             cTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAggg                                          >  3:323179/1‑148 (MQ=255)                                                                                                                 agcagcTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATgg                                     <  7:246953/149‑1 (MQ=255)                                                                                                                  gcagcTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATggg                                    <  8:291419/149‑1 (MQ=255)                                                                                                                      cTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCt                                >  3:265966/1‑149 (MQ=255)                                                                                                                          gcAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTAc                             >  6:38337/1‑148 (MQ=255)                                                                                                                           cAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCAc                                                             >  5:24673/1‑115 (MQ=255)                                                                                                                           cAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCAc                                                             <  6:24673/115‑1 (MQ=255)                                                                                                                                tgtTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCgg                      <  7:212930/149‑1 (MQ=255)                                                                                                                                tgtTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCg                       <  7:214386/148‑1 (MQ=255)                                                                                                                                     gACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACaa                 >  7:30418/1‑149 (MQ=255)                                                                                                                                       gatggCGTCGATCGGCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCt                                >  8:280070/4‑132 (MQ=255)                                                                                                                                       gatggCGTCGATCGGCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCt                                <  7:280070/129‑1 (MQ=255)                                                                                                                                       gatggCGTCGATCGGCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAAcc               >  2:12982/4‑149 (MQ=255)                                                                                                                                        atggCGTCGATCGGCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAAccc              >  8:263719/3‑149 (MQ=255)                                                                                                                                             gcCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACACGGAAATCGCCCTGCTTTCCCCCGTAGTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGGGcc         >  2:242552/1‑149 (MQ=255)                                                                                                                                              ccGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCt        >  3:157179/1‑149 (MQ=255)                                                                                                                                                gATCGACTCCATCACCCCGGACATGATCTCGCAGTACCGCGACTAGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      <  5:270854/149‑1 (MQ=255)                                                                                                                                                gATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      >  8:184722/1‑149 (MQ=255)                                                                                                                                                gATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      >  8:186138/1‑149 (MQ=255)                                                                                                                                                gATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      <  7:184722/149‑1 (MQ=255)                                                                                                                                                gATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      <  3:310770/149‑1 (MQ=255)                                                                                                                                                gATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTgg      <  4:157179/149‑1 (MQ=255)                                                                                                                                                  tCGGCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGc    <  7:263719/149‑1 (MQ=255)                                                                                                                                                  tCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGc    <  6:306494/149‑1 (MQ=255)                                                                                                                                                  tCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGc    <  4:214587/149‑1 (MQ=255)                                                                                                                                                     aCGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTAAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGCGt  >  1:304544/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCTGAAGATCGTCCCAGGCAAGGCCGTCAGGACCCAGAAGGACAACATCTACGAGCTGAAGCAGCTGCGCAATGTGTTCGACTCGGCGCCGATCGACGCCATCACCCCGGCCATGATCGCGCAGTACCGCGACTCGCGGTCGGCGAAGACCAGGGCAAACAGGGAAATCGCCCTGCTTTCCCACGTATTCAACACGGCGAGGGAATGGGGGCTTACAACCCGGGACAACCCGTGCCTGGGCGT  >  NC_002947/4426544‑4426839

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