Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,118	+C	coding (367/396 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,117	1	.	C	97.8%	227.6 / ‑4.3	91	A123G (GCG→GGG)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  major base C (54/33);  minor base A (2/0);  total (58/33)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACC  >  NC_000913/3554043‑3554182                                                                             |                                                                   ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                                   >  1:3245206/1‑76 (MQ=255) ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                                   >  2:3180926/1‑76 (MQ=255)  cgaTTTGCCCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                  <  2:4167793/74‑2 (MQ=255)  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                  <  1:4030910/76‑2 (MQ=255)  cTATTTGCGCCAGCAGCGGCGCCAGCACCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                                  <  1:4139832/76‑2 (MQ=255)   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGt‑gccg‑ccc                                                                 >  2:1476032/1‑74 (MQ=255)        ccGCCAGCAGCGGCTCCAGCACCCGGCTGATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCAc                                                            <  1:1905043/75‑1 (MQ=255)        acGCCAGCAGCGGCTCGAGCGCCCGGTGGATAACATCCGCAGGCGGCTCCGACGTGTTATCGGT‑GCCGCC‑GCTCAc                                                            <  1:4424497/75‑1 (MQ=255)           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCCCACTTc                                                         >  2:1197616/1‑76 (MQ=255)            cagcggGGGCCCCAGCCCCCGGGGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCa                                                        <  2:1474502/76‑1 (MQ=255)            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCa                                                        <  2:1157933/76‑1 (MQ=255)            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCa                                                        >  2:4675287/1‑76 (MQ=255)            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCa                                                        >  1:2015550/1‑76 (MQ=255)              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAac                                                      <  1:1184339/76‑1 (MQ=255)               cagcCGCGCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCAGT‑GCCGCC‑GCTCACTTCAaca                                                     <  1:1521357/76‑1 (MQ=255)                agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCCacac                                                    >  1:5158571/1‑76 (MQ=255)                  gggCTCCGCCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgc                                                  <  1:2450077/75‑1 (MQ=255)                  cGGCCCCAGGCCCCGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCAACAcgc                                                  <  2:4756446/76‑1 (MQ=255)                   ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCCACAcgcg                                                 >  2:4031848/1‑76 (MQ=255)                   ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                 >  1:973893/1‑76 (MQ=255)                   ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                 >  1:4032361/1‑76 (MQ=255)                   ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACAcgcg                                                 >  1:2620736/1‑76 (MQ=255)                   ggCTCCAGCACCCGGCGGATAAAATCCGCAAGCGGCGCCGACGGGCTATCGGT‑GCCGCC‑GCTCTCTCCGACAcgcg                                                 >  1:430747/1‑76 (MQ=255)                     cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGaa                                               <  2:2604784/76‑1 (MQ=255)                      tCCAGCACCGGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAg                                              <  1:3246191/76‑1 (MQ=255)                      tCCAGCACCCAGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAg                                              >  2:2708615/1‑76 (MQ=255)                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAgg                                             >  1:3415682/1‑76 (MQ=255)                       ccAGCACCCGGCGGAGAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑GGCCG‑CAGCTCACCTCAACACGCGAAgg                                             >  1:765461/1‑76 (MQ=255)                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGt                                            >  1:2253356/1‑76 (MQ=255)                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTc                                           >  1:2635715/1‑76 (MQ=255)                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑CGCCGCC‑GCTCCCTTCAACACGCGAAggcc                                           >  2:5465629/1‑74 (MQ=255)                          gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTcc                                          <  1:3948417/76‑1 (MQ=255)                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCac                                        >  2:1476794/1‑75 (MQ=255)                             cgcgggggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATc                                       <  2:2668259/74‑1 (MQ=255)                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                                                 <  1:345115/50‑3 (MQ=255)                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccg                                                                 >  2:345115/1‑48 (MQ=255)                             cccGGCGGATAAAAGCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATc                                       <  1:837798/76‑1 (MQ=255)                              ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCCACACGCGAAGGTCCATCg                                      >  1:3360784/1‑76 (MQ=255)                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCgg                                     <  2:3754572/76‑1 (MQ=255)                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCATccg                                     >  2:2795281/1‑74 (MQ=255)                                ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGc                                    <  1:5001649/76‑1 (MQ=255)                                ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGCCCATCCGc                                    >  1:2817253/1‑76 (MQ=255)                                 gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGG‑GGCCGCC‑GCTCACTTCAACACGCGAAGGGCCATCggga                                   >  1:148181/1‑74 (MQ=255)                                 gcggATAAAAGCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCa                                   <  2:2708224/76‑1 (MQ=255)                                   ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                 >  1:1384032/1‑76 (MQ=255)                                     aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                               <  1:2343869/76‑1 (MQ=255)                                      tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGAC‑GCTCAACTCAACACGCGAAGGTCCATCAGCAAACca                              >  1:2532808/1‑76 (MQ=255)                                       aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                             <  1:3791634/76‑1 (MQ=255)                                       aaaaTCCGCAGGCGGAGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACATCAACACGCGAAGGACCAGCGGCAAACcac                             >  2:3689105/1‑76 (MQ=255)                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                            <  1:2063379/76‑1 (MQ=255)                                        aaaTCCGCAGGCGGCGCCGACAGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACccca                            >  2:3558371/1‑76 (MQ=255)                                         aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCCTCGGCAAACCAAAg                           >  1:5702219/1‑76 (MQ=255)                                         aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAACCCACAg                           >  2:3393108/1‑76 (MQ=255)                                             cgcaGGCGGCGCCGCCGGGTGATCGG‑GGCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                       <  1:1959890/73‑1 (MQ=255)                                             cgcaGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                       <  1:1951158/73‑1 (MQ=255)                                               caGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACATGCGAAGGTCCATCGGCAAACCACAGCGCggg                     <  1:2428698/75‑1 (MQ=255)                                                aGGCGGCGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                    >  1:1397698/1‑76 (MQ=255)                                                 gtcggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                   <  2:4925328/74‑1 (MQ=255)                                                 ggcggcGCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                   <  1:2708615/76‑1 (MQ=255)                                                 ggcggcGCCGACGGGTTATCGAT‑GCCGCC‑GCTCACTTCAACACCCGAAGCTCCACCGGTAAAACACAGCGCGGgca                   >  1:4380447/1‑76 (MQ=255)                                                  ggggcgCCGGCGGGTTAGCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                  <  2:4985708/74‑1 (MQ=255)                                                  gcggcgCCGACGGGTCATCGGT‑GCCGAC‑GCTCAATTAAACACGCGAAGGTCACTCGGCAAACCAAAGTACGGgcag                  >  1:282521/1‑76 (MQ=255)                                                   cggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagc                 <  1:5237032/76‑1 (MQ=255)                                                    ggcgCCGCCGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                <  1:2304804/76‑1 (MQ=255)                                                    ggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                >  1:2042131/1‑76 (MQ=255)                                                    ggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACACCGCCGgcagcc                >  2:4957997/1‑75 (MQ=255)                                                    ggcgCCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCCGCAAACCACAGCGCGGgccgca                >  2:1983737/1‑76 (MQ=255)                                                     gcgcCGAGGGGTTAGCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc               <  2:1914021/76‑1 (MQ=255)                                                     gcgcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGGCCATCGGCAAACCCCAGCGCGGGCAGCAc               >  2:859560/1‑76 (MQ=255)                                                      cgcCGGCGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc              <  1:610000/76‑1 (MQ=255)                                                      cgcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc              >  2:2038761/1‑76 (MQ=255)                                                      cgcCGACGGGGTATAGGT‑GCAG‑CAGCTCACTTCAACACGAGAAGGTACATCGGCAAACCACAGAGCGGGCCGCAcc              >  2:2475317/1‑76 (MQ=255)                                                       gcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg             >  1:4495255/1‑76 (MQ=255)                                                       gcCGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCAGCGGCAAACCACAGCGCGGGCAGCAccc             >  2:3914179/1‑75 (MQ=255)                                                        ccGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCTGCAAACCACAGCGCGGGCAGCACCGt            >  1:1060714/1‑76 (MQ=255)                                                        ccGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt            >  1:1179725/1‑76 (MQ=255)                                                         cgccGGGTTTTTGGGCGCCGCC‑GGTCACTGCAACACGCGAAGGTCCAGCGGCAAACCACAGCGCGGGCAGCCCCGt            <  2:581297/73‑1 (MQ=255)                                                         cGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCgct           >  2:2571927/1‑74 (MQ=255)                                                         cGACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt           >  1:4693144/1‑76 (MQ=255)                                                         cGACGGGGTATCCGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCCCCGtt           <  2:1281511/76‑1 (MQ=255)                                                          gACGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCCAACCACAGCGCGGGCACGACCGttc          >  2:5600180/1‑75 (MQ=255)                                                           aCGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTg         >  1:2047437/1‑76 (MQ=255)                                                           aCGGGTTAGCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTg         <  1:2779434/76‑1 (MQ=255)                                                           aCGGGGTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGGTTg         >  1:2284549/1‑76 (MQ=255)                                                            cGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgc        >  2:2719011/1‑76 (MQ=255)                                                            cGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgc        >  1:1363492/1‑76 (MQ=255)                                                            cGGGTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGCGCAGCCCCGTTTgc        >  1:58355/1‑76 (MQ=255)                                                             gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAACCCCCAGCGCGCGCAGCCCCGTTTgcg       >  2:104855/1‑75 (MQ=255)                                                             gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCCGCACCGTTTgca       >  2:4152974/1‑76 (MQ=255)                                                             gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca       >  2:2109738/1‑76 (MQ=255)                                                             gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca       >  1:375651/1‑76 (MQ=255)                                                             gggTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca       >  1:309657/1‑76 (MQ=255)                                                               gTTATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagc     >  1:5702860/1‑76 (MQ=255)                                                                ttATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca    <  1:4091847/76‑1 (MQ=255)                                                                ttATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca    >  2:4986593/1‑76 (MQ=255)                                                                 tATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCCGCACCGTTTGCAGCAc   >  2:4142438/1‑76 (MQ=255)                                                                 tATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc   <  1:2445/76‑1 (MQ=255)                                                                 tATCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAc   >  1:3982501/1‑76 (MQ=255)                                                                 tATCGGT‑GCCGCC‑GCTAACTTCAACACGCGAAGGTCCATCGGCAAACCAACGCGCGGGCAGCACCGTTCGCAGCAc   >  1:3998456/1‑76 (MQ=255)                                                                  aTCGGT‑GCCGCC‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGCTTGCAGCAcc  >  2:5609412/1‑76 (MQ=255)                                                                             |                                                                   CCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑C‑GCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACC  >  NC_000913/3554043‑3554182

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