Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,027,325	C→T	intergenic (+97/+284)	fre → / ← fadA	flavin reductase/3‑ketoacyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,027,325	0	C	T	100.0%	215.6 / NA	68	intergenic (+97/+284)	fre/fadA	flavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/37);  total (31/37)

GAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGT  >  NC_000913/4027190‑4027449                                                                                                                                        |                                                                                                                             gAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATg                                                                                                                          <  2:1106137/140‑1 (MQ=255)      tCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACg                                                                                                                     >  2:111018/1‑140 (MQ=255)         ccTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTg                                                                                                                  >  2:995130/1‑140 (MQ=255)           tttttGGGGATGCTTTTGCATTTATCTGAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTc                                                                                                                <  2:547375/138‑1 (MQ=255)           tGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTc                                                                                                                >  2:90485/1‑140 (MQ=255)           tGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTc                                                                                                                >  2:508805/1‑140 (MQ=255)             tttGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTcgc                                                                                                              >  2:1059049/1‑140 (MQ=255)                  cGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCCGTCGCGTCtt                                                                                                         >  1:308590/1‑140 (MQ=255)                    aTGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                                                       >  1:547375/1‑140 (MQ=255)                     tGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                                                      >  2:100628/1‑140 (MQ=255)                      gCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCg                                                                                                     <  1:302330/140‑1 (MQ=255)                      gCGTTTGATTTTATCTGAGATATAAAAAAGCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCg                                                                                                     <  1:803530/140‑1 (MQ=255)                        gTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTg                                                                                                   <  2:468838/140‑1 (MQ=255)                          ttGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGcc                                                                                                 <  1:401518/140‑1 (MQ=255)                            gCATTTATCTGAGATAAAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTa                                                                                               <  2:177225/140‑1 (MQ=255)                               tttATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACaa                                                                                            <  1:111018/140‑1 (MQ=255)                                ttATCTGAGATATAAAAAACCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACaaa                                                                                           <  1:529641/140‑1 (MQ=255)                                 tATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACaaaa                                                                                          >  1:313559/1‑139 (MQ=255)                                       agaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCg                                                                                    >  1:422475/1‑140 (MQ=255)                                       agaTATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCg                                                                                    >  1:538405/1‑140 (MQ=255)                                            taAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:90485/140‑1 (MQ=255)                                            aaaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:1143935/139‑1 (MQ=255)                                            aaaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:100628/139‑1 (MQ=255)                                            aaaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:531731/139‑1 (MQ=255)                                            aaaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:508805/139‑1 (MQ=255)                                            aaaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCg                                                                               <  1:44062/139‑1 (MQ=255)                                             aaaaaaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAAAATCCCAAACgc                                                                              >  2:585102/1‑139 (MQ=38)                                                 aaccccGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGc                                                                          <  1:585102/137‑1 (MQ=255)                                                 aaaCCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGc                                                                          <  1:950639/140‑1 (MQ=255)                                                      cGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGAt                                                                     <  1:1059049/140‑1 (MQ=255)                                                                  gCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcg                                                         >  2:129397/1‑140 (MQ=255)                                                                   cGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgc                                                        >  1:1043879/1‑140 (MQ=255)                                                                   cGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcgc                                                        >  2:408366/1‑140 (MQ=255)                                                                        agaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                   >  2:795396/1‑140 (MQ=255)                                                                        agaaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                   <  1:129397/140‑1 (MQ=255)                                                                          aaCGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                 >  2:785233/1‑140 (MQ=255)                                                                            cGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                               <  1:785233/140‑1 (MQ=255)                                                                            cGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                               >  1:575465/1‑140 (MQ=255)                                                                               cAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                            <  2:422475/140‑1 (MQ=255)                                                                                aaCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                           <  2:628276/140‑1 (MQ=255)                                                                                aaCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAt                                           <  1:836558/140‑1 (MQ=255)                                                                                          ttATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCatga                                 >  2:647925/1‑140 (MQ=255)                                                                                               cAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCt                            >  1:1039200/1‑140 (MQ=255)                                                                                                aGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                            <  1:806937/123‑1 (MQ=255)                                                                                                aGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGAGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa                                            >  2:806937/1‑123 (MQ=255)                                                                                                  tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGAt                         >  2:1072700/1‑140 (MQ=255)                                                                                                  tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGAt                         >  1:1090342/1‑140 (MQ=255)                                                                                                  tGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGAt                         >  2:118614/1‑140 (MQ=255)                                                                                                   ggCATTTAGATCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATg                        <  2:239378/140‑1 (MQ=255)                                                                                                            aTCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTgg               <  1:647925/140‑2 (MQ=255)                                                                                                            aTCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTgg               <  1:995130/140‑2 (MQ=255)                                                                                                            aTCTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTgg               <  1:1072700/140‑2 (MQ=255)                                                                                                              cTATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGcg             >  2:265342/1‑140 (MQ=255)                                                                                                               tATGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGgcgc            <  2:575465/140‑1 (MQ=255)                                                                                                                aTGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGcgcg           <  2:171192/140‑1 (MQ=255)                                                                                                                 tGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGt          >  2:302736/1‑140 (MQ=255)                                                                                                                 tGACGTATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCAGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGt          >  2:1036736/1‑140 (MQ=255)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  1:795396/140‑1 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  1:782862/140‑1 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     >  2:830428/1‑140 (MQ=25)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  1:669074/140‑1 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  2:1043879/140‑1 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     >  1:556881/1‑140 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     >  1:1080826/1‑140 (MQ=19)                                                                                                                      tATCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  2:390138/140‑1 (MQ=19)                                                                                                                       aTCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     >  1:243600/1‑139 (MQ=21)                                                                                                                       aTCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTAt     <  2:243600/139‑1 (MQ=21)                                                                                                                        tCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATca   >  2:155239/1‑139 (MQ=12)                                                                                                                        tCTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATca   >  1:127816/1‑139 (MQ=12)                                                                                                                         cTGGCAAACCATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATcag  <  2:312931/140‑3 (MQ=2)                                                                                                                                        |                                                                                                                             GAAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAACCCGCCCCTGACAGGCGGGAAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGT  >  NC_000913/4027190‑4027449

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