Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,566,171	C→T	G324S (GGT→AGT)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,566,171	0	C	T	100.0%	219.0 / NA	61	G324S (GGT→AGT)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/27);  total (34/27)

TCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATT  >  NC_000913/1566038‑1566307                                                                                                                                      |                                                                                                                                         tCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTActggtccg                                                                                                                                    >  1:744053/1‑140 (MQ=255)     tactGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTActggtccgg                                                                                                                                   >  1:370981/1‑137 (MQ=255)     tactGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTActggtccgg                                                                                                                                   <  2:370981/137‑1 (MQ=255)       ctGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGc                                                                                                                              >  2:486041/1‑140 (MQ=255)       ctGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGc                                                                                                                              >  2:842235/1‑140 (MQ=255)          cGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTActggtccgg                                                                                                                                   <  2:225756/132‑1 (MQ=255)          cGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTActggtccgg                                                                                                                                   >  1:225756/1‑132 (MQ=255)          cGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATg                                                                                                                           >  2:550221/1‑140 (MQ=255)            gACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATggg                                                                                                                         >  1:407694/1‑140 (MQ=255)               gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTg                                                                                                                                >  2:728955/1‑130 (MQ=255)               gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTg                                                                                                                                <  1:728955/130‑1 (MQ=255)               gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCaa                                                                                                                      >  1:404010/1‑140 (MQ=255)               gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCaa                                                                                                                      >  1:443198/1‑140 (MQ=255)               gTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCaa                                                                                                                      >  2:287539/1‑140 (MQ=255)                       aaaaGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGt                                                                                                              <  1:550221/140‑1 (MQ=255)                         aaGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGttt                                                                                                            >  1:132173/1‑140 (MQ=255)                              cTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGa                                                                                                       <  2:404010/140‑1 (MQ=255)                                gagagaCTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATa                                                                                                     >  2:51220/1‑140 (MQ=255)                                       tttACGCAGAATTTCATCACCAGTGTTATGGCCCCACATATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTa                                                                                              >  1:188475/1‑140 (MQ=255)                                              aGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcg                                                                                       <  1:740910/140‑1 (MQ=255)                                               gAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcgg                                                                                      <  2:132173/140‑1 (MQ=255)                                                aaTTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcggc                                                                                     >  1:643859/1‑140 (MQ=255)                                                aaTTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcggc                                                                                     >  2:843898/1‑140 (MQ=255)                                                aaTTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcggc                                                                                     >  1:186965/1‑140 (MQ=255)                                                 aTTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAAcggcg                                                                                    <  1:461736/140‑1 (MQ=255)                                                          ccAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTaa                                                                           >  1:582783/1‑140 (MQ=255)                                                             gtgtTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAttt                                                                        >  1:88971/1‑140 (MQ=255)                                                                ttATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGt                                                                     <  1:51220/140‑1 (MQ=255)                                                                ttATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGt                                                                     <  2:407694/140‑1 (MQ=255)                                                                   tGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGt                                                                     >  2:738733/1‑137 (MQ=255)                                                                   tGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGt                                                                     <  1:738733/137‑1 (MQ=255)                                                                   tGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAg                                                                  <  2:643859/140‑1 (MQ=255)                                                                     gCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTa                                                                <  2:88971/140‑1 (MQ=255)                                                                        ccACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACAt                                                             >  1:281277/1‑140 (MQ=255)                                                                        ccACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTTTCAGGACAt                                                             >  1:378686/1‑140 (MQ=255)                                                                          aCGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGAAATTTTGTCAGTACATcc                                                           >  2:255803/1‑140 (MQ=255)                                                                          aCGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATcc                                                           >  1:508426/1‑140 (MQ=255)                                                                           cGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGAAAATTTGTCAGTACATccc                                                          >  1:686553/1‑139 (MQ=255)                                                                           cGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCa                                                          >  2:610492/1‑140 (MQ=255)                                                                                   tGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATCTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTCAGTACATCCCTACCCAcc                                                  >  2:389844/1‑139 (MQ=255)                                                                                          tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGAc                                                   <  1:667607/132‑1 (MQ=255)                                                                                          tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGAc                                                   >  2:667607/1‑132 (MQ=255)                                                                                          tttGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGgc                                           >  1:841017/1‑140 (MQ=255)                                                                                            tGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGgcgc                                         <  2:744053/140‑1 (MQ=255)                                                                                                 ttATCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATAc                                    <  1:888490/140‑1 (MQ=255)                                                                                                    tCAACGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGGAAATTCGTCAGTACCTTCATAACGACATTTTGGGGCGATAattt                                 >  1:125148/1‑139 (MQ=255)                                                                                                        cGTCAATAATCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTCCGTACATCCATACCCACTTTGTGGGGCGATActtcttta                             >  2:719203/1‑138 (MQ=255)                                                                                                            aataatCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACa                         <  2:508426/140‑1 (MQ=255)                                                                                                              taatCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTAAATACATACCGACTTCGTGGCGCGGTACTTCTTCAAACaaa                       >  2:406679/1‑139 (MQ=255)                                                                                                              taatCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAAt                       >  1:511224/1‑140 (MQ=255)                                                                                                                atCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTc                     <  1:486041/140‑1 (MQ=255)                                                                                                                 tCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCa                    <  2:705972/140‑1 (MQ=255)                                                                                                                 tCAGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCa                    <  2:841017/140‑1 (MQ=255)                                                                                                                   aGCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACg                  <  1:828721/140‑1 (MQ=255)                                                                                                                    gCACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATCTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgc                 <  1:389844/140‑1 (MQ=255)                                                                                                                     cACTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATCCATACCGACTTCGTGGCGCGGTACTTCTTTAAACAAATCACgca                >  2:828721/1‑140 (MQ=255)                                                                                                                      aCTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcag               <  2:511224/140‑1 (MQ=255)                                                                                                                      aCTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcag               <  2:186965/140‑1 (MQ=255)                                                                                                                       cTGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcagc              <  1:719203/140‑1 (MQ=255)                                                                                                                       aggACAGTGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcagc              <  1:543173/138‑1 (MQ=255)                                                                                                                        tGACAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcagca             <  2:582783/140‑1 (MQ=255)                                                                                                                          aCAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGGTACTTCTTCAAACAAATTACGCAGCAAg           >  2:617333/1‑140 (MQ=255)                                                                                                                          aCAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGCTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAg           >  1:365915/1‑140 (MQ=255)                                                                                                                          aCAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGAGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcag               <  1:478258/136‑1 (MQ=255)                                                                                                                          aCAGCGGTGTACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGAGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACgcag               >  2:478258/1‑136 (MQ=255)                                                                                                                                   tACTGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAAtt  <  1:842235/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         TCACTACTGCGGACGTTGTCATAAAAGGCCTGAGAGACTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATCCATACCGACTTCGTGGCGCGATACTTCTTCAAACAATTCACGCAGCAAGGTAATAATT  >  NC_000913/1566038‑1566307

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