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	98.9%	288.4 / ‑4.5	88	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (47/39);  total (49/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTT  >  NC_000913/2516593‑2516785                                                                                                 |                                                                                                 gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACACCTTTTCGTTTTCCTGCATCCCCGGCAGCAGCGCAc                                                                                                >  1:2570423/1‑99 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCCTCAACGGCAGCAGCGCACGGc                                                                                             <  2:2042755/100‑1 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                             <  2:307961/100‑1 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                             >  2:2447811/1‑100 (MQ=255)    gCCCTTCCTTTTCCGGCCCGCTCAGCACTATCAGCAGTGCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCATCGGCAGCAGCGCACGGCa                                                                                            <  2:2052765/100‑1 (MQ=255)    gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCCGCACTTGACACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                            >  1:867500/1‑100 (MQ=255)        ttCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAtt                                                                                        >  2:2293506/1‑100 (MQ=255)         tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTGACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                       <  1:1703829/100‑1 (MQ=255)         tCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTg                                                                                       <  1:2114099/100‑1 (MQ=255)          ccGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGa                                                                                      >  2:882400/1‑100 (MQ=255)           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                     >  2:1777389/1‑100 (MQ=255)            gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCCCTACCCGGCAAGTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                    >  2:2770674/1‑100 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                  >  2:1722164/1‑100 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                  <  1:984325/100‑1 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                  <  1:2233039/100‑1 (MQ=255)                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                               >  2:553396/1‑100 (MQ=255)                 cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                               >  1:831071/1‑100 (MQ=255)                  ggCCCACACAGCACTAACAGCAGCACACTACCCGGCAGTTGACACAGCTTTTCGTTTTCCTCCATCAACGCCACCAGCGCACGGAAAAtgaggcaga                                                                                 >  1:1626775/1‑89 (MQ=255)                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                          <  2:2132946/100‑1 (MQ=255)                      cGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTg                                                                          <  2:2378723/100‑1 (MQ=255)                       gCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGa                                                                         >  1:1779447/1‑100 (MQ=255)                       gCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGa                                                                         >  1:830690/1‑100 (MQ=255)                           aGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgc                                                                     <  2:325396/100‑1 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCCATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                  <  2:1043804/100‑1 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAGGGCAAATTGAGCGGATACAGTGAACGCGca                                                                  <  1:1169298/100‑1 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                  >  2:2576181/1‑100 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                  <  2:2741027/100‑1 (MQ=255)                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcat                                                                 >  2:480494/1‑100 (MQ=255)                                 aaCAGCATTTCACTACCCGGCAATTGATACAACTTTCCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                               <  1:1937377/100‑1 (MQ=255)                                  aCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatcat                                                              <  2:2261280/100‑1 (MQ=255)                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                             >  1:1258094/1‑100 (MQ=255)                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCACCAtt                                                             >  1:1201897/1‑100 (MQ=255)                                    agcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                            >  1:978996/1‑100 (MQ=255)                                      cagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAg                                                          >  1:306417/1‑100 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                         <  1:2651855/100‑1 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                         <  1:247995/100‑1 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                         <  2:2146358/100‑1 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                         <  1:458239/100‑1 (MQ=255)                                        gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                        >  2:2504337/1‑100 (MQ=255)                                         tgcACTACCCGGCAATTGATACAACTTTTCGTTTTCCTCAATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                       <  1:1205791/98‑1 (MQ=255)                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                     >  2:457476/1‑100 (MQ=255)                                           cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                     >  2:60427/1‑100 (MQ=255)                                           cACTACCCGACAATAGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCAGTa                                                     >  1:1862422/1‑100 (MQ=255)                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                    <  2:2147980/100‑1 (MQ=255)                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                   <  2:2534144/100‑1 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                                >  1:1512387/1‑100 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTACTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAGGCGCATCATTAAGCCGTAATgac                                                >  2:1676793/1‑100 (MQ=255)                                                 ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                               >  1:1899063/1‑100 (MQ=255)                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                              >  2:2239334/1‑100 (MQ=255)                                                   ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGTTACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                             >  2:2812339/1‑100 (MQ=255)                                                    gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                            <  2:324220/100‑1 (MQ=255)                                                     cAATTGATACAACTTGTCGTTTTCCCGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGTATCATTAAGCCTTAATGACGACt                                            <  1:1343522/99‑1 (MQ=255)                                                      aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCa                                          >  2:1510328/1‑100 (MQ=255)                                                       aTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCAt                                         >  2:2905160/1‑100 (MQ=255)                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                        >  2:2136245/1‑100 (MQ=255)                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                        <  2:2255129/100‑1 (MQ=255)                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                        <  1:2735413/100‑1 (MQ=255)                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                      >  2:2665653/1‑100 (MQ=255)                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                      >  2:2865942/1‑100 (MQ=255)                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                      >  2:941203/1‑100 (MQ=255)                                                           aTACAACTTTTCTTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                     >  1:1933094/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                     >  2:2405631/1‑99 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                  <  2:918274/100‑1 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                  >  2:1356800/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                  >  1:2832678/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                  >  1:40614/1‑100 (MQ=255)                                                               aaCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAATCaa                                 >  1:1919008/1‑100 (MQ=255)                                                                  ttttCGTTTTCCTGCATCAACGGCTGCAGCGCCCGGTAAATTGAGCGGTTACAGTGAACGCGCATCATTAAGCCGTAGTGACGACTCATAAACTCAAGa                               <  1:299465/99‑1 (MQ=255)                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                              <  1:1720853/100‑1 (MQ=255)                                                                  ttttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAt                              <  1:749725/100‑1 (MQ=255)                                                                   tttCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAtt                             >  2:1554927/1‑100 (MQ=255)                                                                    ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTa                            <  1:867777/100‑1 (MQ=255)                                                                     tCGTTTTCCTGCATCAACGGCAGCGGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                           <  1:1288325/100‑1 (MQ=255)                                                                     tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                           <  1:308496/100‑1 (MQ=255)                                                                      cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                           >  1:710807/1‑99 (MQ=255)                                                                      cGTTTTCCTGCATCAACGGCAGCAGAGCAGGGCAAATTGAGAGGATACAGTGAACGAGCAGCATTAAGCCGGAATGACGACTCATAAACCCAAGATTATc                          >  1:214883/1‑100 (MQ=255)                                                                       gTTTTCCTGCATTAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGTGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                         >  2:566598/1‑100 (MQ=255)                                                                         tttCCTGCATCAACCGGAGAAGCAAACGCCAAATTGACCGGATACAGCAAACACGCATCATTAAGCCGTAATAACGACTCATAAACTCAAGATTATCaa                        >  1:726998/1‑99 (MQ=255)                                                                          ttCCTGCATCAACGGCAGCAGCGCACTGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                      <  1:2864669/100‑1 (MQ=255)                                                                             cTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcgc                   >  2:1268411/1‑100 (MQ=255)                                                                              tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                  <  1:2125659/100‑1 (MQ=255)                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc               <  2:21055/100‑1 (MQ=255)                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc               <  2:1536048/100‑1 (MQ=255)                                                                                 aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACGCATAACCTCAAGATTATCAATGCgcaggc               >  1:2571398/1‑100 (MQ=255)                                                                                        gcagcaGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaa        <  1:2259770/100‑1 (MQ=255)                                                                                         cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa       >  1:1041306/1‑100 (MQ=255)                                                                                           gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct     >  2:1061912/1‑100 (MQ=255)                                                                                           gcagcGCACGGCAAATGGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct     <  1:1982013/100‑1 (MQ=255)                                                                                              gcgcACTGCAAATTGAGCGGATACAGTAAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCtt  <  1:770061/100‑1 (MQ=255)                                                                                                 |                                                                                                 GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCTT  >  NC_000913/2516593‑2516785

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