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	92.5%	265.7 / 6.7	89	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base C (35/47);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516593‑2516784                                                                                                 |                                                                                                gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              <  1:1947744/100‑1 (MQ=255) gTCGCCCTTCCGTTTCCGGCCCGCGCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAca                                                                                              <  1:2336234/100‑2 (MQ=255)   cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  1:969181/1‑100 (MQ=255)     cccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCCAGTGATACAACTTTTCGTTTTCCTGCATCAACAGCAGCAGCGCACGGCaa                                                                                          >  1:2413511/1‑100 (MQ=255)      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         >  2:2309743/1‑100 (MQ=255)      ccTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCaaa                                                                                         >  2:1291827/1‑100 (MQ=255)       cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        >  1:878824/1‑100 (MQ=255)          ccGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGa                                                                                     <  2:422890/100‑1 (MQ=255)           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGa                                                                                     >  1:1044664/1‑99 (MQ=255)           cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  1:2190403/1‑100 (MQ=255)              ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                 <  1:1076282/100‑1 (MQ=255)                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGa                                                                                >  2:2160320/1‑99 (MQ=255)                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               >  2:2492871/1‑100 (MQ=255)                ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               <  2:602625/100‑1 (MQ=255)                 cGGCCCGCACAGCACTAACAGCAGTTCACTACACGGCAATTGATAAAACTTTTAGTTTTCCTGCATCAACGGCAGCCGCGCACGGCAAATTGCGCGGa                                                                                >  1:107820/1‑98 (MQ=255)                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCgg                                                                                 >  1:1857235/1‑96 (MQ=255)                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             >  2:1164351/1‑100 (MQ=255)                  ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             <  2:882785/100‑1 (MQ=255)                   gCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                            <  2:296238/100‑1 (MQ=255)                    cccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACa                                                                            <  1:2401713/99‑1 (MQ=255)                     ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGt                                                                          <  2:1270967/100‑1 (MQ=255)                     ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGt                                                                          <  2:1187305/100‑1 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   >  2:1092588/1‑100 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  1:2029732/100‑1 (MQ=255)                            gCACTAACAGCAGTTCACTACCCGGCAATTGATACAACCTTTCGTTTTCCTGCATCAATGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgcg                                                                   <  1:1408846/100‑1 (MQ=255)                             cACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACgcgc                                                                  <  1:275968/100‑1 (MQ=255)                              aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 <  1:2197965/100‑1 (MQ=255)                               cTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcat                                                                <  2:871545/100‑1 (MQ=255)                                tAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                               <  1:2205940/100‑1 (MQ=255)                                 aaCAGCAGTTCACTACCCGGCAATTTATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCtcatca                                                              <  1:1739713/100‑1 (MQ=255)                                 aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatc                                                               >  2:113160/1‑99 (MQ=255)                                  aCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatcat                                                             <  1:435711/100‑1 (MQ=255)                                   cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGAATACAGTGAACGCGCATCAtt                                                            >  2:571189/1‑100 (MQ=255)                                       agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  1:2347014/100‑1 (MQ=255)                                            aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   <  1:405439/100‑1 (MQ=255)                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  <  2:1013350/100‑1 (MQ=255)                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  <  2:1321684/100‑1 (MQ=255)                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  <  2:2475788/100‑1 (MQ=255)                                             cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  >  1:865582/1‑100 (MQ=255)                                               aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATga                                                >  2:1082908/1‑100 (MQ=255)                                               aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATga                                                >  1:2443270/1‑100 (MQ=255)                                                cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  2:2592802/1‑100 (MQ=255)                                                  cGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacga                                             <  1:2254011/100‑1 (MQ=255)                                                     cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                          <  2:1031666/100‑1 (MQ=255)                                                     cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                          >  2:1453638/1‑100 (MQ=255)                                                       aTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCAt                                        <  1:743275/100‑1 (MQ=255)                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                       <  2:2118251/100‑1 (MQ=255)                                                        ttGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATa                                       <  2:435857/100‑1 (MQ=255)                                                         tGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCAGATACAGTGAACGCGCATCATTAAGCCGTAATTACGACTCATaa                                      <  1:353324/100‑1 (MQ=255)                                                          gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                     >  1:2349221/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   >  2:101688/1‑100 (MQ=255)                                                            tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACt                                   <  1:709763/100‑1 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  1:260416/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  2:1878186/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCa                                 >  1:531038/1‑100 (MQ=255)                                                              cAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATAATTAAGCCGTAATGACGACTCATAAACTCa                                 >  1:2580733/1‑100 (MQ=255)                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  1:51616/100‑1 (MQ=255)                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  1:2381530/100‑1 (MQ=255)                                                                aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAg                               <  1:575664/100‑1 (MQ=255)                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  1:370262/100‑1 (MQ=255)                                                                 cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  2:2306006/100‑1 (MQ=255)                                                                   tttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAtt                            <  1:1883691/100‑1 (MQ=255)                                                                   tttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAtt                            >  1:733473/1‑100 (MQ=255)                                                                   tttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGAtt                            >  2:833175/1‑100 (MQ=255)                                                                     tCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          <  2:2491494/100‑1 (MQ=255)                                                                      cggtttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCTTTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  1:1387283/97‑1 (MQ=255)                                                                       gTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCa                        <  1:2510165/100‑1 (MQ=255)                                                                          ttCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  1:1082904/100‑1 (MQ=255)                                                                           tCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgc                    >  2:841625/1‑100 (MQ=255)                                                                           tCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgc                    <  2:2238809/100‑1 (MQ=255)                                                                            ccTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcg                   >  1:1532194/1‑100 (MQ=255)                                                                              tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 <  1:101687/100‑1 (MQ=255)                                                                              tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 <  1:2047790/100‑1 (MQ=255)                                                                               gCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcag                >  2:1178922/1‑100 (MQ=255)                                                                                  tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             >  2:1463900/1‑100 (MQ=255)                                                                                   cAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  2:810101/100‑1 (MQ=255)                                                                                    aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  1:1887680/1‑100 (MQ=255)                                                                                     aCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagca          >  2:1626813/1‑100 (MQ=255)                                                                                     aCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagca          <  2:1001610/100‑1 (MQ=255)                                                                                      cGGCAGCAGCGTACGGCAAATCGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggcag            <  1:2084905/97‑1 (MQ=255)                                                                                       ggCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATTCGCAGGCAGCaan        <  2:2517290/100‑2 (MQ=255)                                                                                         cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCCTTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaa       <  2:1532203/99‑1 (MQ=255)                                                                                         cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      >  1:389037/1‑100 (MQ=255)                                                                                          agcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAc     <  1:128272/100‑1 (MQ=255)                                                                                           gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    >  2:594697/1‑100 (MQ=255)                                                                                           gcagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAAct    <  1:640335/100‑1 (MQ=255)                                                                                            cagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActc   <  2:416993/100‑1 (MQ=255)                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  1:2586239/100‑1 (MQ=255)                                                                                             agcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  >  2:1543245/1‑100 (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: ‑