Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,333,691	A→G	pseudogene (712/4605 nt)	yfaS ←	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,333,691	0	A	G	100.0%	238.8 / NA	63	pseudogene (712/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/34);  total (29/34)

CTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCA  >  NC_000913/2333547‑2333834                                                                                                                                                 |                                                                                                                                                ctGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGt                                                                                                                                             <  2:581047/149‑1 (MQ=255) ctGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGt                                                                                                                                             <  2:169868/149‑1 (MQ=255) ctGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGt                                                                                                                                             <  1:765239/149‑1 (MQ=255)  tGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTccc                                                                                                                                          >  1:302705/1‑151 (MQ=255)   gAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCg                                                                                                                                         >  1:547689/1‑151 (MQ=255)    aGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCgg                                                                                                                                        <  1:381684/151‑1 (MQ=255)        gggCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGctct                                                                                                                                    <  2:585290/151‑1 (MQ=255)        gggCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGctct                                                                                                                                    <  2:282237/151‑1 (MQ=255)                   cACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCg                                                                                                                           >  1:120698/1‑149 (MQ=255)                     cccTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                         >  2:566478/1‑149 (MQ=255)                           tttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAatat                                                                                                                 >  1:122187/1‑151 (MQ=255)                           tttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAata                                                                                                                  <  1:578137/150‑1 (MQ=255)                             tttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAAtata                                                                                                                >  1:748630/1‑150 (MQ=255)                              ttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATaa                                                                                                               <  2:764028/150‑1 (MQ=255)                              ttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATaa                                                                                                               >  1:572987/1‑150 (MQ=255)                                aCCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATg                                                                                                            >  2:479431/1‑151 (MQ=255)                                 cccGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGt                                                                                                           <  1:680083/151‑1 (MQ=255)                                  ccGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGtt                                                                                                          <  2:248242/151‑1 (MQ=255)                                  ccGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGtt                                                                                                          >  1:613444/1‑151 (MQ=255)                                   cgcgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTa                                                                                                         <  1:364300/151‑1 (MQ=255)                                     cgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGCATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTAcc                                                                                                       >  2:587968/1‑151 (MQ=255)                                      ggTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTAccc                                                                                                      >  1:530103/1‑151 (MQ=255)                                       gTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCg                                                                                                     >  1:679716/1‑151 (MQ=255)                                         ccACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCgg                                                                                                    <  2:461719/150‑1 (MQ=255)                                            cacCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgc                                                                                                <  1:382993/151‑1 (MQ=255)                                              ccAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                              <  2:381158/151‑1 (MQ=255)                                                aGAAGCTCTTTGCCTGACCCTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcggcg                                                                                            <  2:343539/151‑1 (MQ=255)                                                aGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcggcg                                                                                            <  2:105097/151‑1 (MQ=255)                                                       ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGa                                                                                      >  1:451375/1‑150 (MQ=255)                                                        tttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGaa                                                                                     >  2:301906/1‑150 (MQ=255)                                                          tGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAAt                                                                                   <  1:394207/150‑1 (MQ=255)                                                            ccTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGGGATGAAATTg                                                                                 >  1:676837/1‑150 (MQ=255)                                                               gACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTgg                                                                             >  1:29392/1‑151 (MQ=255)                                                                acacTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGa                                                                            <  1:302520/151‑1 (MQ=255)                                                                 cacTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGAt                                                                           <  2:400254/151‑1 (MQ=255)                                                                    tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                                        >  2:699706/1‑151 (MQ=255)                                                                          aaGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTcc                                                                  >  2:364521/1‑151 (MQ=255)                                                                          aaGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTc                                                                   <  2:472932/150‑1 (MQ=255)                                                                              gcCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtt                                                              <  1:229872/151‑1 (MQ=255)                                                                                cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtttt                                                            >  1:316520/1‑151 (MQ=255)                                                                                 aCGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtttt                                                            >  1:134978/1‑150 (MQ=255)                                                                                       tCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCtt                                                     >  1:347506/1‑151 (MQ=255)                                                                                        cGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCtt                                                     <  2:419000/150‑1 (MQ=255)                                                                                          gAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgct                                                  >  1:538169/1‑151 (MQ=255)                                                                                          gAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgct                                                  <  2:535102/151‑1 (MQ=255)                                                                                           aaacaaacACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgctg                                                 <  1:336562/151‑1 (MQ=255)                                                                                            aacaaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTtgctgc                                                <  2:192194/151‑1 (MQ=255)                                                                                              caaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCgg                                              >  1:369799/1‑151 (MQ=255)                                                                                                aacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCt                                            <  1:364536/151‑1 (MQ=255)                                                                                                 acacCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTc                                           <  2:357509/151‑1 (MQ=255)                                                                                                 acacCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTc                                           >  2:65218/1‑151 (MQ=255)                                                                                                           cgtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGgtt                                  >  2:493363/1‑150 (MQ=255)                                                                                                                   gATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCa                         <  2:732356/151‑1 (MQ=255)                                                                                                                    aTACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCAt                        <  2:49589/151‑1 (MQ=255)                                                                                                                     tACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCAt                        >  1:372825/1‑150 (MQ=255)                                                                                                                          accaaccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTgg                   >  1:507489/1‑150 (MQ=255)                                                                                                                           ccaaccgGCGCCTCGACGAGGTGCATGCCTGGCGTTTGTTGGCCGAGAGGAATATAAATGTTTCCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACGCATTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTg                    <  1:82275/148‑1 (MQ=255)                                                                                                                           ccaaccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGAt                 <  2:575459/151‑1 (MQ=255)                                                                                                                            caaccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGCGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGGACGGTTTAGGCTGCCATAGTgggt                 >  2:580485/1‑148 (MQ=255)                                                                                                                               ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACg             <  2:393334/151‑1 (MQ=255)                                                                                                                                 aTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGa           <  1:289494/151‑1 (MQ=255)                                                                                                                                         gacgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCCTAGTGGATAACGGAACTGttt    >  1:73053/1‑150 (MQ=255)                                                                                                                                          acgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCa  <  1:730197/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                CTGAGCCGGGCTTCGCTTCACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGAACTGTTTCA  >  NC_000913/2333547‑2333834

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