Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,217,003	A→G	intergenic (+7/‑325)	ymgC → / → ycgG	Blue light, low temperature and stress induced protein/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	1,217,003	0	A	G	100.0%	243.1 / NA	63	intergenic (+7/‑325)	ymgC/ycgG	Blue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/34);  total (29/34)

CAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGC  >  NC_000913/1216856‑1217149                                                                                                                                                    |                                                                                                                                                   cAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGt                                                                                                                                                  >  2:206388/1‑150 (MQ=255)  aGGCTGTGCGCGTGCGTTTAAAAGATGAAACCAGTGTGATTCTTCCTGAGCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGGCTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTa                                                                                                                                                 <  2:879023/150‑1 (MQ=255)     cTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGt                                                                                                                                             >  2:611695/1‑151 (MQ=255)        tACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGttt                                                                                                                                           <  1:491083/150‑1 (MQ=255)          cGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGtttt                                                                                                                                          >  1:366303/1‑149 (MQ=255)           gCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGa                                                                                                                                       <  2:278060/151‑1 (MQ=255)           gCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGa                                                                                                                                       <  1:288504/151‑1 (MQ=255)            cTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGa                                                                                                                                       >  1:326419/1‑150 (MQ=255)             ttGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAAt                                                                                                                                     >  1:399008/1‑151 (MQ=255)              tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAAt                                                                                                                                     <  1:1651/150‑1 (MQ=255)                 gTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAAc                                                                                                                                 <  1:343253/151‑1 (MQ=255)                  tttAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAAc                                                                                                                                 >  2:410590/1‑150 (MQ=255)                     aaaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATc                                                                                                                              >  1:531859/1‑150 (MQ=255)                     aaaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATc                                                                                                                              <  1:490417/150‑1 (MQ=255)                     aaaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCa                                                                                                                             <  2:531809/151‑1 (MQ=255)                             aaaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATc                                                                                                                     >  1:109076/1‑151 (MQ=255)                              aaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATc                                                                                                                     >  2:605617/1‑150 (MQ=255)                              aaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATc                                                                                                                     >  1:569141/1‑150 (MQ=255)                                      tgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTg                                                                                                            >  2:116854/1‑151 (MQ=255)                                       gATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGa                                                                                                           <  1:270475/151‑1 (MQ=255)                                          tctttctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGc                                                                                                        >  1:647473/1‑151 (MQ=255)                                             ttctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCaa                                                                                                      >  2:800312/1‑150 (MQ=255)                                                     tctATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAt                                                                                             >  1:356331/1‑151 (MQ=255)                                                        aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAttt                                                                                           >  2:804302/1‑150 (MQ=255)                                                        aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAt                                                                                             <  1:557922/148‑1 (MQ=255)                                                        aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAt                                                                                             <  2:876642/148‑1 (MQ=255)                                                           tGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAgg                                                                                       >  1:817319/1‑151 (MQ=255)                                                            gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAg                                                                                        <  2:751372/149‑1 (MQ=255)                                                            gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGc                                                                                      <  1:435070/151‑1 (MQ=255)                                                            gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGc                                                                                      <  2:427093/151‑1 (MQ=255)                                                             cAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                                                     >  1:523437/1‑151 (MQ=255)                                                               aaTATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                                                   <  1:538266/151‑1 (MQ=255)                                                                atatTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                                                   >  2:171822/1‑150 (MQ=255)                                                                atatTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTa                                                                                  <  1:6100/151‑1 (MQ=255)                                                                    tGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                                              >  1:492632/1‑151 (MQ=255)                                                                     gCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                                              <  2:459540/150‑1 (MQ=255)                                                                               gAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATg                                                                   >  2:607060/1‑151 (MQ=255)                                                                                 aaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTc                                                                 <  2:148662/151‑1 (MQ=255)                                                                                 aaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTc                                                                 >  2:566174/1‑151 (MQ=255)                                                                                    gCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCg                                                               <  1:233404/150‑1 (MQ=255)                                                                                    gCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGc                                                              <  1:566224/151‑1 (MQ=255)                                                                                    gCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGc                                                              <  2:426402/151‑1 (MQ=255)                                                                                        aaGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAata                                                          <  2:324902/151‑1 (MQ=255)                                                                                        aaGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAat                                                           <  2:414824/150‑1 (MQ=255)                                                                                         aGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAatat                                                         <  1:835667/151‑1 (MQ=255)                                                                                           aCTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCa                                                       <  2:578129/151‑1 (MQ=255)                                                                                           aCTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCa                                                       <  1:842885/151‑1 (MQ=255)                                                                                              acacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgt                                                    <  1:489659/151‑1 (MQ=255)                                                                                                acTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                                                  >  1:153370/1‑151 (MQ=255)                                                                                                  tAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTc                                                 >  1:111734/1‑150 (MQ=255)                                                                                                   aaCGGATATCCTTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                                                >  2:681808/1‑150 (MQ=255)                                                                                                   aaCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                                                <  1:514287/150‑1 (MQ=255)                                                                                                       gATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAAt                                           >  2:557546/1‑151 (MQ=255)                                                                                                             attattAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgag                                       >  2:899532/1‑149 (MQ=255)                                                                                                                       aCAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct                           >  2:813212/1‑151 (MQ=255)                                                                                                                         aaggAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTGATTATGTCCGCAATATCAGGTGGAAAGAATGGAGAGTTCTCGCgctc                          <  2:821044/148‑1 (MQ=255)                                                                                                                            gAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAtt                      >  2:736152/1‑151 (MQ=255)                                                                                                                                    gCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACg               <  2:610848/150‑1 (MQ=255)                                                                                                                                         ctTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt            <  1:866280/148‑1 (MQ=255)                                                                                                                                         ctTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt            <  1:171846/148‑1 (MQ=255)                                                                                                                                           tAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGata        <  2:734259/150‑1 (MQ=255)                                                                                                                                            aGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGatat       <  1:718323/150‑1 (MQ=255)                                                                                                                                                cTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGc  >  2:196301/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGC  >  NC_000913/1216856‑1217149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑