Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,082,696	G→C	14.7%	Q22E (CAG→GAG)	fdoH ←	formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,082,696	0	G	C	14.7%	223.2 / 31.8	75	Q22E (CAG→GAG)	fdoH	formate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base C (7/4);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

GGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATG  >  NC_000913/4082554‑4082842                                                                                                                                               |                                                                                                                                                   gggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                                                                                                              <  4:413208/149‑1 (MQ=255) gggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGgcgc                                                                                                                                              >  3:759231/1‑149 (MQ=255)   gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                            <  3:768227/149‑1 (MQ=255)   gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg                                                                                                                                            <  4:759231/149‑1 (MQ=255)    ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                           <  2:755680/149‑1 (MQ=255)         gTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                                                                                                      <  4:572543/149‑1 (MQ=255)           acacCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCAACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                                                                                                    >  1:662501/1‑149 (MQ=255)               cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                >  4:544136/1‑149 (MQ=255)               cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                >  2:755061/1‑149 (MQ=255)               cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                <  1:755061/149‑1 (MQ=255)                 cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCTCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                                                                                                 >  3:754071/1‑146 (MQ=255)                 cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                                                                                                 <  4:754071/146‑1 (MQ=255)                 cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                                                                              >  4:311039/1‑149 (MQ=255)                  aaTGATATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                                                                                                             <  3:426740/149‑1 (MQ=255)                   aTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                                                                                                            >  1:689411/1‑149 (MQ=255)                         ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                                  >  4:1854/1‑121 (MQ=255)                         ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                                  <  3:1854/121‑1 (MQ=255)                         ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGa                                                                                                                      >  1:532201/1‑149 (MQ=255)                              gACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                                                                                                 >  3:536850/1‑149 (MQ=255)                                  gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                                                                                                          <  2:772150/120‑1 (MQ=255)                                  gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgggg                                                                                                                                          >  1:772150/1‑120 (MQ=255)                                  gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                                                                                                             <  3:26033/149‑1 (MQ=255)                                     tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGat                                                                                                          <  3:544136/149‑1 (MQ=255)                                     tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGat                                                                                                          >  4:318158/1‑149 (MQ=255)                                       gcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                        >  3:445131/1‑149 (MQ=255)                                       gcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGGTTCGCCACTTCTTCCTGGAAGTCCCGCGCCCGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                        >  1:668677/1‑149 (MQ=255)                                        cgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATc                                                                                                       <  4:644765/149‑1 (MQ=255)                                         gAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCt                                                                                                      <  4:170383/149‑1 (MQ=255)                                          aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                     <  2:532201/149‑1 (MQ=255)                                          aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCATGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                     <  4:329030/149‑1 (MQ=255)                                            tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgc                                                                                                   <  2:689411/149‑2 (MQ=255)                                            tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                                                                                                   <  3:318158/149‑1 (MQ=255)                                              tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                                 <  4:385370/149‑1 (MQ=255)                                              tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                                 <  4:756470/149‑1 (MQ=255)                                                    cACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGAt                                                                                            >  1:354910/1‑148 (MQ=255)                                                    cACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGAt                                                                                            <  2:354910/148‑1 (MQ=255)                                                     aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATaa                                                                                          <  2:662501/149‑1 (MQ=255)                                                     aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATaa                                                                                          >  3:180634/1‑149 (MQ=255)                                                     aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATaa                                                                                          >  1:249690/1‑149 (MQ=255)                                                     aCTCTGAACAAGCCACCTGGCAGGCTTTACCGCCGATACAGGTGGTAACGTCGATTAGTTTTGCCACCTCTTCCTGGGAGTCCAGCGGCTGAGGCGCGGGGGTCAGACCGTTAGTCGCAGAACGACGAATGATATCTTGAGATTGATaa                                                                                          <  2:161041/149‑1 (MQ=255)                                                              acGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                     >  2:280231/1‑129 (MQ=255)                                                              acGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                     <  1:280231/129‑1 (MQ=255)                                                              acGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgt                                                                                 >  3:696689/1‑149 (MQ=255)                                                                  cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgtc                                                                             <  2:668677/149‑1 (MQ=255)                                                                  cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgtc                                                                             <  4:445131/149‑1 (MQ=255)                                                                        aCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGtt                                                                       >  4:507552/1‑149 (MQ=255)                                                                         cAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTa                                                                      >  1:124869/1‑149 (MQ=255)                                                                                  cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                             <  4:456969/149‑1 (MQ=255)                                                                                            aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac                                                          <  3:307045/142‑1 (MQ=255)                                                                                            aGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac                                                          >  4:307045/1‑142 (MQ=255)                                                                                             ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaa                                                  <  3:311039/149‑1 (MQ=255)                                                                                                 gtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                         <  3:767888/74‑1 (MQ=255)                                                                                                 gtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                         >  4:767888/1‑74 (MQ=255)                                                                                                    aCGTCGATGAGTTTCGTCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAaggaagga                                           <  3:507552/149‑1 (MQ=255)                                                                                                      gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTaca                                                                    <  3:85250/122‑1 (MQ=255)                                                                                                      gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTaca                                                                    >  4:85250/1‑122 (MQ=255)                                                                                                          aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                             <  1:222695/125‑1 (MQ=255)                                                                                                          aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                             >  2:222695/1‑125 (MQ=255)                                                                                                          aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaaa                                     <  2:249690/149‑1 (MQ=255)                                                                                                               tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCg                                <  2:124869/149‑1 (MQ=255)                                                                                                                   gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGt                            >  4:326008/1‑149 (MQ=255)                                                                                                                     cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCt                          <  4:180634/149‑1 (MQ=255)                                                                                                                         tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                      <  2:747388/149‑1 (MQ=255)                                                                                                                                ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCa               >  3:164204/1‑149 (MQ=255)                                                                                                                                   acccccGCGCCTCAGGCGCGGGAGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTaca                                                                    >  2:727569/4‑93 (MQ=255)                                                                                                                                   acccccGCGCCTCAGGCGCGGGAGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTaca                                                                    <  1:727569/90‑1 (MQ=255)                                                                                                                                     ccccGCGCCTCGGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTACGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTGCGGCGTCTGCGTGTTCGCATCAcc          >  3:527258/2‑149 (MQ=255)                                                                                                                                     ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                      >  3:596806/2‑137 (MQ=255)                                                                                                                                     ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                      <  4:596806/136‑1 (MQ=255)                                                                                                                                     ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc          >  2:106287/2‑149 (MQ=255)                                                                                                                                       ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa        >  2:527094/1‑149 (MQ=255)                                                                                                                                       ccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa        >  3:261397/1‑149 (MQ=255)                                                                                                                                       ccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa        >  1:122646/1‑149 (MQ=255)                                                                                                                                             cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg  >  2:629053/1‑149 (MQ=255)                                                                                                                                             cTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg  >  3:65824/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                   GGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATG  >  NC_000913/4082554‑4082842

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