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%	308.2 / 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)

GAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCT  >  NC_000913/4485672‑4485964                                                                                                                                                    |                                                                                                                                                  gAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGgtgc                                                                                                                                                 <  2:871345/150‑1 (MQ=255)     aaCCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATg                                                                                                                                            >  2:185205/1‑151 (MQ=255)        cTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAc                                                                                                                                         >  1:138911/1‑151 (MQ=255)          gCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAcc                                                                                                                                        <  1:1114279/150‑1 (MQ=255)           cTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTAccc                                                                                                                                       >  2:846184/1‑150 (MQ=255)            ttGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCAt                                                                                                                                     >  2:9373/1‑151 (MQ=255)             tGTACTGACGCTCATCCAGCTCTCGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGTAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCAAt                                                                                                                                     <  2:1005857/150‑1 (MQ=255)                   gACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTg                                                                                                                                <  2:1168930/149‑1 (MQ=255)                    tcGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtt                                                                                                                             >  2:574259/2‑151 (MQ=255)                         caaccAGCGCTCGTTCTTTGCCACAACCAATAAGGAGAATTCGCGCGGAAAGTACATTCGGACTATGGTGCCGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGAtttt                                                                                                                           <  2:1233211/145‑1 (MQ=255)                          aTCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCg                                                                                                                       <  1:402976/151‑1 (MQ=255)                            ccAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAg                                                                                                                     <  2:994840/151‑1 (MQ=255)                              aGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCt                                                                                                                   >  1:350154/1‑151 (MQ=255)                               gCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTg                                                                                                                  >  1:573943/1‑151 (MQ=255)                                  cACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGtt                                                                                                                >  2:79449/1‑150 (MQ=255)                                   aCGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCt                                                                                                              >  1:250763/1‑151 (MQ=255)                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                             <  1:902626/151‑1 (MQ=255)                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                             >  2:903301/1‑151 (MQ=255)                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTg                                                                                                             <  1:1240852/151‑1 (MQ=255)                                    cGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTATTCt                                                                                                              <  2:802221/150‑1 (MQ=255)                                         tttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATc                                                                                                        >  1:1036975/1‑151 (MQ=255)                                          ttGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCaa                                                                                                          <  1:526880/148‑1 (MQ=255)                                           tGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                      >  2:504995/1‑151 (MQ=255)                                           tGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                      >  2:216758/1‑151 (MQ=255)                                             ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                      >  1:557182/1‑149 (MQ=255)                                             ccGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCgg                                                                                                      >  2:290753/1‑149 (MQ=255)                                                cAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                  >  2:1247529/1‑150 (MQ=255)                                                cAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                  <  2:29415/150‑1 (MQ=255)                                                 aaCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAACGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGaa                                                                                                  >  1:43680/1‑149 (MQ=255)                                                  aCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                                               <  2:949667/151‑1 (MQ=255)                                                  aCCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAgg                                                                                               <  1:191681/151‑1 (MQ=255)                                                   ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                              <  1:373104/151‑1 (MQ=255)                                                   ccAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGc                                                                                              >  1:1225141/1‑151 (MQ=255)                                                    cAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCg                                                                                             >  1:336917/1‑151 (MQ=255)                                                        aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACg                                                                                          <  2:6735/150‑1 (MQ=255)                                                        aaGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGt                                                                                         >  2:739289/1‑151 (MQ=255)                                                         aGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTg                                                                                        >  2:58704/1‑151 (MQ=255)                                                            agaATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGtt                                                                                     >  1:696428/1‑151 (MQ=255)                                                                  cGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAg                                                                                <  2:656970/150‑1 (MQ=255)                                                                  cGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAg                                                                                <  2:1101647/150‑1 (MQ=255)                                                                     tCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGc                                                                            >  1:1031503/1‑151 (MQ=255)                                                                     tCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGc                                                                            >  2:1132468/1‑151 (MQ=255)                                                                      cGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACg                                                                             <  1:36161/149‑1 (MQ=255)                                                                      cGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGcc                                                                           >  2:242261/1‑151 (MQ=255)                                                                          aaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacg                                                                       >  2:433961/1‑151 (MQ=255)                                                                          aaGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacg                                                                       >  2:929919/1‑151 (MQ=255)                                                                              aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCacgac                                                                     >  1:480907/1‑149 (MQ=255)                                                                                aTTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATg                                                                 >  1:499798/1‑151 (MQ=255)                                                                                aTTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATg                                                                 >  2:70284/1‑151 (MQ=255)                                                                                  tCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGc                                                                >  2:568982/1‑150 (MQ=255)                                                                                   cGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCa                                                               <  2:779258/150‑1 (MQ=255)                                                                                     gAACCTGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGGCCCGACGTAGCAGGGTGCTGATGTTCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCa                                                               <  1:402147/148‑1 (MQ=255)                                                                                       aCATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCAc                                                          >  1:610496/1‑151 (MQ=255)                                                                                               cagcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGt                                                  <  1:96320/151‑1 (MQ=255)                                                                                                 gcaACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGttt                                                >  2:71612/1‑151 (MQ=255)                                                                                                     caaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTcc                                            >  1:458526/1‑151 (MQ=255)                                                                                                       aTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCgg                                          <  2:961021/151‑1 (MQ=255)                                                                                                         gTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCggg                                         <  1:574282/150‑1 (MQ=255)                                                                                                          tCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                       >  1:439073/1‑151 (MQ=255)                                                                                                           cTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                       >  2:184207/1‑150 (MQ=255)                                                                                                           cTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTCTTTCTCCGGGCt                                       >  2:226918/1‑150 (MQ=255)                                                                                                            tGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAc                                     >  1:557701/1‑151 (MQ=255)                                                                                                             gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTa                                      <  1:30729/149‑1 (MQ=255)                                                                                                                ccGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACctct                                 <  2:1013508/151‑1 (MQ=255)                                                                                                                  ggTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttt                               >  2:40090/1‑151 (MQ=255)                                                                                                                     tttCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt                             >  2:460271/1‑150 (MQ=255)                                                                                                                     tttCCTTCCAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt                             >  1:504744/1‑150 (MQ=255)                                                                                                                      ttccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTac                           >  1:1013343/1‑151 (MQ=255)                                                                                                                        ccttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                         >  2:795746/1‑151 (MQ=255)                                                                                                                         cttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                        >  2:29288/1‑151 (MQ=255)                                                                                                                          ttccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                        >  2:514979/1‑150 (MQ=255)                                                                                                                           tccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACt                        <  2:458515/149‑1 (MQ=255)                                                                                                                           tccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTa                       <  1:558270/150‑1 (MQ=255)                                                                                                                            ccAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa                      <  1:605472/150‑1 (MQ=255)                                                                                                                             cAGTTCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAAc                    <  1:846187/151‑1 (MQ=255)                                                                                                                                 tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACtac                 <  2:1219025/150‑3 (MQ=255)                                                                                                                                 tCGCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTcc                 >  1:663441/1‑150 (MQ=255)                                                                                                                                   gCCCCGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATg              <  1:700337/151‑1 (MQ=255)                                                                                                                                     cccGACGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCa            <  2:1036980/151‑1 (MQ=255)                                                                                                                                          cGTAGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACg       <  1:573583/151‑1 (MQ=255)                                                                                                                                             aGCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCt     <  2:260566/150‑1 (MQ=255)                                                                                                                                              gCAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTcc   <  1:566531/151‑1 (MQ=255)                                                                                                                                               cAGGGTGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCt  >  2:136064/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GAATAACCTGCTTGTACTGACGCTCATCCAGCTCACGTTCTTTGCCGCAACCAATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCT  >  NC_000913/4485672‑4485964

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