Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,872,148	T→C	S36S (AGT→AGC)	dgcJ →	putative diguanylate cyclase DgcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,872,148	0	T	C	100.0%	232.7 / NA	63	S36S (AGT→AGC)	dgcJ	putative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/40);  total (23/40)

TTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGG  >  NC_000913/1872013‑1872280                                                                                                                                        |                                                                                                                                     ttCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCt                                                                                                                                  >  2:492335/1‑140 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  2:315832/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  2:1136023/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  2:157341/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  1:492335/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  1:1309345/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  1:471624/140‑1 (MQ=255)   cACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTAt                                                                                                                                <  2:585118/140‑1 (MQ=255)        aCCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTgcg                                                                                                                           <  1:679453/140‑1 (MQ=255)        aCCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTgcg                                                                                                                           >  1:884310/1‑140 (MQ=255)              ggAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTatat                                                                                                                     >  1:361398/1‑140 (MQ=255)                  tGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTg                                                                                                                 >  2:664330/1‑140 (MQ=255)                  tGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTg                                                                                                                 >  2:66574/1‑140 (MQ=255)                  tGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTg                                                                                                                 >  2:850565/1‑140 (MQ=255)                  tGCCATTTATCACGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTg                                                                                                                 >  1:884915/1‑140 (MQ=255)                    ccATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCa                                                                                                               <  2:1065466/140‑1 (MQ=255)                       tttATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGaa                                                                                                            >  2:48047/1‑140 (MQ=255)                       tttATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGaa                                                                                                            >  1:95438/1‑140 (MQ=255)                          aTCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGc                                                                                                         >  2:1250886/1‑140 (MQ=255)                           tCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCa                                                                                                        <  1:412093/140‑1 (MQ=255)                             aTGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGa                                                                                                      <  1:1145304/140‑1 (MQ=255)                             aTGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGa                                                                                                      <  2:1320400/140‑1 (MQ=255)                             aTGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGa                                                                                                      <  1:624155/140‑1 (MQ=255)                             aTGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGa                                                                                                      <  1:226494/140‑1 (MQ=255)                              tGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGAc                                                                                                     >  1:523991/1‑140 (MQ=255)                                aaaTTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGActc                                                                                                   <  1:285747/140‑1 (MQ=255)                                 aaTTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGActcc                                                                                                  >  2:329229/1‑140 (MQ=255)                                  aTTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGActcct                                                                                                 <  2:451647/140‑1 (MQ=255)                                     gCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTcctcct                                                                                              <  2:361398/140‑1 (MQ=255)                                           tAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCctttcttt                                                                                        <  1:711612/140‑1 (MQ=255)                                           tAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCctttcttt                                                                                        >  1:1200096/1‑140 (MQ=255)                                           tAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTACTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGCCAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCa                                                                                                        <  2:1217301/124‑1 (MQ=255)                                           tAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTACTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCa                                                                                                        >  1:1217301/1‑124 (MQ=255)                                                tGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGAt                                                                                   <  2:95438/140‑1 (MQ=255)                                                 gCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATa                                                                                  >  2:72374/1‑140 (MQ=255)                                                     cGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTa                                                                              >  2:778362/1‑140 (MQ=255)                                                       gCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatc                                                                            <  1:1250886/140‑1 (MQ=255)                                                       gCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatc                                                                            <  1:733410/140‑1 (MQ=255)                                                        cATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatca                                                                           <  2:1315600/140‑1 (MQ=255)                                                          ttttATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATACCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatcaga                                                                         >  1:746047/1‑140 (MQ=255)                                                            ttATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTGCTTTCTTTATGATAAGTatcagaat                                                                       <  1:1252416/140‑1 (MQ=255)                                                            ttATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatcagaat                                                                       <  1:329229/140‑1 (MQ=255)                                                            ttATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTatcagaat                                                                       <  1:48047/140‑1 (MQ=255)                                                             tATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATaa                                                                                 <  1:853434/129‑1 (MQ=255)                                                             tATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATaa                                                                                 >  2:853434/1‑129 (MQ=255)                                                                 gccgcAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCagag                                                                  <  1:77032/140‑1 (MQ=255)                                                                 gccgcAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCagag                                                                  <  1:437904/140‑1 (MQ=255)                                                                  ccgcAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGt                                                                 >  2:935365/1‑140 (MQ=255)                                                                       aGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGc                                                            >  1:534848/1‑140 (MQ=255)                                                                         tgtCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCcg                                                          <  2:523991/140‑1 (MQ=255)                                                                         tgtCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCcg                                                          <  1:351034/140‑1 (MQ=255)                                                                                gtgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCAtgt                                                   <  2:1181425/140‑1 (MQ=255)                                                                                 tgCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCAtgtg                                                  <  2:884915/140‑1 (MQ=255)                                                                                    tGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATg                                               <  1:850565/140‑1 (MQ=255)                                                                                     gACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATACCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgc                                              <  2:746047/140‑1 (MQ=255)                                                                                      aCATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcg                                             <  2:1200096/140‑1 (MQ=255)                                                                                        aTCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATgcgcg                                           >  2:560516/1‑140 (MQ=255)                                                                                           ttccttccTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCcg                                                          >  1:1109656/1‑122 (MQ=255)                                                                                           ttccttccTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCcg                                                          <  2:1109656/122‑1 (MQ=255)                                                                                                tccttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCgct                                   >  2:566974/1‑140 (MQ=255)                                                                                                   ttATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCgctgct                                <  1:566974/140‑1 (MQ=255)                                                                                                       ttttGAACTTGTCGCCAGCGACAGAGCAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGc                            <  1:560516/140‑1 (MQ=255)                                                                                                                                 cAATGAGCGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAgcgg  <  1:935365/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     TTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCGCCGCAAGTGTCATTGTGCTGACATCTTCCTTCCTTATTTTTGAACTTGTCGCCAGCGACAGAGCAATGAGTGCCTATCTGCGCTATATCGTGCAGAAAGCAGACTCCTCCTTTCTTTATGATAAGTATCAGAATCAGAGTATTGCCGCGCATGTGATGCGCGCTCTCGCTGCTGAGCAGTCGGAAGTGTCGCCAGAACAGCGG  >  NC_000913/1872013‑1872280

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