Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,435	T→C	T395A (ACT→GCT)	ylaB ←	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,435	0	T	C	100.0%	282.3 / NA	78	T395A (ACT→GCT)	ylaB	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/42);  total (36/42)

CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATG  >  NC_000913/477300‑477549                                                                                                                                        |                                                                                                                   cGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAg                                                                                                                 >  2:1100426/1‑139 (MQ=255)        cAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgc                                                                                                          <  2:873486/139‑1 (MQ=255)        cAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAgc                                                                                                          <  2:848547/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  1:385286/139‑1 (MQ=255)             ccTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  1:480759/139‑1 (MQ=255)               tGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtct                                                                                                   >  2:105664/1‑139 (MQ=255)                gTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctg                                                                                                  >  1:859679/1‑139 (MQ=255)                 tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAtctgt                                                                                                 <  2:678852/139‑1 (MQ=255)                 tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:627705/135‑1 (MQ=255)                 tAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:627705/1‑135 (MQ=255)                       aCTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:1341196/129‑1 (MQ=255)                       aCTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:1341196/1‑129 (MQ=255)                        cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          <  2:859679/139‑1 (MQ=255)                        cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          <  1:798494/139‑1 (MQ=255)                        cTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          <  1:1528939/139‑1 (MQ=255)                           aaaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:652060/1‑125 (MQ=255)                           aaaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:652060/125‑1 (MQ=255)                            aaCTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAGATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAAc                                                                                      <  2:177231/139‑1 (MQ=255)                                tGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                  <  2:1667527/139‑1 (MQ=255)                                 gAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTgtt                                                                                 >  1:1406538/1‑137 (MQ=255)                                  aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:1585405/1‑118 (MQ=255)                                  aaTACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:1585405/118‑1 (MQ=255)                                       cccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgat                                                                           >  2:1389036/1‑139 (MQ=255)                                        ccGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtgatt                                                                          >  1:131884/1‑139 (MQ=255)                                         cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          <  1:871404/122‑1 (MQ=255)                                         cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          >  2:871404/1‑122 (MQ=255)                                         cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                  <  2:613862/130‑1 (MQ=255)                                         cGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTAAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                  >  1:613862/1‑130 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAg                                                                                                           <  1:249849/103‑1 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAg                                                                                                           >  2:249849/1‑103 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGa                                                                                                      >  2:1387479/1‑108 (MQ=255)                                           tACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGa                                                                                                      <  1:1387479/108‑1 (MQ=255)                                                 aaTCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:1077111/103‑1 (MQ=255)                                                 aaTCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:1077111/1‑103 (MQ=255)                                                  atcatcAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGTTAGTGATTGATCATGTCa                                                                <  2:1406538/139‑1 (MQ=255)                                                     atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGt                                                                  >  2:1617498/1‑134 (MQ=255)                                                     atcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGt                                                                  <  1:1617498/134‑1 (MQ=255)                                                      tcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          >  1:1258937/1‑109 (MQ=255)                                                      tcaAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGAt                                                                                          <  2:1258937/109‑1 (MQ=255)                                                                     tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  2:268183/1‑83 (MQ=255)                                                                     tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  1:268183/83‑1 (MQ=255)                                                                     tGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAttt                                             <  1:1100426/139‑1 (MQ=255)                                                                          cGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                        >  1:331099/1‑139 (MQ=255)                                                                            ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAt                                               <  2:1502925/130‑1 (MQ=255)                                                                            ccTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGAt                                               >  1:1502925/1‑130 (MQ=255)                                                                                      gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:55859/66‑1 (MQ=255)                                                                                      gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:55859/1‑66 (MQ=255)                                                                                      gCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                            <  2:688040/139‑1 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGTGTTCAAGCGCGATCTGTCTGGGATTAACCTGa                                                                                  <  1:992502/84‑1 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCa                                                                                                             >  2:127329/1‑57 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCa                                                                                                             <  1:127329/57‑1 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  1:1175758/65‑1 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  2:1175758/1‑65 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgg                           >  1:608483/1‑139 (MQ=255)                                                                                       cGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCAATCTGTCTGGGATTAACCTGa                                                                                  >  2:992502/1‑84 (MQ=255)                                                                                        gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  2:1146402/1‑64 (MQ=255)                                                                                        gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  1:1146402/64‑1 (MQ=255)                                                                                        gagaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggg                          <  1:1389036/139‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:1000363/63‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     <  2:1625930/63‑1 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:1625930/1‑63 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGAt                                                                                                     >  1:1000363/1‑63 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgggg                         >  1:756988/1‑139 (MQ=255)                                                                                         agaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACgggg                         >  1:63805/1‑139 (MQ=255)                                                                                          gaAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACggggg                        <  2:608483/139‑1 (MQ=255)                                                                                                     ggCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTg             >  2:612732/1‑139 (MQ=255)                                                                                                         cTGGTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCg         <  2:756988/139‑1 (MQ=255)                                                                                                           ggTTTTCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGaa       <  1:612732/139‑1 (MQ=255)                                                                                                               ttCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAatat   >  2:1442200/1‑139 (MQ=255)                                                                                                                tCGGATCGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATAtg  <  1:105664/139‑1 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                        <  1:48388/95‑1 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTcc                                        >  2:48388/1‑95 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                            >  2:1550358/1‑107 (MQ=255)                                                                                                                      cGGCAAACTCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACg                            <  1:1550358/107‑1 (MQ=255)                                                                                                                              tCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtg                                                                             >  2:1316957/1‑50 (MQ=255)                                                                                                                              tCGCGTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGtg                                                                             <  1:1316957/50‑1 (MQ=255)                                                                                                                                  gTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACgcag                                                           >  1:185109/1‑64 (MQ=255)                                                                                                                                  gTTCAGCGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACgcag                                                           <  2:185109/64‑1 (MQ=255)                                                                                                                                        |                                                                                                                   CGACATCCAGATCCTGTAAATAACTTAAACTTGAATACCCCGTACCAAAATCATCAAGATAAATTTCATGGCCCGCCTCCCGGTAGCGAGAAATTATCGGGGCGCTGGTTTTCGGATCGGCAAACTCGCGTTCAGTGAGTTCAAGCGCGATCTGTCTGGGATTAACCTGATAGTGATTGATCATGTCACGCAGCAATTGCGGGATTTTTTCCGAGGTGAGCACGGGGGATTCAAGATTGATCGAAATATG  >  NC_000913/477300‑477549

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: ‑