Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,152,771	T→G	K110Q (AAA→CAA)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,152,771	0	T	G	100.0%	242.4 / NA	64	K110Q (AAA→CAA)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/30);  total (34/30)

TAAACAGGCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCTCGTCGTTCGACAAATTT  >  NC_000913/4152625‑4152916                                                                                                                                                   |                                                                                                                                                  tAAACAGGCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACg                                                                                                                                                <  1:426103/150‑1 (MQ=255)  aaaCAGGCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGc                                                                                                                                               >  1:396501/1‑150 (MQ=255)     cAGGCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAggg                                                                                                                                           >  1:110698/1‑151 (MQ=255)        gCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGc                                                                                                                                        <  2:446845/151‑1 (MQ=255)         cGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGc                                                                                                                                        >  1:358839/1‑150 (MQ=255)         cGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCg                                                                                                                                       <  1:87249/151‑1 (MQ=255)             cACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCg                                                                                                                                   <  1:744901/151‑1 (MQ=255)             cACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCg                                                                                                                                   >  1:62265/1‑151 (MQ=255)                ttCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATc                                                                                                                                >  2:464250/1‑151 (MQ=255)                  ccaccaTTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCa                                                                                                                               <  2:758112/150‑1 (MQ=255)                    accaTTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACtt                                                                                                                            <  1:206088/151‑1 (MQ=255)                     ccaTTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTc                                                                                                                           >  2:595798/1‑151 (MQ=255)                         tttGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCggg                                                                                                                       >  2:39551/1‑151 (MQ=255)                          ttGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGt                                                                                                                      >  1:348142/1‑151 (MQ=255)                           tgtgGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                                                                                     >  2:915580/1‑151 (MQ=255)                           tgtgGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGt                                                                                                                      <  2:637178/150‑1 (MQ=255)                             tgGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGCCAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGtt                                                                                                                     <  2:861473/149‑1 (MQ=255)                               gATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTg                                                                                                                 >  2:325721/1‑151 (MQ=255)                                 tCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGc                                                                                                               >  2:311419/1‑151 (MQ=255)                                      gtACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCtt                                                                                                          <  1:246833/151‑1 (MQ=255)                                          gacgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGcc                                                                                                      >  2:248784/1‑151 (MQ=255)                                            cgGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCtt                                                                                                    <  2:463537/151‑1 (MQ=255)                                                  aTTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCg                                                                                              >  2:246570/1‑151 (MQ=255)                                                      cGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGaaa                                                                                           <  1:781914/150‑1 (MQ=255)                                                         ttGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCt                                                                                       >  2:638749/1‑151 (MQ=255)                                                           gggTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtg                                                                                      <  1:681308/150‑1 (MQ=255)                                                            ggTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCtgtg                                                                                    <  1:543562/151‑1 (MQ=255)                                                               gAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGt                                                                                  <  2:823901/150‑1 (MQ=255)                                                               gAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGt                                                                                  <  2:423852/150‑1 (MQ=255)                                                                  ccGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTg                                                                              >  1:190283/1‑151 (MQ=255)                                                                     tGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTc                                                                           >  2:356049/1‑151 (MQ=255)                                                                        ggACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcggc                                                                        >  1:107923/1‑151 (MQ=255)                                                                          aCCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTcggcgg                                                                      <  1:482794/151‑1 (MQ=255)                                                                             aGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGgtgt                                                                   <  2:918455/151‑1 (MQ=255)                                                                             aGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGgtgt                                                                   <  1:258930/151‑1 (MQ=255)                                                                             aGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGgtgt                                                                   <  1:708813/151‑1 (MQ=255)                                                                                tCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTgg                                                                >  2:160359/1‑151 (MQ=255)                                                                                 ccAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTgg                                                                <  2:259529/150‑1 (MQ=255)                                                                                 ccAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGc                                                               <  2:453889/151‑1 (MQ=255)                                                                                    gcgacagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAg                                                            <  2:601595/151‑1 (MQ=255)                                                                                       acagcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGtt                                                         <  2:382912/151‑1 (MQ=255)                                                                                         agcgaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTc                                                        >  1:431514/1‑150 (MQ=255)                                                                                            gaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAgg                                                     >  1:260518/1‑150 (MQ=255)                                                                                            gaTTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGCTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAgg                                                     >  2:666094/1‑150 (MQ=255)                                                                                              ttCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGaa                                                   >  2:103837/1‑150 (MQ=255)                                                                                              ttCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAAc                                                  >  2:779814/1‑151 (MQ=255)                                                                                                ccACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACt                                                 <  2:408702/150‑1 (MQ=255)                                                                                                   cTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACt                                             >  2:136348/1‑151 (MQ=255)                                                                                                     gCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTa                                            >  2:593391/1‑150 (MQ=255)                                                                                                      cTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTaa                                           >  1:858573/1‑150 (MQ=255)                                                                                                             gATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAacgcacgc                                   <  1:539029/151‑1 (MQ=255)                                                                                                             gATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAacgcacgc                                   >  1:520369/1‑151 (MQ=255)                                                                                                              aTGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCAcgcg                                  >  1:73354/1‑151 (MQ=255)                                                                                                                ggTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCa                                <  2:298595/151‑1 (MQ=255)                                                                                                                     cTTCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACggg                           >  2:255722/1‑151 (MQ=255)                                                                                                                      ttCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGGTCAGGAACTGACTAAACGCACGCGCAACggg                           >  2:630275/1‑150 (MQ=255)                                                                                                                       tCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACgg                            <  2:390637/148‑1 (MQ=255)                                                                                                                       tCGCTCAGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGgc                          >  2:574392/1‑150 (MQ=255)                                                                                                                             aGTTCCGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCt                   <  1:595813/151‑1 (MQ=255)                                                                                                                                 ccGGCTGGTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCtcgtc               <  2:345808/151‑1 (MQ=255)                                                                                                                                       ggTTTTTCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCTCGTCGTTc           <  1:311427/149‑1 (MQ=255)                                                                                                                                          ttttCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACT‑AACGCACGCGCAACGGGCAGCAGCTCGTCGTTCGACaa      >  1:76385/1‑150 (MQ=255)                                                                                                                                           tttCAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCTCGTCGTTCGACaa      >  2:579176/1‑150 (MQ=255)                                                                                                                                              cAGTTGACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCTCGTCGTTCGACAAAttt  >  2:432970/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TAAACAGGCGTTCACTTCCACCATTTTGTGGATCAGTGTACGACGGGTAATTTCGGTTGGGTGAGCCGTGAGGACCAGTTCCAGCGACAGCGATTCCACTGCTTTTTTGATGGTGTCTTCGCTCAGTTCCGGCTGGTTTTTCAGTTTACGCAGGGTGCGGGCGATCACTTCCGGGTTGCTGGCAGCTTCGCCTTTCGGCGAAATGCTGTGGTATTGCTCGGCGGTGTTGGCCAGGTTCAGGAACTGACTAAACGCACGCGCAACGGGCAGCAGCTCGTCGTTCGACAAATTT  >  NC_000913/4152625‑4152916

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