Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	4,015,535	Δ1 bp	coding (635/816 nt)	ysgA ←	putative carboxymethylenebutenolidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	4,015,534	0	A	.	96.8%	265.4 / 1.9	63	coding (636/816 nt)	ysgA	putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (33/27);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGCGGCT  >  NC_000913_3_pae_pgi/4015389‑4015678                                                                                                                                                  |                                                                                                                                                 aTTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCacgc                                                                                                                                               >  2:113687/1‑146 (MQ=255)  ttCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCacgct                                                                                                                                              >  2:302143/1‑145 (MQ=255)     cAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCtctc                                                                                                                                          >  2:69373/1‑149 (MQ=255)          tatgCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCg                                                                                                                                     <  1:263465/146‑1 (MQ=255)                       aTCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATc                                                                                                                        <  1:32073/149‑1 (MQ=255)                             tGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGCGTAAATGTTAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGacc                                                                                                                  >  1:910978/1‑149 (MQ=255)                              gCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCAGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGacca                                                                                                                 >  2:573580/1‑149 (MQ=255)                                        catgcatgCTCGGGCGATAATCAGCCGTGAAGGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGTCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGc                                                                                                       <  2:320412/142‑1 (MQ=255)                                             tAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCgaga                                                                                                  >  1:477693/1‑149 (MQ=255)                                                 tCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATc                                                                                                                        <  1:374078/123‑1 (MQ=255)                                                 tCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATc                                                                                                                        >  2:374078/1‑123 (MQ=255)                                                     gCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACGACCGTATAAGCCGAGAATCGgcgc                                                                                          >  1:736310/1‑149 (MQ=255)                                                     gCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGgcgc                                                                                          <  1:113687/149‑1 (MQ=255)                                                       gATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGt                                                                                        >  2:871633/1‑149 (MQ=255)                                                        atcttCAGCGTTGAATACATGCCCGGCGTCCGGGTACGCGGTAATCTCTGCTTTCGCATCAGCTGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCTTATAAGCCGAGAATCGGCGCttt                                                                                       <  2:893496/145‑1 (MQ=255)                                                        aTAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGtt                                                                                       >  1:623237/1‑149 (MQ=255)                                                         tAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTACGGAATGCTGTTATCCTGACCACCATAGAAGCCGAGAATTGgctctgta                                                                                      >  2:321902/1‑143 (MQ=255)                                                             cAGCGTTGAATGCATGCCAGGCGACCGGGTACACGGTAATCTCTGCTTTCGCATTCGCAGCCCGCAGCGCCTGGCGCAGGGTTTC‑ACGCTCTCCTGAGGAATGCTGTTATCCTGACCACCATATGAGCCGACAATCGGCGCGTTAAGAt                                                                                  >  1:893496/1‑149 (MQ=255)                                                              aGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATc                                                                                 >  2:809022/1‑149 (MQ=255)                                                               gCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTg                                                                                                                             <  2:522992/104‑1 (MQ=255)                                                               gCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTg                                                                                                                             >  1:522992/1‑104 (MQ=255)                                                                   tGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGTGCGTTAAGATCGGTTg                                                                            <  1:488402/149‑1 (MQ=255)                                                                          atgcCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                                     <  2:175457/133‑1 (MQ=255)                                                                          atgcCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTaa                                                                                     >  1:175457/1‑133 (MQ=255)                                                                            gcCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGCATCGGCGCGTTAAGATCGGTTGCGATATCaa                                                                   >  1:365275/1‑149 (MQ=255)                                                                             cccGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAAc                                                                  <  1:378046/149‑1 (MQ=255)                                                                             cccGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAAc                                                                  <  2:602681/149‑1 (MQ=255)                                                                               cGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAg                                                                <  1:867316/149‑1 (MQ=255)                                                                                  cggccGGGGACACGACAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                             <  2:418194/146‑1 (MQ=255)                                                                                  cGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgtt                                                             <  1:573580/149‑1 (MQ=255)                                                                                    tCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGgttgt                                                           <  1:113742/149‑1 (MQ=255)                                                                                         gTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATCTAAGCCGCGAATCGGCGCGTTAAGATCGGGTGCACTATCAACAGGTTGTTTCgg                                                      >  2:807424/1‑149 (MQ=255)                                                                                         gTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGCTCGGTTGCGATATCGACAGGTTGtctccg                                                      >  1:126423/1‑147 (MQ=255)                                                                                          tACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                                     <  2:949437/149‑1 (MQ=255)                                                                                          tACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                                     <  1:302143/149‑1 (MQ=255)                                                                                          tACACGATAATCGCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGt                                                     <  2:3678/149‑1 (MQ=255)                                                                                            cacGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGa                                                   >  1:436135/1‑149 (MQ=255)                                                                                              ccacaacctctGCTTTCGCATTAGCCGCCCGCAGCGCCTGGCGCATGG‑TTCCCCGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAt                                                 <  1:362721/142‑1 (MQ=255)                                                                                               gATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                                >  2:104675/1‑149 (MQ=255)                                                                                                 tAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAGCAGGTTGTTTCGGTGAATTCa                                              >  1:196813/1‑149 (MQ=255)                                                                                                   ctccctGCTTTCCTATTAGCAGCCCGCAGCGCCGGGCGTATGGTTTC‑CCGCTCGCCGGCGGAATGGTGTTATCCTTACCACCATATAATCCGATAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGc                                            <  2:471275/145‑1 (MQ=255)                                                                                                    tctcTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                           >  1:849360/1‑149 (MQ=255)                                                                                                    tctcTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCg                                           >  1:196565/1‑149 (MQ=255)                                                                                                    tctcTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCACCg                                           >  1:121372/1‑149 (MQ=255)                                                                                                      tctGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGAc                                         >  2:110561/1‑149 (MQ=255)                                                                                                           tttCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                                <  1:632572/137‑1 (MQ=255)                                                                                                           tttCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAAtt                                                >  2:632572/1‑137 (MQ=255)                                                                                                           tttCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTc                                    >  1:169735/1‑149 (MQ=255)                                                                                                            ttCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCg                                   <  1:227976/149‑1 (MQ=255)                                                                                                               gCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATCAGCCGAAAATCGGCGCGTTAAGATCGGTTGTGATATCAACAGTTTGTTTCGGGGAATTCAGCTACTTGTCGccc                                >  2:330474/1‑149 (MQ=255)                                                                                                                cATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCg                               >  1:245557/1‑149 (MQ=255)                                                                                                                 aTTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                             <  1:141447/134‑1 (MQ=255)                                                                                                                 aTTAGCAGCCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAg                                             >  2:141447/1‑134 (MQ=255)                                                                                                                   tAGCAGCCCGCAGCGCCTGGCCCATGGTCGC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCa                            <  1:563595/149‑1 (MQ=255)                                                                                                                      cagcCCGCAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCTGCGACTTGTCGCCCGTCAGtt                         >  2:379196/1‑149 (MQ=255)                                                                                                                        tcccGCAGCGCCTGTCGCATGGATTC‑ACGCTCACCTGCGGAATGCTGTTATCCTTACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTTAATTCAGCGACTTGTCGCCCGTCAGTTTg                       <  1:883704/148‑1 (MQ=255)                                                                                                                             cAGCGCCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTa                  <  2:436135/149‑1 (MQ=255)                                                                                                                              aGCGCCTGGCGCATGGTTTC‑ACGCTCCCCTGCGGAATGCTGTGACCCGCACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACGTGTCGCCCGTCAGTTTGCCGTAc                 <  2:546704/149‑1 (MQ=255)                                                                                                                                 gcCTGGCGCATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCAc              <  1:361732/149‑1 (MQ=255)                                                                                                                                   cTGGCGGATGGTTT‑AACGCACTCCAGAGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGc            <  1:267768/149‑1 (MQ=255)                                                                                                                                         cATGGTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCCCTGc      >  2:825395/1‑149 (MQ=255)                                                                                                                                          atgcattCAACGGTGTCCTGCGGAACGCTGTTATCCTGACCACCAGATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGc      <  2:947330/144‑1 (MQ=255)                                                                                                                                            ggTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGCGGc   >  1:437872/1‑149 (MQ=255)                                                                                                                                            ggTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGCGGc   >  1:459616/1‑149 (MQ=255)                                                                                                                                             gTTTC‑ACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGCGGCt  >  2:680984/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 ATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAAGCCGAGAATCGGCGCGTTAAGATCGGTTGCGATATCAACAGGTTGTTTCGGTGAATTCAGCGACTTGTCGCCCGTCAGTTTGCCGTACCACGCCACTGCGGCT  >  NC_000913_3_pae_pgi/4015389‑4015678

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