Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,275,076	A→G	intergenic (‑35/‑395)	yjcB ← / → pdeC	uncharacterized protein YjcB/c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,275,076	0	A	G	100.0%	239.5 / NA	64	intergenic (‑35/‑395)	yjcB/pdeC	uncharacterized protein YjcB/c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/37);  total (27/37)

ACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAG  >  NC_000913/4274942‑4275203                                                                                                                                       |                                                                                                                                aCCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAgcg                                                                                                                            >  1:480424/1‑140 (MQ=255)  ccGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAgcgc                                                                                                                           >  2:249806/1‑140 (MQ=255)  ccGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGCGCAgcgc                                                                                                                           >  2:341610/1‑140 (MQ=255)     acgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCgcggc                                                                                                                        >  2:246341/1‑140 (MQ=255)     acgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCgcggc                                                                                                                        >  2:681852/1‑140 (MQ=255)     acgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCATCTGAAAAGAGCCGATGCCCTGGGCAGCgcggc                                                                                                                        >  2:844671/1‑140 (MQ=255)      cgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGTCGATGCCCTGGGCAGCgcggcg                                                                                                                       <  2:245595/140‑1 (MQ=255)      cgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCgcggcg                                                                                                                       <  2:914517/140‑1 (MQ=255)      cgaAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCgcggcg                                                                                                                       <  1:1005225/140‑1 (MQ=255)           aCGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCgttgtt                                                                                                                  <  1:681852/140‑1 (MQ=255)                 gTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACggg                                                                                                            <  2:875077/140‑1 (MQ=255)                   tGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAg                                                                                                          <  1:945647/140‑1 (MQ=255)                   tGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAg                                                                                                          <  2:802230/140‑1 (MQ=255)                   tGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAg                                                                                                          <  1:109043/140‑1 (MQ=255)                     aGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGc                                                                                                        <  1:715312/140‑1 (MQ=255)                      gtgtGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCg                                                                                                       <  2:480424/140‑1 (MQ=255)                           cTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTa                                                                                                  >  2:596572/1‑140 (MQ=255)                           cTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTa                                                                                                  <  1:329821/140‑1 (MQ=255)                           cTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTa                                                                                                  >  2:747221/1‑140 (MQ=255)                              gCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTacac                                                                                               >  2:1048291/1‑140 (MQ=255)                              gCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTacac                                                                                               >  1:951893/1‑140 (MQ=255)                              gCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTacac                                                                                               >  1:842124/1‑140 (MQ=255)                               ccAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACAct                                                                                              <  1:1048291/140‑1 (MQ=255)                                        tcatTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATc                                                                                     >  1:46892/1‑140 (MQ=255)                                          atTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATcgc                                                                                   <  2:493526/140‑1 (MQ=255)                                          atTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATcgc                                                                                   <  1:870169/140‑1 (MQ=255)                                           ttACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATcgcg                                                                                  <  1:249806/140‑1 (MQ=255)                                               ggCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGc                                                                              <  2:155678/140‑1 (MQ=255)                                                 cACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCtt                                                                            <  2:46892/140‑1 (MQ=255)                                                 cACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCtt                                                                            >  2:894500/1‑140 (MQ=255)                                                      agaagTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGa                                                                       >  2:785432/1‑140 (MQ=255)                                                             gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAacac                                                                   <  1:205595/137‑1 (MQ=255)                                                             gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAacac                                                                   >  2:205595/1‑137 (MQ=255)                                                             gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACaaa                                                                >  1:529243/1‑140 (MQ=255)                                                             gCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACaaa                                                                >  1:40593/1‑140 (MQ=255)                                                              ccAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACaaaa                                                               >  1:834615/1‑140 (MQ=255)                                                                aGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAAt                                                             >  2:1028181/1‑140 (MQ=255)                                                                  cGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGa                                                           >  2:951671/1‑140 (MQ=255)                                                                   gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGaa                                                          <  1:246341/140‑1 (MQ=255)                                                                   gaagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGaa                                                          <  2:40593/140‑1 (MQ=255)                                                                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGaaa                                                         <  1:341255/140‑1 (MQ=255)                                                                    aagaagaaCCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGaaa                                                         <  2:529243/140‑1 (MQ=255)                                                                              aCGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                               <  2:951893/140‑1 (MQ=255)                                                                              aCGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgttt                                               >  2:439286/1‑140 (MQ=255)                                                                                gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCtgtttgt                                             <  1:414417/140‑1 (MQ=255)                                                                                      ggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCAtt                                       >  1:643720/1‑140 (MQ=255)                                                                                       gtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTa                                      <  1:785432/140‑1 (MQ=255)                                                                                       gtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTa                                      <  2:894019/140‑1 (MQ=255)                                                                                       gtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTa                                      >  1:31050/1‑140 (MQ=255)                                                                                        tgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTaa                                     <  2:842124/140‑1 (MQ=255)                                                                                               aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                              >  2:634369/1‑140 (MQ=255)                                                                                               aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                              >  1:233707/1‑140 (MQ=255)                                                                                               aGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTtgtg                              <  1:128993/140‑1 (MQ=255)                                                                                                gCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGa                             >  2:82548/1‑140 (MQ=255)                                                                                                 ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGaa                            <  1:439286/140‑1 (MQ=255)                                                                                                 ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGaa                            <  1:4860/140‑1 (MQ=255)                                                                                                 ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGaa                            <  1:634369/140‑1 (MQ=255)                                                                                                    tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACta                         <  2:233707/140‑1 (MQ=255)                                                                                                        aCAATTCCTCTGAAAAGAGCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATc                     <  1:124067/140‑1 (MQ=255)                                                                                                                       agagCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtttg      <  2:371842/140‑1 (MQ=255)                                                                                                                        gagCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCCCAATTGAttgtttgt     <  2:310710/140‑1 (MQ=255)                                                                                                                         agCCGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt    >  2:972136/1‑140 (MQ=255)                                                                                                                           ccGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtttg      <  2:786905/136‑1 (MQ=255)                                                                                                                           ccGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgtttg      >  1:786905/1‑136 (MQ=255)                                                                                                                           ccGATGCCCTGGGCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAg  <  2:31050/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                ACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACGGCACTAAGAAGTTGCCAGCGAAGAAGAACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAG  >  NC_000913/4274942‑4275203

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