Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,472,037	A→G	pseudogene (1521/3497 nt)	ydbA →	pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,472,037	0	A	G	100.0%	291.9 / NA	81	pseudogene (1521/3497 nt)	ydbA	pseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/42);  total (39/42)

AACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCT  >  NC_000913/1471905‑1472172                                                                                                                                     |                                                                                                                                        aacaCCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAtt                                                                                                                                   <  2:1545408/139‑1 (MQ=255) aacaCCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAtt                                                                                                                                   <  2:1005629/139‑1 (MQ=255)     ccGGGATGCTGGAAGTTAATAATTATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGGTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAttcttc                                                                                                                               >  1:1177064/1‑139 (MQ=255)     ccGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAttcttc                                                                                                                               >  2:339770/1‑139 (MQ=255)      cGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTcc                                                                                                                              <  2:1596013/139‑1 (MQ=255)        ggATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCca                                                                                                                            >  2:508862/1‑139 (MQ=255)             cTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACg                                                                                                                       <  2:577097/139‑1 (MQ=255)             cTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACg                                                                                                                       <  1:508862/139‑1 (MQ=255)                gAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:1831704/139‑1 (MQ=255)                gAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  1:215228/139‑1 (MQ=255)                      aataataatTCAGCTTTCAATAACCGCGGCGAGTTTATTATTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                              >  2:917318/1‑139 (MQ=255)                                   cTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAtt                                                                                                                                   >  2:409077/1‑105 (MQ=255)                                   cTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAtt                                                                                                                                   <  1:409077/105‑1 (MQ=255)                                   cTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCgg                                                                                                 >  2:438606/1‑139 (MQ=255)                                        aaTAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATg                                                                                            <  1:339770/139‑1 (MQ=255)                                              cgcgGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGtt                                                                                      <  1:257168/139‑1 (MQ=255)                                                  gcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  1:1829538/105‑1 (MQ=255)                                                  gcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  1:1711941/1‑105 (MQ=255)                                                  gcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:1711941/105‑1 (MQ=255)                                                  gcgAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  2:1829538/1‑105 (MQ=255)                                                     aGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATAtg                                                                                                         <  1:22193/113‑1 (MQ=255)                                                     aGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATAtg                                                                                                         >  2:22193/1‑113 (MQ=255)                                                           ttCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACgg                                                                         <  2:1177064/139‑1 (MQ=255)                                                               tGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  1:1596998/1‑92 (MQ=255)                                                               tGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:1596998/92‑1 (MQ=255)                                                               tGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACg                                                                     <  2:1669520/139‑1 (MQ=255)                                                                gATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAttct                                                                                                                                 >  2:1114142/1‑78 (MQ=255)                                                                gATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAAttct                                                                                                                                 <  1:1114142/78‑1 (MQ=255)                                                                     cGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  1:859468/86‑1 (MQ=255)                                                                     cGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  2:859468/1‑86 (MQ=255)                                                                                  gtgCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGt                                                  <  2:815995/139‑1 (MQ=255)                                                                                   tgCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      <  1:324584/134‑1 (MQ=255)                                                                                   tgCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      >  2:324584/1‑134 (MQ=255)                                                                                     cacaTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGAt                                                                   >  2:535214/1‑119 (MQ=255)                                                                                     cacaTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGAt                                                                   <  1:535214/119‑1 (MQ=255)                                                                                        aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGtt                                                                                      >  2:311278/1‑97 (MQ=255)                                                                                        aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGtt                                                                                      <  1:311278/97‑1 (MQ=255)                                                                                        aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                            >  2:680241/1‑139 (MQ=255)                                                                                        aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                            >  1:288808/1‑139 (MQ=255)                                                                                        aTCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                            >  1:769975/1‑139 (MQ=255)                                                                                          cAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                              >  2:1586242/1‑71 (MQ=255)                                                                                          cAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTa                                                                                                              <  1:1586242/71‑1 (MQ=255)                                                                                             ccAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTc                                       >  2:299895/1‑139 (MQ=255)                                                                                             ccAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTc                                       <  2:288808/139‑1 (MQ=255)                                                                                             ccAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTc                                       <  1:680241/139‑1 (MQ=255)                                                                                                  ccGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGa                                                                            <  1:685284/97‑1 (MQ=255)                                                                                                  ccGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGa                                                                            >  2:685284/1‑97 (MQ=255)                                                                                                  ccGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      <  1:210020/119‑1 (MQ=255)                                                                                                  ccGGTACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      >  2:210020/1‑119 (MQ=255)                                                                                                      tACGCTTTATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCtaat                              >  1:58420/1‑139 (MQ=255)                                                                                                           tttATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      >  2:1817320/1‑110 (MQ=255)                                                                                                           tttATAATACCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      <  1:1817320/110‑1 (MQ=255)                                                                                                               taataCCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGc                     >  2:41899/1‑139 (MQ=255)                                                                                                                aataCCGGTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCg                    <  1:41899/139‑1 (MQ=255)                                                                                                                   aCCGGTCACATGAACATCGGTAATTCTTCCCACAAC‑GAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGAtt                                                                  >  1:1661873/1‑89 (MQ=255)                                                                                                                   aCCGGTCACATGAACATCGGTAATTCTTCCCACAAC‑GAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGAtt                                                                  <  2:1661873/89‑1 (MQ=255)                                                                                                                    ccGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  1:152315/1‑39 (MQ=255)                                                                                                                    ccGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:152315/39‑1 (MQ=38)                                                                                                                     cGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  1:836642/38‑1 (MQ=255)                                                                                                                     cGGTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  2:836642/1‑38 (MQ=255)                                                                                                                     cGGTCACATGAACATCGGTAATTCTTCCCACAAC‑GAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      >  2:1408357/1‑99 (MQ=255)                                                                                                                     cGGTCACATGAACATCGGTAATTCTTCCCACAAC‑GAGCCATTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCg                                                      <  1:1408357/99‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:355156/36‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  1:355156/1‑36 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    >  1:733401/1‑36 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAg                                                                                                                    <  2:733401/36‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGaa                                                                       <  2:1275642/81‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGaa                                                                       >  2:152978/1‑81 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGaa                                                                       <  1:152978/81‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGaa                                                                       >  1:1275642/1‑81 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGaa                                                    <  1:800081/100‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGaa                                                    >  2:800081/1‑100 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGt                                                  >  2:636499/1‑102 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGt                                                  <  1:636499/102‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                            >  2:67698/1‑108 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTgg                                            <  1:67698/108‑1 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGaa             >  1:1207733/1‑139 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGaa             >  2:1822398/1‑139 (MQ=255)                                                                                                                       gTCACATGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGaa             >  1:1374386/1‑139 (MQ=255)                                                                                                                            aTGAACATCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCct        <  2:763456/139‑1 (MQ=255)                                                                                                                                  aTCGGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTtcttct  <  1:299895/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        AACACCGGGATGCTGGAAGTTAATAATAATTCAGCTTTCAATAACCGCGGCGAGTTTATTCTTGATAACGACAAGAATGCTGTGCACATCAACCAGTCCGGTACGCTTTATAATACCGGTCACATGAACATCAGTAATTCTTCCCACAACGGAGCCGTTAATATGTGGGGCGGAAATGGTCGTTTTATCAATGACGGAACGATTGATGTTTCTGCGAAGTCACTGGTAGTCAGCGCTAATAATGCCGGCGATCAGAATGCCTTCTTCT  >  NC_000913/1471905‑1472172

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