Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,819	C→T	A45T (GCC→ACC)	pepA ←	multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,819	0	C	T	100.0%	314.0 / NA	82	A45T (GCC→ACC)	pepA	multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/35);  total (47/35)

GCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGA  >  NC_000913/4485681‑4485966                                                                                                                                           |                                                                                                                                                    gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAccc                                                                                                                                         >  2:796292/1‑151 (MQ=255) gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAcc                                                                                                                                          <  1:857356/150‑1 (MQ=255)         tGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTg                                                                                                                                  >  1:988868/1‑150 (MQ=255)         tGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTg                                                                                                                                  >  2:15368/1‑150 (MQ=255)          gACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTTATGTACCCATCGCTg                                                                                                                                  <  1:380878/149‑1 (MQ=255)          gACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAt                                                                                                                                >  1:526355/1‑151 (MQ=255)             gCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtttt                                                                                                                             >  1:40212/1‑151 (MQ=255)              cTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtttt                                                                                                                             >  2:357178/1‑150 (MQ=255)              cTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtttt                                                                                                                             >  1:373468/1‑150 (MQ=255)                cATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTAt                                                                                                                           <  1:655035/150‑1 (MQ=255)                cATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATc                                                                                                                          <  2:953363/151‑1 (MQ=255)                     aGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCt                                                                                                                     >  2:60702/1‑151 (MQ=255)                      gCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTg                                                                                                                    >  2:993457/1‑151 (MQ=255)                       cTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGt                                                                                                                   >  1:269989/1‑151 (MQ=255)                         cACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGtt                                                                                                                  <  2:973285/150‑1 (MQ=255)                           cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCt                                                                                                                <  2:892876/150‑1 (MQ=255)                                tttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                                          >  2:723071/1‑151 (MQ=255)                                 ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAAt                                                                                                           <  1:609217/149‑1 (MQ=255)                                 ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAAt                                                                                                           <  2:702977/149‑1 (MQ=255)                                  tGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        >  1:559091/1‑151 (MQ=255)                                   gCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                        <  2:276977/150‑1 (MQ=255)                                      gCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                    <  1:880078/151‑1 (MQ=255)                                       cAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                    >  1:509565/1‑150 (MQ=255)                                        aaCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                    >  2:580613/1‑149 (MQ=255)                                        aaCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                    >  2:85405/1‑149 (MQ=255)                                        aaCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAg                                                                                                  <  2:166385/151‑1 (MQ=255)                                          ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                                <  1:344875/151‑1 (MQ=255)                                          ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                                >  2:241866/1‑151 (MQ=255)                                           cAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                               <  1:273761/151‑1 (MQ=255)                                             aTAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGAc                                                                                             >  2:819714/1‑151 (MQ=255)                                               aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGt                                                                                           >  2:539599/1‑151 (MQ=255)                                               aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTA‑CCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTg                                                                                          >  1:305441/1‑151 (MQ=255)                                                aGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTg                                                                                          <  2:951899/151‑1 (MQ=255)                                                 ggAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTgg                                                                                         >  2:996651/1‑151 (MQ=255)                                                        tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGaa                                                                                   <  2:148325/150‑1 (MQ=255)                                                        tCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGaa                                                                                   <  1:752668/150‑1 (MQ=255)                                                               gAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGccc                                                                            >  2:579800/1‑150 (MQ=255)                                                               gAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCa                                                                           <  2:885109/151‑1 (MQ=255)                                                                aaaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCac                                                                          <  1:150949/151‑1 (MQ=255)                                                                aaaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCac                                                                          <  2:299317/151‑1 (MQ=255)                                                                  aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacga                                                                        <  1:995606/151‑1 (MQ=255)                                                                  aGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacg                                                                         >  1:880778/1‑150 (MQ=255)                                                                    tACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgacg                                                                      >  1:139297/1‑151 (MQ=255)                                                                          cGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGc                                                                  >  2:901719/1‑149 (MQ=255)                                                                           ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAgg                                                               >  1:998430/1‑151 (MQ=255)                                                                             aaCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCa                                                             <  1:574425/151‑1 (MQ=255)                                                                              aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                            >  1:868383/1‑151 (MQ=255)                                                                               cATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACt                                                           >  1:351669/1‑151 (MQ=255)                                                                               cATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACt                                                           <  2:567869/151‑1 (MQ=255)                                                                                aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTc                                                          >  2:88080/1‑151 (MQ=255)                                                                                aTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTc                                                          <  2:696090/151‑1 (MQ=255)                                                                                  ggTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCg                                                        >  1:664645/1‑151 (MQ=255)                                                                                      cagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGt                                                    <  1:426053/151‑1 (MQ=255)                                                                                      cagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGt                                                    >  1:788823/1‑151 (MQ=255)                                                                                      cagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGt                                                    <  2:977810/151‑1 (MQ=255)                                                                                       agcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGtt                                                   >  1:796578/1‑151 (MQ=255)                                                                                       agcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGtt                                                   <  1:968299/151‑1 (MQ=255)                                                                                          aacaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtct                                                >  2:264007/1‑151 (MQ=255)                                                                                             aaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCg                                             >  2:669601/1‑151 (MQ=255)                                                                                              aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCgg                                            >  1:234574/1‑151 (MQ=255)                                                                                              aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCgg                                            >  2:323832/1‑151 (MQ=255)                                                                                                 tCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                         >  2:62387/1‑151 (MQ=255)                                                                                                    gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTa                                        <  1:637055/149‑1 (MQ=255)                                                                                                        cGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtt                                  <  1:714900/151‑1 (MQ=255)                                                                                                        cGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCtt                                  >  1:520071/1‑151 (MQ=255)                                                                                                         ggTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttt                                 >  2:738749/1‑151 (MQ=255)                                                                                                         ggTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttt                                 >  2:929949/1‑151 (MQ=255)                                                                                                              tccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTCCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac                             >  2:53941/1‑150 (MQ=255)                                                                                                               ccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                           >  2:870722/1‑151 (MQ=255)                                                                                                                cttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                          >  2:570366/1‑151 (MQ=255)                                                                                                                   ccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa                       <  2:140299/151‑1 (MQ=255)                                                                                                                   ccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa                       <  1:99419/151‑1 (MQ=255)                                                                                                                     aGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACt                     <  2:309595/151‑1 (MQ=255)                                                                                                                        tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                   >  1:452466/1‑150 (MQ=255)                                                                                                                        tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                   >  1:552289/1‑150 (MQ=255)                                                                                                                           ccccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg                >  1:340860/1‑150 (MQ=255)                                                                                                                            cccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACCACGATGCAGGCACTCCGCTGCTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc               >  1:64359/1‑150 (MQ=255)                                                                                                                             ccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc               >  2:261784/1‑149 (MQ=255)                                                                                                                             ccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCAc             <  2:941345/151‑1 (MQ=255)                                                                                                                                      cAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCt    <  2:471755/151‑1 (MQ=255)                                                                                                                                        gggTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGa  <  1:540316/151‑1 (MQ=255)                                                                                                                                        gggTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGa  <  2:137513/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                    GCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGA  >  NC_000913/4485681‑4485966

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