Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,907,398	C→G	6.1%	A182P (GCC→CCC)	eno ←	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,907,398	0	C	G	6.1%	234.1 / 6.2	66	A182P (GCC→CCC)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

GCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGC  >  NC_000913/2907256‑2907546                                                                                                                                               |                                                                                                                                                     gCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCttt                                                                                                                                                <  4:179024/149‑1 (MQ=255)       gCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGt                                                                                                                                          <  1:187833/149‑1 (MQ=255)        cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                                         <  3:306515/149‑1 (MQ=255)        cGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGtt                                                                                                                                         >  1:410851/1‑149 (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                                                                                                                                              >  2:3984/1‑112 (MQ=255)                                        aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTca                                                                                                                                              <  1:3984/112‑1 (MQ=255)                                        aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                                         <  3:277401/149‑1 (MQ=255)                                        aGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTg                                                                                                         <  4:245561/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                                                                                                     <  2:255376/109‑1 (MQ=255)                                                                                    gCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGata                                                                                                     >  1:255376/1‑109 (MQ=255)                                                                                        ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                                         <  2:438341/149‑1 (MQ=255)                                                                                        ttCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTtcat                                                         <  3:282076/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                                       <  1:58273/149‑1 (MQ=255)                                                                                                          gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                       <  3:250259/149‑1 (MQ=255)                                                                                                          gCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                       >  4:439442/1‑149 (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                                                                                          <  2:273792/78‑1 (MQ=37)                                                                                                                            aaaGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACgttgtt                                                                                          >  1:273792/3‑80 (MQ=37)                                                                                                                              aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                       <  3:386668/129‑1 (MQ=255)                                                                                                                              aGAAGCCATCCGCATGGGTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATaga                                       >  4:386668/1‑129 (MQ=255)                                                                                                                               gAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgtt                  >  4:11296/1‑149 (MQ=255)                                                                                                                                   ccATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAACTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagt              >  1:225716/1‑149 (MQ=255)                                                                                                                                     aTGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACGGGCTGAATCTTGAATTCCTGGATATCAACGTTGTTGCCGGCTTGCTCACCACCGTGGATGATGTCCACCATCGGAACCGGCATAGAGTTTTTGCCCGGAGTACCgttcagttc            >  2:405608/1‑149 (MQ=255)                                                                                                                                     aTGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCgttcagttc            >  3:128501/1‑149 (MQ=255)                                                                                                                                        cGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGc         <  1:345528/149‑1 (MQ=255)                                                                                                                                           aTGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGAt      <  3:27927/149‑1 (MQ=255)                                                                                                                                             ggCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGAtgt    <  2:91801/149‑1 (MQ=255)                                                                                                                                               cttcttTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGc  >  1:370185/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                     GCTTCAGCGATAACAGCCAGAGCTTCAGCGTTGGAACCCAGGTTCGGCGCATAGCCACCTTCGTCACCAACAGCAGTGTTCATGCCTTTCGCTTTCAGAACTTTTGCCAGGTGATGGAAAACTTCAGAACCCATGCGGATGGCTTCTTTCACAGTTTTCGCGCCAACCGGCTGAATCATGAATTCCTGGATATCAACGTTGTTGTCAGCGTGCTCACCACCGTTGATGATGTTCATCATCGGAACCGGCATAGAGTATTTGCCCGGAGTACCGTTCAGTTCAGCGATGTGC  >  NC_000913/2907256‑2907546

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