Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,516,689	T→C	100%	T382A (ACG→GCG)	yfeA ←	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,516,689	0	T	C	96.5%	263.8 / ‑0.3	86	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (45/36);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516593‑2516784                                                                                                 |                                                                                                gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              >  2:1475930/1‑100 (MQ=255) gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              >  2:187969/1‑100 (MQ=255) gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              >  2:2286422/1‑100 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  2:902983/1‑100 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  1:2170297/1‑100 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  1:1262996/1‑100 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            <  2:2710454/100‑1 (MQ=255)    gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                           >  2:1622631/1‑100 (MQ=255)     cccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaa                                                                                          <  1:1598353/100‑1 (MQ=255)       cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        >  2:2481634/1‑100 (MQ=255)        ttCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAtt                                                                                       <  1:2643983/100‑1 (MQ=255)         tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                      <  1:1485330/100‑1 (MQ=255)         tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                      <  1:2548445/100‑1 (MQ=255)           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTTCTGTTTCAAAAACAGCAGCGACCGGCAAACTGAg                                                                                    >  1:1288798/1‑100 (MQ=255)           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  2:2616635/1‑100 (MQ=255)            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTATCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  1:940586/1‑99 (MQ=255)            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTGCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAACTTGAGc                                                                                   >  1:982879/1‑100 (MQ=255)            gTTTCCGGACAGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaacag                                                                                      >  1:981015/1‑94 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                 <  1:1366695/100‑1 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                 <  2:611792/100‑1 (MQ=255)                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               <  1:976568/100‑1 (MQ=255)                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               >  1:1066816/1‑100 (MQ=255)                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              <  2:1698521/100‑1 (MQ=255)                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             >  2:1362937/1‑100 (MQ=255)                    cccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAg                                                                           <  1:1117815/100‑1 (MQ=255)                     ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGt                                                                          <  2:2402185/100‑1 (MQ=255)                        cTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                       >  1:2032829/1‑100 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   >  1:1882962/1‑100 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  2:1125478/100‑1 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  1:2357097/100‑1 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 >  2:1420093/1‑100 (MQ=255)                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                               <  1:1558023/100‑1 (MQ=255)                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                               >  1:2854605/1‑100 (MQ=255)                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                               <  1:2136433/100‑1 (MQ=255)                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  2:2607340/1‑100 (MQ=255)                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  1:1417579/1‑100 (MQ=255)                                 aaCAGCAATTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  2:1523011/1‑100 (MQ=255)                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            >  2:1220003/1‑100 (MQ=255)                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            >  2:145410/1‑100 (MQ=255)                                     gcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTaa                                                          <  1:2178572/100‑1 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  2:1896936/100‑1 (MQ=255)                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  2:2485400/1‑100 (MQ=255)                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  2:2273428/1‑100 (MQ=255)                                        gTTAACTACCCGGCAATAGATAAAACTTTTCGGTTTCCTGCATCGACGGCAGCAGCGCACGGCAAATAGAGCGAAAACAGTGAACGCGCATCATTAAGc                                                        >  1:2400762/1‑99 (MQ=255)                                         ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:1935886/100‑1 (MQ=255)                                         ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      >  2:1659345/1‑100 (MQ=255)                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    >  1:1711387/1‑100 (MQ=255)                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   >  2:2570238/1‑100 (MQ=255)                                              tACCCGGCAATTGATACAACTTTTCGTGTTCCTGCATCAACGGCAGCAGCACACGGCAACTTGAGCGGATACCGTCAACGCGCATCATTAAGCCGTaa                                                   >  1:2753626/1‑98 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  1:2825460/1‑100 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               <  1:187967/100‑1 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAGTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               <  1:178502/100‑1 (MQ=255)                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                             <  1:350710/100‑1 (MQ=255)                                                    gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           >  2:1865225/1‑100 (MQ=255)                                                      aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                         <  1:1869691/100‑1 (MQ=255)                                                         tGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaa                                      <  1:1221250/100‑1 (MQ=255)                                                           aTACCACTTTCCGTTCTCCTGCATCAACGGCAGCAGCGCCCGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaa                                      <  1:1046009/98‑1 (MQ=255)                                                           aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCACACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    >  1:1604058/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  2:1524044/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:666290/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:1604509/1‑100 (MQ=255)                                                            tACAACTTTTCGTATTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  1:1800289/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  2:1904263/1‑100 (MQ=255)                                                               aaCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCaa                                <  1:112155/100‑1 (MQ=255)                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACCGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  1:658941/100‑1 (MQ=255)                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              >  2:413172/1‑100 (MQ=255)                                                                     tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          <  2:1996677/100‑1 (MQ=255)                                                                     tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          <  1:1475927/100‑1 (MQ=255)                                                                      cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  2:2670050/100‑1 (MQ=255)                                                                       gTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        >  2:1724848/1‑100 (MQ=255)                                                                        ttttCCTACATCAACGGCAGCAGCGAAAGGCGGATTGAAAGGATACAGTGAGAGCGCTTCATTCGGCCGTAAAGATGGCTGATAATCTCTAGCTTATCa                        >  1:2319717/1‑99 (MQ=255)                                                                         tttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAAt                      >  2:1811607/1‑100 (MQ=255)                                                                          ttCCTGCATCAACGGCAGCAGCGCACGTCAAATTGAGCTGATACATTAAACGCGCATCATTAAGCCTTAATGACGACTCATAAACTCAATATTATCAAt                      <  1:1283669/99‑1 (MQ=255)                                                                          ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  2:1882970/100‑1 (MQ=255)                                                                            ccTGCATCAACGGCAGCAGCGCACGGCAACTTGAGCGGATACAGTGAACGCGAATCATTAAGCCGTAATGCCGACTCATAAACTCACGATTATCAATgc                    >  1:662752/1‑99 (MQ=255)                                                                               gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAATGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                <  1:2176507/100‑1 (MQ=255)                                                                               gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                <  2:2748979/100‑1 (MQ=255)                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  1:1869299/1‑100 (MQ=255)                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             >  1:464669/1‑100 (MQ=255)                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           <  1:413170/100‑1 (MQ=255)                                                                                     aCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagca          >  2:1303306/1‑100 (MQ=255)                                                                                      cGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaa         <  2:1579589/100‑1 (MQ=255)                                                                                       ggCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCACATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaa        >  1:2160867/1‑100 (MQ=255)                                                                                          agcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAc     <  1:2288698/100‑1 (MQ=255)                                                                                          agcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAc     <  2:175252/100‑1 (MQ=255)                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  2:1891500/100‑1 (MQ=255)                                                                                                 |                                                                                                GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516593‑2516784

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