Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,907,385	C→G	6.2%	G186A (GGT→GCT)	eno ←	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,907,385	0	C	G	6.2%	228.1 / 5.6	64	G186A (GGT→GCT)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (2/2);  total (30/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 = 9.28e-01

CAGCAGCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTT  >  NC_000913/2907242‑2907530                                                                                                                                                |                                                                                                                                                  cagcagCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAAccc                                                                                                                                                 <  2:338458/146‑1 (MQ=255) cagcagCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAAccc                                                                                                                                                 >  1:338458/1‑146 (MQ=255)      gcTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGAt                                                                                                                                         <  2:20676/149‑1 (MQ=255)        tttAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATgg                                                                                                                                       >  1:114851/1‑149 (MQ=255)               gCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCttt                                                                                                                                <  4:179024/149‑1 (MQ=255)                     gCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGt                                                                                                                          <  1:187833/149‑1 (MQ=255)                      cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                         >  1:410851/1‑149 (MQ=255)                      cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                         <  3:306515/149‑1 (MQ=255)                             aGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGcttctt                                                                                                                                 >  4:378182/1‑134 (MQ=255)                             aGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGcttctt                                                                                                                                 <  3:378182/134‑1 (MQ=255)                              gCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGcc                                                                                                                 >  3:267500/1‑149 (MQ=255)                                      ttCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCt                                                                                                         >  4:439529/1‑149 (MQ=255)                                         aGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGACCTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                      <  4:109096/149‑1 (MQ=255)                                                  aCCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAAtt                                                                                             >  3:284635/1‑149 (MQ=255)                                                     cAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCt                                                                                          >  4:282076/1‑149 (MQ=255)                                                      aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTca                                                                                                                              <  1:3984/112‑1 (MQ=255)                                                      aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTca                                                                                                                              >  2:3984/1‑112 (MQ=255)                                                      aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                         <  4:245561/149‑1 (MQ=255)                                                      aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                         <  3:277401/149‑1 (MQ=255)                                                      aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                         <  3:15901/149‑1 (MQ=255)                                                        gTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGa                                                                                       >  4:81401/1‑149 (MQ=255)                                                             gcgcATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCa                                                                                  >  2:465210/1‑149 (MQ=255)                                                               gcATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAAc                                                                                >  4:192538/1‑149 (MQ=255)                                                                cATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACg                                                                               <  1:94880/149‑1 (MQ=255)                                                                    gCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgt                                                                           >  2:58273/1‑149 (MQ=255)                                                                     ccACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                          >  3:168294/1‑149 (MQ=255)                                                                       aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGaa                                                                                               >  1:61393/1‑126 (MQ=255)                                                                       aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGaa                                                                                               <  2:61393/126‑1 (MQ=255)                                                                       aCCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTtgttgt                                                                        >  3:341801/1‑149 (MQ=255)                                                                       aCCTTCGGCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCGGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACTGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTtgttgt                                                                        >  4:402583/1‑149 (MQ=255)                                                                                    aCAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTcacc                                                           >  1:76413/1‑149 (MQ=255)                                                                                         agTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGt                                                      <  3:81401/149‑1 (MQ=255)                                                                                          gtgtTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGtt                                                     >  1:22795/1‑149 (MQ=255)                                                                                            gtTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTtga                                                   <  4:326435/149‑1 (MQ=255)                                                                                               cATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTtgatga                                                <  4:267500/149‑1 (MQ=255)                                                                                                aTGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTgatgat                                               >  2:270523/1‑149 (MQ=255)                                                                                                 tGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGatgatg                                              <  4:284635/149‑1 (MQ=255)                                                                                                  gCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGata                                                                                     >  1:255376/1‑109 (MQ=255)                                                                                                  gCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGata                                                                                     <  2:255376/109‑1 (MQ=255)                                                                                                      ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                         <  3:282076/149‑1 (MQ=255)                                                                                                      ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                         <  2:438341/149‑1 (MQ=255)                                                                                                        cGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcatca                                       >  3:186803/1‑149 (MQ=255)                                                                                                                gAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAAcc                               >  3:427152/1‑149 (MQ=255)                                                                                                                   cTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGc                            >  1:91801/1‑149 (MQ=255)                                                                                                                      ttGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATa                         <  3:439529/149‑1 (MQ=255)                                                                                                                        gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       >  4:439442/1‑149 (MQ=255)                                                                                                                        gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       <  1:58273/149‑1 (MQ=255)                                                                                                                        gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       <  3:250259/149‑1 (MQ=255)                                                                                                                        gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       <  2:76413/149‑1 (MQ=255)                                                                                                                           aGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACCTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTa                    >  2:191660/1‑149 (MQ=255)                                                                                                                            ggTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAt                   <  2:114851/149‑1 (MQ=255)                                                                                                                              tGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTAttt                 <  2:22795/149‑1 (MQ=255)                                                                                                                                 tGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                      >  2:272483/1‑61 (MQ=255)                                                                                                                                 tGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGaa                                                                                                      <  1:272483/61‑1 (MQ=255)                                                                                                                                   gAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCg            <  2:410851/149‑1 (MQ=255)                                                                                                                                     aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa          >  1:292971/1‑149 (MQ=255)                                                                                                                                     aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa          <  3:402583/149‑1 (MQ=255)                                                                                                                                     aaaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGa          <  1:270523/149‑1 (MQ=255)                                                                                                                                      aaCTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAg         <  2:292971/149‑1 (MQ=255)                                                                                                                                        cTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTa       <  4:168294/149‑1 (MQ=255)                                                                                                                                          tCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTAcc     <  1:465210/149‑1 (MQ=255)                                                                                                                                          aaaGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                          >  1:273792/3‑80 (MQ=37)                                                                                                                                          aaaGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                          <  2:273792/78‑1 (MQ=37)                                                                                                                                            aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       >  4:386668/1‑129 (MQ=255)                                                                                                                                            aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                       <  3:386668/129‑1 (MQ=255)                                                                                                                                             gAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgtt  >  4:11296/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                  CAGCAGCTTTAACAGCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTT  >  NC_000913/2907242‑2907530

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