Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,404,872	T→C	pseudogene (141/300 nt)	yrdE →	pseudogene, yedZ paralog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,404,872	0	T	C	100.0%	268.5 / NA	68	pseudogene (141/300 nt)	yrdE	pseudogene, yedZ paralog
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/37);  total (31/37)

CTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAA  >  NC_000913/3404731‑3405015                                                                                                                                              |                                                                                                                                                cTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATcccc                                                                                                                                          >  2:312959/1‑149 (MQ=255)   ggTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCa                                                                                                                                         <  2:670141/148‑1 (MQ=255)      aCCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAAcc                                                                                                                                   >  1:959725/1‑151 (MQ=255)       ccATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCg                                                                                                                                  >  2:397362/1‑151 (MQ=255)        cATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGc                                                                                                                                 <  2:877232/151‑1 (MQ=255)              gTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGt                                                                                                                            >  1:10904/1‑150 (MQ=255)               tGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTa                                                                                                                           >  1:564279/1‑150 (MQ=255)                gAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTaa                                                                                                                          <  2:69548/150‑1 (MQ=255)                gAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTaa                                                                                                                          <  1:823138/150‑1 (MQ=255)                gAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTaa                                                                                                                          <  2:56689/150‑1 (MQ=255)                gAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAAt                                                                                                                         <  1:555186/151‑1 (MQ=255)                  ataatCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATgt                                                                                                                       >  1:761345/1‑151 (MQ=255)                   taatCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATgtc                                                                                                                      >  2:519052/1‑151 (MQ=255)                      tCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATgtcgtc                                                                                                                   >  2:109120/1‑151 (MQ=255)                         gATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTAt                                                                                                                <  2:213873/151‑1 (MQ=255)                             ttgtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                                                                             <  2:3859/150‑1 (MQ=255)                             ttgtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                                                                             <  1:184351/150‑1 (MQ=255)                             ttgtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCtt                                                                                                             <  2:303255/150‑1 (MQ=255)                               gtttgACTACAAATTAATCACTCGAGCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGt                                                                                                           <  2:793593/150‑1 (MQ=255)                               gtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGt                                                                                                           <  1:544894/150‑1 (MQ=255)                               gtttgACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTc                                                                                                          >  2:222182/1‑151 (MQ=255)                                   gACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTcgcg                                                                                                       <  2:953195/150‑1 (MQ=255)                                    aCTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGAt                                                                                                     <  1:395405/151‑1 (MQ=255)                                      tACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATcc                                                                                                   >  2:437564/1‑151 (MQ=255)                                        cAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTg                                                                                                 <  2:565961/151‑1 (MQ=255)                                         aaaTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTgg                                                                                                >  2:679705/1‑151 (MQ=255)                                           aTTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGc                                                                                               >  1:806249/1‑150 (MQ=255)                                             tAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATc                                                                                            >  1:161906/1‑151 (MQ=255)                                              aaTCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATcc                                                                                           <  2:264541/151‑1 (MQ=255)                                               aTCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATccc                                                                                          >  2:29784/1‑151 (MQ=255)                                                 cACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCt                                                                                         <  2:633320/150‑1 (MQ=255)                                                     cGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatt                                                                                    >  2:147266/1‑151 (MQ=255)                                                      gAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                                                                                   <  2:15644/151‑1 (MQ=255)                                                      gAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatta                                                                                   <  2:376698/151‑1 (MQ=255)                                                      gAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACatt                                                                                    >  2:64654/1‑150 (MQ=255)                                                       aaCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattat                                                                                  >  2:83296/1‑151 (MQ=255)                                                        aCCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACattatt                                                                                 >  1:636286/1‑151 (MQ=255)                                                          cTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTTGTCTATCTTGTCGCGATCCTGGCATCACTACattatt                                                                                 >  2:436971/1‑149 (MQ=255)                                                          cTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                                                                                >  2:96523/1‑150 (MQ=255)                                                          cTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTAttt                                                                                <  2:269526/150‑1 (MQ=255)                                                            aTTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtg                                                                               <  1:523280/149‑1 (MQ=255)                                                            aTTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTtg                                                                               <  2:135221/149‑1 (MQ=255)                                                              ttAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTgg                                                                            <  2:560490/150‑1 (MQ=255)                                                                 aTGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCtg                                                                        <  1:784349/151‑1 (MQ=255)                                                                     tGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGa                                                                     >  2:224642/1‑150 (MQ=255)                                                                       aGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGa                                                                  <  2:775871/151‑1 (MQ=255)                                                                                   cGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTg                                                       <  2:44018/150‑1 (MQ=255)                                                                                             tGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCa                                             <  2:564221/150‑1 (MQ=255)                                                                                               cGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCt                                          <  1:271509/151‑1 (MQ=255)                                                                                               cGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCt                                          >  1:325786/1‑151 (MQ=255)                                                                                                gTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCt                                          <  1:31851/150‑1 (MQ=255)                                                                                                  gCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTc                                        <  2:769760/150‑1 (MQ=255)                                                                                                  gCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTc                                        <  2:105215/150‑1 (MQ=255)                                                                                                     ttAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGc                                      <  1:853853/149‑1 (MQ=255)                                                                                                          aTTCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctg                               <  1:303609/151‑1 (MQ=255)                                                                                                           ttCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctgg                              >  2:506386/1‑151 (MQ=255)                                                                                                            tCACATCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGggctggc                             >  1:241804/1‑151 (MQ=255)                                                                                                              acaTCTATCCAGACGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGgctt                           >  1:53018/1‑151 (MQ=255)                                                                                                                           cGATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCAc              >  2:734701/1‑151 (MQ=255)                                                                                                                            gATGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCAcc             >  1:732470/1‑151 (MQ=255)                                                                                                                              tGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt           >  1:718857/1‑151 (MQ=255)                                                                                                                              tGCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt           >  1:857854/1‑151 (MQ=255)                                                                                                                               gCAGTGAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCtt           >  1:662592/1‑150 (MQ=255)                                                                                                                                    gAAAATTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTta     <  1:312996/151‑1 (MQ=255)                                                                                                                                     aaaaTTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat    <  1:621647/151‑1 (MQ=255)                                                                                                                                     aaaaTTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat    <  1:29788/151‑1 (MQ=255)                                                                                                                                     aaaaTTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTtat    >  1:784679/1‑151 (MQ=255)                                                                                                                                       aaTTGGGCAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATaa  <  2:77587/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                CTGGTACCATGGGGTGAATAATCTGATTTTGTTTGACTACAAATTAATCACTCGAACCTATTTAATGCTGAGCATTGTCAATCGGTTAATTTTGCGTGCTTTAGCATTCACATCTATCCAGACGATGCAGTGAAAATTGGGTAATCCCCAGCAACCGCTGCGTAATGTCGTCTATCTTGTCGCGATCCTGGCATCCCTACATTATTTGTGGTCTGTGAAGATTATCTCATTGCAGCCCCTCATCTTCGCAGGGCTGGCTTTTCAGCTTTTCACCTTACGTTATAA  >  NC_000913/3404731‑3405015

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