Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	415,717	A→G	S13P (TCA→CCA)	sbcC ←	ATP dependent, structure specific DNA nuclease ‑ SbcC subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	415,717	0	A	G	100.0%	222.1 / NA	60	S13P (TCA→CCA)	sbcC	ATP dependent, structure specific DNA nuclease ‑ SbcC subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/35);  total (25/35)

GTTTCGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGAGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGA  >  NC_000913/415581‑415852                                                                                                                                         |                                                                                                                                        gTTTCGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGtt                                                                                                                                      <  1:1291896/140‑1 (MQ=255) gTTTCGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGtt                                                                                                                                      <  2:459627/140‑1 (MQ=255)   ttCGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCa                                                                                                                                    <  2:686862/140‑1 (MQ=255)     cGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAgg                                                                                                                                  <  2:193800/140‑1 (MQ=255)            cAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCa                                                                                                                           >  1:288052/1‑140 (MQ=255)               cgcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgc                                                                                                                        >  2:69917/1‑140 (MQ=255)               cgcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgc                                                                                                                        >  1:1340465/1‑140 (MQ=255)               cgcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgc                                                                                                                        >  1:1510879/1‑140 (MQ=255)               cgcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgc                                                                                                                        >  1:176296/1‑140 (MQ=255)                gcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGGAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgcg                                                                                                                       >  1:1292649/1‑140 (MQ=255)                gcCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGgcg                                                                                                                       <  1:69917/140‑1 (MQ=255)                    gACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAgg                                                                                                                   >  1:1175693/1‑140 (MQ=255)                    gACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAgg                                                                                                                   <  1:466761/140‑1 (MQ=255)                         aTGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTgag                                                                                                              <  2:949735/140‑1 (MQ=255)                               tCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAAtttt                                                                                                        >  2:80283/1‑140 (MQ=255)                               tCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAAtttt                                                                                                        >  1:661138/1‑140 (MQ=255)                                       ggtggtTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTc                                                                                                <  1:561876/140‑1 (MQ=255)                                         tggtTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCgt                                                                                              >  2:492683/1‑140 (MQ=255)                                           gtTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCgtgt                                                                                            >  2:327123/1‑140 (MQ=255)                                                  ccGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCgg                                                                                     <  1:865519/140‑1 (MQ=255)                                                  ccGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCgg                                                                                     <  1:619941/140‑1 (MQ=255)                                                   cGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGc                                                                                    >  1:778038/1‑140 (MQ=255)                                                   cGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGc                                                                                    <  1:395557/140‑1 (MQ=255)                                                    gCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCg                                                                                   <  1:1102327/140‑1 (MQ=255)                                                     cACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGa                                                                                  <  1:754679/140‑1 (MQ=255)                                                     cACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGa                                                                                  >  1:65080/1‑140 (MQ=255)                                                      aCCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAg                                                                                 <  1:1467283/140‑1 (MQ=255)                                                       ccTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAgg                                                                                <  2:176296/140‑1 (MQ=255)                                                         tGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGt                                                                              <  1:1359942/140‑1 (MQ=255)                                                                ccGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAAc                                                                       >  2:762218/1‑140 (MQ=255)                                                                  ggTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAAcgt                                                                     >  1:825693/1‑140 (MQ=255)                                                                   gTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCCATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAAcgtc                                                                    >  2:1438738/1‑140 (MQ=255)                                                                   gTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGAGGGTGAAATCAATGTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGAGAGGGTCTGCAAcgac                                                                    >  2:576484/1‑138 (MQ=255)                                                                    tAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAAcgtcg                                                                   >  1:853069/1‑140 (MQ=255)                                                                       tAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCgtgg                                                                <  2:1174854/140‑1 (MQ=255)                                                                            aaCAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGaaa                                                           <  2:1310880/140‑1 (MQ=255)                                                                             aCAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGaaaa                                                          >  1:791540/1‑140 (MQ=255)                                                                              cAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCTAGGGTATGCAACGTCGTGGTGAAAAg                                                         >  1:555027/1‑140 (MQ=255)                                                                              cAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAg                                                         <  2:5864/140‑1 (MQ=255)                                                                                gCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGAt                                                       >  2:790063/1‑140 (MQ=255)                                                                                   cGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGAtgct                                                    <  2:1175693/140‑1 (MQ=255)                                                                                    gTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGAtgctg                                                   <  2:1535793/140‑1 (MQ=255)                                                                                    gTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGAtgctg                                                   <  1:80283/140‑1 (MQ=255)                                                                                     ttGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGTAACGTCGTGGTGAAAAGAtgctgc                                                  <  1:630518/140‑1 (MQ=255)                                                                                     ttGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGAtgctgc                                                  >  2:454639/1‑140 (MQ=255)                                                                                      ggCTGGGGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCACGTGATTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGGTGCa                                                 <  2:698722/139‑1 (MQ=255)                                                                                         tGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGAc                                              >  1:896249/1‑140 (MQ=255)                                                                                          ggCGAACGGCTCGCGGGGGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACg                                             >  1:1158001/1‑140 (MQ=255)                                                                                                 ggCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTtgctgc                                      <  2:288052/140‑1 (MQ=255)                                                                                                  gCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTtgttgct                                     <  1:1996/140‑1 (MQ=255)                                                                                                  gCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTtgctgct                                     <  2:682/140‑1 (MQ=255)                                                                                                       cgGGTGAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGAt                                <  2:661138/140‑1 (MQ=255)                                                                                                            gAAATCAATCTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATc                           <  2:778038/140‑1 (MQ=255)                                                                                                                 caatcTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGtt                      <  1:886740/140‑1 (MQ=255)                                                                                                                     cTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                  <  2:1292649/140‑1 (MQ=255)                                                                                                                     cTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                  <  2:7196/140‑1 (MQ=255)                                                                                                                     cTTCCATTCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGttcttc                  <  2:896249/140‑1 (MQ=255)                                                                                                                           ttCGCCTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGc            <  2:65080/140‑1 (MQ=255)                                                                                                                                cTTTTAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGAc       >  1:479072/1‑140 (MQ=255)                                                                                                                                   ctAATGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGCATGCAACGTCGCGCTGAAAAGATGCTGCAGACTTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGc    >  2:696744/2‑140 (MQ=255)                                                                                                                                     aaTGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa  <  2:1510879/140‑1 (MQ=255)                                                                                                                                     aaTGGGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGa  <  1:454639/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        GTTTCGTGATACAGCGCCAGACAAATGGCGTCCAGCAGGGTGGTTTTCCCCGCACCTGTTGGGCCGGTAATAGCAAACAGCCCGTTGCTGGCGAACGGCTCGCGGGTGAAATCAATCTTCCATTCGCCTTTTAATGAGTTCAGGTTTTTCAGGCGCAGGCTGAGAATTTTCATGCTTCGTGTTCTCCGGCGAGGGTATGCAACGTCGTGGTGAAAAGATGCTGCAGACGTTGCTGCTGCGATTCATCCAGTTCTTCCAGTGCCAGACGGCGA  >  NC_000913/415581‑415852

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