Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,275,912	A→G	M148V (ATG→GTG)	yjcC →	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	4,275,912	0	A→G	100.0%	226.6	60	M148V (ATG→GTG)	yjcC	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (31/29):  ref base (0/0):  total (31/29)

AATGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAT  >  NC_000913/4275716‑4276106                                                                                                                                                                                                     |                                                                                                                                                                                                   aaTGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGc                                                                                                                                                                                                 <  1:662448/200‑1 (MQ=255)    gTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTTTGGCTATAAAATGACCTATGTGCAGc                                                                                                                                                                                              >  1:572131/1‑200 (MQ=255)                          cTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTAt                                                                                                                                                                        >  1:223358/1‑200 (MQ=255)                                      gAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTgg                                                                                                                                                                               >  1:310982/1‑181 (MQ=255)                                      gAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCtctt                                                                                                                                                            >  1:269188/1‑200 (MQ=255)                                      gAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCcctt                                                                                                                                                            >  1:43274/1‑200 (MQ=255)                                            aTCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGag                                                                                                                                                      >  1:22670/1‑200 (MQ=255)                                            aTCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGag                                                                                                                                                      >  1:291039/1‑200 (MQ=255)                                                   cccGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCGTTTTTCTCTGGCGATACAATGACCTATGTGc                                                                                                                                                                                                 >  2:373103/1‑150 (MQ=255)                                                   ccagTGATAACCATTTATTATGCGCATCACTGATAGCGCCTGTAAGCGGCTATACGATTGCAACGGCAGATTATAAGCGTGAACCTGACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGc                                                                                                                                                                                                 <  2:242345/147‑1 (MQ=255)                                                     cGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGAt                                                                                                                                             >  1:199920/1‑200 (MQ=255)                                                         aTAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCt                                                                                                                                         <  1:782839/200‑1 (MQ=255)                                                             ccATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGaaa                                                                                                                                                                                       >  2:413394/1‑150 (MQ=255)                                                             ccATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCtgatg                                                                                                                                     <  1:707921/200‑1 (MQ=255)                                                              cATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGTaact                                                                                                                                                                                      >  2:1773/1‑148 (MQ=255)                                                                 tttttATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCgg                                                                                                                                                                                            >  4:82753/1‑141 (MQ=255)                                                                     tATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTgg                                                                                                                                                                               <  2:310957/150‑1 (MQ=255)                                                                       tGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAg                                                                                                                                                                                               <  3:383498/132‑1 (MQ=255)                                                                           cATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTAtg                                                                                                                                                                                  >  4:260157/1‑141 (MQ=255)                                                                             tCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTAtgtt                                                                                                                                                                                >  4:160640/1‑140 (MQ=255)                                                                                cTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAAc                                                                                                                                                                    <  2:596252/150‑1 (MQ=255)                                                                                       cgcCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTgg                                                                                                                                                                               <  3:82748/132‑1 (MQ=255)                                                                                       cgcATGTAAACGGCTATACGATTGCACCGGACGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCctctct                                                                                                                                                             <  2:271579/150‑1 (MQ=255)                                                                                           tGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAAccc                                                                                                                                                                  <  4:448866/141‑1 (MQ=255)                                                                                            gTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      <  1:514821/200‑1 (MQ=255)                                                                                            gTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      >  1:353117/1‑200 (MQ=255)                                                                                            gTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTACGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCct                                                                                                                                                                 >  4:526396/1‑141 (MQ=255)                                                                                             tAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCC‑TTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGag                                                                                                                                                      <  2:167086/150‑1 (MQ=255)                                                                                                 cGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGagtga                                                                                                                                                   >  2:450555/1‑150 (MQ=255)                                                                                                 cGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      <  1:450606/195‑1 (MQ=255)                                                                                                     tataCGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTgg                                                                                                                                                        <  4:531765/141‑1 (MQ=255)                                                                                                     tataCGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTgg                                                                                                                                                        <  4:444886/141‑1 (MQ=255)                                                                                                         cGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTgg                                                                                                                                                        >  1:605583/1‑137 (MQ=255)                                                                                                         cGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTgg                                                                                                                                                        <  2:605523/137‑1 (MQ=255)                                                                                                                aCCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTGTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGag                                                                                                                                                      <  3:472993/132‑1 (MQ=255)                                                                                                                      cGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGAc                                                                                                                              <  2:572074/150‑1 (MQ=255)                                                                                                                       gATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGagtgaagtg                                                                                                                                               >  3:373185/1‑132 (MQ=255)                                                                                                                        aTTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCtgatg                                                                                                                                     >  4:83582/1‑141 (MQ=255)                                                                                                                         ttATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCa                                                                         >  1:146675/1‑200 (MQ=255)                                                                                                                         ttATAAGCGTGAACCTAACGTTTCTATCCATTATTACCGCGCTACGCCTTTTTTCTCTCGCTCTAATATCACCTATGTGCAGCGGGTAAATTATGTGTCGGTTATCACCCCTCTCTTCTTGATCGCAGTGAt                                                                                                                                             >  3:595148/1‑132 (MQ=255)                                                                                                                                  tGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATgggg                                                                                                                  <  2:712454/150‑1 (MQ=255)                                                                                                                                    aaCCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGc                                                                                                                                                                                              >  1:624121/1‑72 (MQ=255)                                                                                                                                    aaCCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGc                                                                                                                                                                                              <  2:624058/72‑1 (MQ=255)                                                                                                                                           cGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGa                                                                                                         >  2:737803/1‑150 (MQ=255)                                                                                                                                           cGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      <  1:737871/153‑1 (MQ=255)                                                                                                                                           cGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCaa                                                       >  1:356182/1‑200 (MQ=255)                                                                                                                                            gTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGAt                                                                                                        <  2:559402/150‑1 (MQ=255)                                                                                                                                            gTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGAt                                                                                                        <  2:223123/150‑1 (MQ=255)                                                                                                                                             tttCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGc                                                                                                                         <  3:260146/132‑1 (MQ=255)                                                                                                                                              ttCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      <  2:199908/150‑1 (MQ=255)                                                                                                                                                     tattatTACCGCGATACGACTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGAt                                                                                                        <  4:432495/141‑1 (MQ=255)                                                                                                                                                      attattACTGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTTTGATGACCCGACATTGCAATGGGgtg                                                                                                                >  3:356870/1‑132 (MQ=255)                                                                                                                                                      attattACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGa                                                                                              >  2:513646/1‑150 (MQ=255)                                                                                                                                                           taCCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTAt                                                                                                           <  3:619247/132‑1 (MQ=255)                                                                                                                                                                cgATACGCCTTTTTTCGCTGGCTATAAAATGACCTATGTGAAGCGGGGAAATTATGTGGCGGTTAGCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATAc                                                                                                      <  3:577630/132‑1 (MQ=255)                                                                                                                                                                   tACGCCTTTTTTCTCTGGCTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAAc                                                                                          <  4:619036/141‑1 (MQ=255)                                                                                                                                                                                     cTATAAAATGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTc                                                               >  2:91670/1‑150 (MQ=255)                                                                                                                                                                                         aaaaTGACCTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCTGACATTGTAATGGGGTGTGTATGATTCGGTGACGAAAACCTTTTTCTCGTTAAGCAATGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGCTTTAACAGTACAAAGaaa         >  1:18623/1‑200 (MQ=255)                                                                                                                                                                                                ccTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAag                                                                      >  3:56920/1‑132 (MQ=255)                                                                                                                                                                                                ccTATGTGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATAGCTAt  >  1:561621/1‑200 (MQ=255)                                                                                                                                                                                                     |                                                                                                                                                                                                   AATGTTGAATATTATTCGTGGCTATCTTTATATTAATGAATTGATCTATGCCCGTGATAACCATTTTTTATGCTCATCGCTGATAGCGCCTGTAAACGGCTATACGATTGCACCGGCCGATTATAAGCGTGAACCTAACGTTTCTATCTATTATTACCGCGATACGCCTTTTTTCTCTGGCTATAAAATGACCTATATGCAGCGGGGAAATTATGTGGCGGTTATCAACCCTCTCTTCTGGAGTGAAGTGATGTCTGATGACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTAT  >  NC_000913/4275716‑4276106

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.