Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,206,742	T→C	N264S (AAC→AGC)	purH ←	bifunctional AICAR transformylase/IMP cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,206,742	0	T	C	100.0%	246.6 / NA	71	N264S (AAC→AGC)	purH	bifunctional AICAR transformylase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/41);  total (30/41)

GCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCT  >  NC_000913/4206613‑4206868                                                                                                                                  |                                                                                                                               gcgcgATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAaga                                                                                                                      <  1:351048/140‑1 (MQ=255)   gcgATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAagaga                                                                                                                    <  2:1471116/140‑1 (MQ=255)   gcgATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAagaga                                                                                                                    >  1:33346/1‑140 (MQ=255)        cGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCt                                                                                                               <  2:481457/140‑1 (MQ=255)         gTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCtt                                                                                                              <  2:691334/140‑1 (MQ=255)         gTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCtt                                                                                                              <  1:913386/140‑1 (MQ=255)              cATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTAccc                                                                                                         <  2:1225783/140‑1 (MQ=255)              cATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTAccc                                                                                                         <  2:541290/140‑1 (MQ=255)              cATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTAccc                                                                                                         <  2:240190/140‑1 (MQ=255)                 cAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGa                                                                                                      >  1:586739/1‑140 (MQ=255)                 cAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGa                                                                                                      >  1:283100/1‑140 (MQ=255)                 cAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGa                                                                                                      >  1:1169134/1‑140 (MQ=255)                 cAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGa                                                                                                      >  1:116192/1‑140 (MQ=255)                    gAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAAcc                                                                                                   >  1:448273/1‑140 (MQ=255)                    gAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAAcc                                                                                                   <  1:645856/140‑1 (MQ=255)                        ggAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTggg                                                                                               >  1:736218/1‑140 (MQ=255)                         gAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGt                                                                                              >  2:755968/1‑140 (MQ=255)                            ttGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGc                                                                                           >  1:1306023/1‑140 (MQ=255)                             tGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGACCTCGTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCg                                                                                          <  2:1207602/140‑1 (MQ=255)                             tGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCg                                                                                          <  1:180278/140‑1 (MQ=255)                              gCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCgg                                                                                         >  1:849404/1‑140 (MQ=255)                              gCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCgg                                                                                         >  1:1109719/1‑140 (MQ=255)                               ccGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGt                                                                                        <  2:1243669/140‑1 (MQ=255)                                    aGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCaa                                                                                   <  2:33346/140‑1 (MQ=255)                                         cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACgg                                                                                <  1:1207132/138‑1 (MQ=255)                                         cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACgg                                                                                >  2:1207132/1‑138 (MQ=255)                                         cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAg                                                                              >  2:848657/1‑140 (MQ=255)                                         cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAg                                                                              >  1:179194/1‑140 (MQ=255)                                         cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAg                                                                              >  2:1189915/1‑140 (MQ=255)                                         cgccgcAAGGGTTGACGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAg                                                                              <  2:465365/140‑1 (MQ=255)                                             gcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGctt                                                                          <  2:179194/140‑1 (MQ=255)                                                  ggTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTc                                                                     <  2:284861/140‑1 (MQ=255)                                                     tGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                  <  2:363693/140‑1 (MQ=255)                                                     tGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                  <  2:1097880/140‑1 (MQ=255)                                                      ggCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAt                                                                 >  1:96586/1‑140 (MQ=255)                                                      ggCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAt                                                                 >  2:834823/1‑140 (MQ=255)                                                       gCGTGCTTCGCAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                <  2:1306023/140‑1 (MQ=255)                                                       gCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                <  2:586739/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCTGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATtct                                                              <  1:1380510/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAGCGGAGGCTTCTTTCACATtct                                                              <  2:1048347/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATtct                                                              <  1:755968/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATtct                                                              <  1:848657/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATtct                                                              <  1:834823/140‑1 (MQ=255)                                                         gTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATtct                                                              <  1:750336/140‑1 (MQ=255)                                                            cTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtctt                                                           >  2:1364163/1‑140 (MQ=255)                                                            cTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtctt                                                           >  2:1011132/1‑140 (MQ=255)                                                            cTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtctt                                                           <  1:1011132/140‑1 (MQ=255)                                                            cTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtctt                                                           <  1:1364163/140‑1 (MQ=255)                                                             ttCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                >  2:986588/1‑134 (MQ=255)                                                             ttCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                <  1:986588/134‑1 (MQ=255)                                                              tCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtcttct                                                         <  2:116192/140‑1 (MQ=255)                                                                  aaTCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCtatat                                                     <  1:1125155/140‑1 (MQ=255)                                                                       cacaTGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAg                                                >  1:1162415/1‑140 (MQ=255)                                                                         caTGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGc                                              >  1:520465/1‑140 (MQ=255)                                                                            gCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGc                                           <  2:1169134/140‑1 (MQ=255)                                                                            gCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGc                                           <  2:736218/140‑1 (MQ=255)                                                                                  tCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAgg                                               >  1:891875/1‑130 (MQ=255)                                                                                  tCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAgg                                               <  2:891875/130‑1 (MQ=255)                                                                                  tCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCctgctg                                     >  2:49504/1‑140 (MQ=255)                                                                                          ctctTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGt                             >  1:248372/1‑140 (MQ=255)                                                                                             tttCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTtct                          <  2:1109719/140‑1 (MQ=255)                                                                                              ttCACGCACTCCAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTtctc                         >  1:520188/1‑140 (MQ=255)                                                                                                         cAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCa              <  2:520188/140‑1 (MQ=255)                                                                                                         cAGCGCCGCATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCa              >  1:639760/1‑140 (MQ=255)                                                                                                               cgcATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCa        <  1:49504/140‑1 (MQ=255)                                                                                                               cgcATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCa        <  2:283100/140‑1 (MQ=255)                                                                                                                 cATCGGTATCGGCGATGCTGTTATAAGAGAGGGTTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc      >  1:1202613/1‑140 (MQ=255)                                                                                                                 cATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc      >  1:324867/1‑140 (MQ=255)                                                                                                                 cATCGGTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc      >  2:889082/1‑140 (MQ=255)                                                                                                                     ggTATCGGCGATGCTGTTATAAGAGAGGGTTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGcttct  <  2:1202613/140‑1 (MQ=255)                                                                                                                     ggTATCGGCGATGCTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGcttct  <  2:639760/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                               GCGCGATCGTAAGCATCAAGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCT  >  NC_000913/4206613‑4206868

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