Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,355,453	T→C	L6L (TTA→TTG)	sapB ←	putrescine ABC exporter membrane subunit SapB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,355,453	0	T	C	100.0%	136.6 / NA	41	L6L (TTA→TTG)	sapB	putrescine ABC exporter membrane subunit SapB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/21);  total (20/21)

AAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGTAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAAGAGAGCGCCTTACGCAATACGCTG  >  NC_000913/1355210‑1355695                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   aaaaCCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTg                                                                                                                                                                                                                                              >  1:367533/1‑250 (MQ=255)           gCGAGGATGCACAGTTCCATCGTCGCAGGGAAGACCTCTGTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCTCAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGACAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAAtaat                                                                                                                                                                                                                                   <  2:110317/251‑4 (MQ=255)             gAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATgg                                                                                                                                                                                                                                 >  2:351818/1‑251 (MQ=255)                 aTGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGtttt                                                                                                                                                                                                                             <  2:285035/251‑1 (MQ=255)                   gCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGtttttt                                                                                                                                                                                                                           <  2:143538/251‑1 (MQ=255)                   gCACAGTTCCATCGTCGCCGGGAAGAACTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCCAGGTGAGGATAATCATGGtttttt                                                                                                                                                                                                                           >  2:176250/1‑251 (MQ=255)                   gCACAGTTCCATCGTCGCAGGGAAGACCTCTTTTAACTGTTCGGCAATCGTTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGAGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGtttttt                                                                                                                                                                                                                           >  2:176253/1‑251 (MQ=255)                             aTCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATcc                                                                                                                                                                                                                 >  2:63102/1‑251 (MQ=255)                                    ccGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTtct                                                                                                                                                                                                          >  2:29394/1‑251 (MQ=255)                                      gggAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTcg                                                                                                                                                                                                        <  2:110841/251‑1 (MQ=255)                                                     aaCTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCa                                                                                                                                                                                         <  1:249079/251‑1 (MQ=255)                                                                aaTCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGtt                                                                                                                                                                              >  1:287505/1‑251 (MQ=255)                                                                         aCCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACgg                                                                                                                                                                     >  2:250527/1‑251 (MQ=255)                                                                           cGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGc                                                                                                                                                                   >  1:17293/1‑251 (MQ=255)                                                                                  gCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTAc                                                                                                                                                            <  1:227832/251‑1 (MQ=255)                                                                                                            tCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAgg                                                                                                                                  >  2:172672/1‑251 (MQ=255)                                                                                                             cAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAgg                                                                                                                                  <  1:350459/250‑1 (MQ=255)                                                                                                              aGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACTGATCGGcc                                                                                                                                >  2:11899/1‑251 (MQ=255)                                                                                                                     tAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCaaa                                                                                                                          <  2:320719/250‑1 (MQ=255)                                                                                                                     tAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGAAATCGTACCGATAGGCCTGCAAAc                                                                                                                         >  1:112770/1‑251 (MQ=255)                                                                                                                      aaaCCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACg                                                                                                                        >  2:61391/1‑251 (MQ=255)                                                                                                                                        aTTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCgg                                                                                                      <  1:221592/251‑1 (MQ=255)                                                                                                                                                gTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATgg                                                                                              <  1:116785/251‑1 (MQ=255)                                                                                                                                                            gTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTgc                                                                                   <  2:303155/250‑1 (MQ=255)                                                                                                                                                                     cATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATAc                                                                         >  2:140952/1‑251 (MQ=255)                                                                                                                                                                                aaaTAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCg                                                              <  2:314008/251‑1 (MQ=255)                                                                                                                                                                                    aGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCAt                                                          <  2:78580/251‑1 (MQ=255)                                                                                                                                                                                     gCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGc                                                                                            >  2:217070/1‑216 (MQ=255)                                                                                                                                                                                     gCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGc                                                                                            <  1:217070/216‑1 (MQ=255)                                                                                                                                                                                              gAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCa                                                  <  2:280142/249‑1 (MQ=255)                                                                                                                                                                                              gAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAAt                                                <  2:118265/251‑1 (MQ=255)                                                                                                                                                                                               aaGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATa                                               <  1:310996/251‑1 (MQ=255)                                                                                                                                                                                                              agcagGAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAgctgc                                >  1:266159/1‑251 (MQ=255)                                                                                                                                                                                                                   gAACAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCg                            >  1:359921/1‑250 (MQ=255)                                                                                                                                                                                                                    aaCAAGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGaa                          <  1:107590/251‑1 (MQ=255)                                                                                                                                                                                                                       aaGGTGACAATCAATAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGAAATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAagag                       <  2:112770/251‑1 (MQ=255)                                                                                                                                                                                                                                atcaatAACAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAAGAGAGCGCCTTa              <  1:140952/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACTTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGt                                                                                                                                           <  2:181316/118‑1 (MQ=255)                                                                                                                                                                                                                                        cAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGt                                                                                                                                           >  1:181316/1‑118 (MQ=255)                                                                                                                                                                                                                                        cAAAATGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAAGAGAGCGCCTTACGCAATAc      >  1:206481/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTGCGGCGCAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAAGAGAGCGCCTTACGCAATACGCTg  >  1:1253/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   AAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAACACCCAGGCATTCCACAGTGACGCGCCTTGTAACGGCGCATGAGGGGTGAAATAGCTCAAGCTGAAGCCAACGAAGGTCAGCAGGAACAAGGTGACAATCAATAACAAAATGCGGCGTAAGGTGAAGATAATCATGGTTTTTTCACCTCATCCTGTTTCTCGCGATACACCCCAGCAAAGGAGGCGTTACCAAACGGGCTAAGTACCAGACCTTTGATATCGTACCGATAGGCCTGCAAACGCAATGACGACGCCAGCGGCAAAATGGGCAATTCCTGCGCCAGAATACTCTGCGCTTCGTCATAGGCTTCAATACGCGCCGCCAGCTGCTGCGAAGAGAGCGCCTTACGCAATACGCTG  >  NC_000913/1355210‑1355695

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