Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,113,731	C→T	D208N (GAT→AAT)	pppA ←	prepilin peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,113,731	0	C	T	100.0%	133.2 / NA	40	D208N (GAT→AAT)	pppA	prepilin peptidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/12);  total (28/12)

AGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCA  >  NC_000913/3113507‑3113973                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   agTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGc                                                                                                                                                                                                                           >  1:468456/1‑250 (MQ=255)    tAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAAAGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCGTCTTTACGCAGa                                                                                                                                                                                                                        >  2:538893/1‑250 (MQ=255)       aaTGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACt                                                                                                                                                                                                                     >  1:469526/1‑250 (MQ=255)        aTGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTAt                                                                                                                                                                                                                   <  2:215101/251‑1 (MQ=255)          gTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGCGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTc                                                                                                                                                                                                                 >  2:249047/1‑251 (MQ=255)                  gATAATGGCTTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTTGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGACCTATTCCGTCTATc                                                                                                                                                                                                         <  2:26488/251‑1 (MQ=255)                    tAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATcc                                                                                                                                                                                                        >  1:523464/1‑250 (MQ=255)                              gTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCagag                                                                                                                                                                                               >  1:481544/1‑249 (MQ=255)                                    aTTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGCGTaaaa                                                                                                                                                                                        >  2:566900/1‑250 (MQ=255)                                            aaTGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATaa                                                                                                                                                                                >  2:443324/1‑250 (MQ=255)                                                   gTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCccacc                                                                                                                                                                        >  1:607531/1‑251 (MQ=255)                                                        tAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGAc                                                                                                                                                                   >  2:71811/1‑251 (MQ=255)                                                             tgtTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCgg                                                                                                                                                              >  1:217920/1‑251 (MQ=255)                                                                    ataCCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCGAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACt                                                                                                                                                         <  1:509200/249‑1 (MQ=255)                                                                     taCCGCCCCGACTTAAGCACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAGAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCTGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGc                                                                                                                                                       >  2:34615/1‑250 (MQ=255)                                                                             aGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGt                                                                                                                                                >  1:121166/1‑249 (MQ=255)                                                                                      cacGGTCCAAAAGGCAGTGTGGTTGATACTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCgg                                                                                                                                     <  1:194995/251‑1 (MQ=255)                                                                                            ccAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTctg                                                                                                                               >  1:534029/1‑251 (MQ=255)                                                                                                   gCAGTGTGGTTGATCCTCTGTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGcc                                                                                                                        <  1:538893/251‑1 (MQ=255)                                                                                                                                            aGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGGGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGGCCTTTCCACAATACGCCCTGAGTaaa                                                                               <  2:401420/251‑1 (MQ=255)                                                                                                                                               ccGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAAc                                                                            >  1:197515/1‑251 (MQ=255)                                                                                                                                                  cagcaTGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATc                                                                         >  2:206716/1‑251 (MQ=255)                                                                                                                                                  cagcaTGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGAGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATc                                                                         >  1:355256/1‑251 (MQ=255)                                                                                                                                                     caTGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCGCCCAGCCTCCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCggg                                                                      <  1:429564/251‑1 (MQ=255)                                                                                                                                                       tGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGc                                                                     >  1:211831/1‑250 (MQ=255)                                                                                                                                                            cGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACCATACGCCCTGAGTAAAAACATCGGGCAGCCAt                                                               >  2:189314/1‑251 (MQ=255)                                                                                                                                                              aTTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTATCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTg                                                             >  1:343547/1‑251 (MQ=255)                                                                                                                                                                          ttGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTTCGCCCATGCCGCAATCa                                                                                                           <  2:611831/193‑1 (MQ=255)                                                                                                                                                                          ttGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCa                                                                                                           >  1:611831/1‑193 (MQ=255)                                                                                                                                                                               cAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATg                                             >  1:88739/1‑250 (MQ=255)                                                                                                                                                                                      aaTGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTcgc                                     <  1:71811/251‑1 (MQ=255)                                                                                                                                                                                       aTGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTcgcg                                    <  2:228109/251‑1 (MQ=255)                                                                                                                                                                                            cccacccaGCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTATGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAAt                               >  1:96488/1‑251 (MQ=255)                                                                                                                                                                                                   aGCNACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGTCTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGGAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACGCTCGCGTCAATCAGCCAg                        <  2:595566/251‑1 (MQ=255)                                                                                                                                                                                                    gCCACCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAgg                       >  1:74802/1‑251 (MQ=255)                                                                                                                                                                                                       aCCTAACGCGGCGAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCgg                    >  2:71284/1‑251 (MQ=255)                                                                                                                                                                                                                  cgAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGCACTATTCAGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGCGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATatcatc         >  1:598579/1‑251 (MQ=255)                                                                                                                                                                                                                   gAAAAGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCACGGCAATCAGCCAGGCGGATAATatcatc         >  1:265527/1‑250 (MQ=255)                                                                                                                                                                                                                     aaaGTAATACATTGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCAcc      >  1:103376/1‑251 (MQ=255)                                                                                                                                                                                                                          aaTACATTGCCCATGCCTAATGCTTCTTTACGCAAAACTATTCCGGCTATCCAACGCAGAAAGTAAAAAGTGATAATGCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACGCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCa  <  2:355256/250‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   AGTTAAAATGTCACTTTGATAATGACGTGGTTATCATTAAAACAATGCCTGTAGATAAAGTGTTGCTATACCGCCCAGACTTAAACACGGTCCAAAAGGCAGTGTGGTTGATCCTCTTTTTGTAATAACGGCATATATCAGGCCGCAGCATGAGGCGATTAAAGCAACATTGGGCAGCGACAATGCCCCCACCCAGCCACCTAACGCGGCGAAAAGTAATACATCGCCCATGCCTAATGCTTCTTTACGCAGAACTATTCCGGCTATCCAACGCAGAGAGTAAAAAGTGATAAATCCCACCAGGACGCCGGTGACTGCATCTTGTAGCGTGAGCGGACTCTGCTGCGCCCATGCCGCAATCAGTCCTGTCCACAATACGCCCTGAGTAAAAACATCGGGCAGCCATTGGTGATCGAGGTCGATGACACTCGCGGCAATCAGCCAGGCGGATAATATCATCACCGCCA  >  NC_000913/3113507‑3113973

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