Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	353,717	C→G	6.1%	R338G (CGC→GGC)	prpE →	propionate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	353,717	0	C	G	6.1%	213.3 / 6.1	66	R338G (CGC→GGC)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base G (2/2);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

TTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAACTA  >  NC_000913/353572‑353862                                                                                                                                                  |                                                                                                                                                  ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTcgcg                                                                                                                                                <  3:155682/149‑1 (MQ=255)    gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTg                                                                                                                                              <  6:50013/148‑1 (MQ=255)    gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGc                                                                                                                                             >  7:266817/1‑149 (MQ=255)    gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGc                                                                                                                                             >  3:56129/1‑149 (MQ=255)      atatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCt                                                                                                                                            >  2:78757/1‑148 (MQ=255)       tatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTg                                                                                                                                           >  3:159341/1‑148 (MQ=255)       tatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTg                                                                                                                                           >  3:95692/1‑148 (MQ=255)           gTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGaaaa                                                                                                                                       >  6:6357/1‑148 (MQ=255)            tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGaaaaaa                                                                                                                                     >  7:326048/1‑149 (MQ=255)               aCGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTc                                                                                                                                  >  8:118593/1‑149 (MQ=255)               aCGCGCCGCTGCTGGCGGGGATGGCGACTAACTTATACGAAGTATTGCCGACATGGCCGGACTGCGGCGTGAGGTGGACAATCGTCGAGAAATATCAGGTTAGCCCGATGATCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAAtta                                                                                                                                  >  4:271449/1‑148 (MQ=255)                 gcgcCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                                                >  4:207110/1‑149 (MQ=255)                 gcgcCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTAGCGTGCTGAAAAAATTccc                                                                                                                                >  3:208743/1‑149 (MQ=255)                 gcgcCGCTGCTGGCGGGGATGGCGACTAGCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                                                <  8:92207/149‑1 (MQ=255)                    ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGCCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGTTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTAc                                                                                                                             >  6:290491/1‑149 (MQ=255)                    ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTAc                                                                                                                             <  5:217140/149‑1 (MQ=255)                               ggggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  <  3:102359/149‑1 (MQ=255)                                gggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                                                 <  4:56129/149‑1 (MQ=255)                                gggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                                                 >  4:33467/1‑149 (MQ=255)                                 ggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGc                                                                                                                <  1:31561/149‑1 (MQ=255)                                      gCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAaca                                                                                                           >  5:81834/1‑149 (MQ=255)                                           tATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtc                                                                                                      <  3:126249/149‑1 (MQ=255)                                             tCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    >  4:168778/1‑149 (MQ=255)                                             tCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    >  4:164544/1‑149 (MQ=255)                                               gTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                                                  >  2:7277/1‑149 (MQ=255)                                               gTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctct                                                                                                   <  7:98332/148‑1 (MQ=255)                                                  tACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCtcgt                                                                                                <  3:207110/148‑1 (MQ=255)                                                    cGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAac                                                                                                            >  2:50180/1‑134 (MQ=255)                                                    cGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAac                                                                                                            <  1:50180/134‑1 (MQ=255)                                                     gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  >  5:319936/1‑127 (MQ=255)                                                     gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  <  6:319936/127‑1 (MQ=255)                                                     gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCt                                                                                            >  7:40207/1‑149 (MQ=255)                                                         gATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                                                        <  1:239030/149‑1 (MQ=255)                                                              ccGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGct                                                                                   <  1:286802/149‑1 (MQ=255)                                                               cGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGctc                                                                                  >  4:67043/1‑149 (MQ=255)                                                                   cTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATc                                                                              <  5:11582/149‑1 (MQ=255)                                                                        cGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggct                                                                         >  7:15287/1‑149 (MQ=255)                                                                         ggACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggctg                                                                        <  8:9691/149‑1 (MQ=255)                                                                           aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGga                                                                      <  3:168778/149‑1 (MQ=255)                                                                           aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGga                                                                      <  3:164544/149‑1 (MQ=255)                                                                           aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggctg                                                                        >  3:120914/1‑147 (MQ=255)                                                                            cTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                                                     <  3:271449/149‑1 (MQ=255)                                                                            cTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                                                     <  7:118593/149‑1 (MQ=255)                                                                               cggcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGaa                                                                   <  1:78757/148‑1 (MQ=255)                                                                               cggcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAAc                                                                  >  1:56782/1‑149 (MQ=255)                                                                                 gcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCg                                                                <  8:326048/149‑1 (MQ=255)                                                                                  cgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGc                                                               <  4:159341/149‑1 (MQ=255)                                                                                          ggACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAg                                                       >  1:130052/1‑149 (MQ=255)                                                                                                tcgtcgAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGAcc                                                 <  6:49102/149‑1 (MQ=255)                                                                                                  gtcgAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCt                                                                                            <  6:276071/104‑1 (MQ=255)                                                                                                     gagaAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCa                                             <  2:130052/148‑1 (MQ=255)                                                                                                          atatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTggg                                       <  1:39419/149‑1 (MQ=255)                                                                                                          atatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCGCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg                                        <  3:33467/148‑1 (MQ=255)                                                                                                           tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    <  2:143217/87‑1 (MQ=255)                                                                                                           tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      <  1:7277/149‑1 (MQ=255)                                                                                                           tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      <  8:266817/149‑1 (MQ=255)                                                                                                             tCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      >  3:178241/1‑147 (MQ=255)                                                                                                               aGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAg                                   <  4:120914/148‑1 (MQ=255)                                                                                                                    aGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATa                              <  3:67043/148‑1 (MQ=255)                                                                                                                      ccGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGc                           <  5:6357/149‑1 (MQ=255)                                                                                                                          aTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGCGGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGa                       <  4:95692/149‑1 (MQ=255)                                                                                                                               ctcAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAACTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGATAACCGCTGGACGAGCCGACCGCCACTTGGATGAGCAATACGCTGGATGTGc                  >  4:184332/1‑149 (MQ=255)                                                                                                                                  aGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGc                  <  7:37558/146‑1 (MQ=255)                                                                                                                                       cGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATc          <  7:72484/149‑1 (MQ=255)                                                                                                                                          ccGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGAc       <  8:40207/149‑1 (MQ=255)                                                                                                                                            gCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACaa     <  3:205415/149‑1 (MQ=255)                                                                                                                                               aTTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAacta  >  5:217880/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAACTA  >  NC_000913/353572‑353862

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