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	100.0%	306.9 / NA	91	T382A (ACG→GCG)	yfeA	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/45);  total (46/45)

AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784                                                                                                  |                                                                                                aGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATCCAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCAc                                                                                               >  2:1530881/1‑100 (MQ=255)  gTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACg                                                                                              <  1:1886170/100‑1 (MQ=255)    ngCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGc                                                                                            >  2:193047/2‑100 (MQ=255)    cGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACGACTTTTCGTTTTCCTGCATCAACGGCAGCACCGCACGGc                                                                                            >  2:2141653/1‑100 (MQ=255)     gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                           >  1:244892/1‑100 (MQ=255)     gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCa                                                                                           >  2:63484/1‑100 (MQ=255)     gCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTAGTTTTCATGAATCAACGGAAGCAGCGCACGGCa                                                                                           >  1:397283/1‑100 (MQ=255)        cTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAAt                                                                                        >  1:2140681/1‑100 (MQ=255)            cGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAg                                                                                    >  2:2152490/1‑100 (MQ=255)             gTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGc                                                                                   <  2:1163950/100‑1 (MQ=255)              tttCCGTCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  2:529274/100‑1 (MQ=255)              tgtCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  <  2:529273/98‑1 (MQ=255)               ttCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCgg                                                                                 >  1:2658151/1‑100 (MQ=255)                tCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGa                                                                                >  2:1873219/1‑100 (MQ=255)                 ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCg                                                                                  >  2:1832112/1‑97 (MQ=255)                 ccGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGAt                                                                               >  2:2620281/1‑100 (MQ=255)                  cGGCCCGCTCAGCACTAACAGCAGTTCATTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:2657128/1‑100 (MQ=255)                  cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  2:1732307/1‑100 (MQ=255)                  cGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATa                                                                              >  1:1406837/1‑100 (MQ=255)                   ggCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTCCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATAc                                                                             >  2:1867607/1‑100 (MQ=255)                      ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGt                                                                          <  1:123730/100‑1 (MQ=255)                      ccGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGt                                                                          >  1:628434/1‑100 (MQ=255)                         cTCAGCACTAACAGCAGTTCACTACCCGGCACTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGaa                                                                       >  1:482941/1‑100 (MQ=255)                          tCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAc                                                                      <  1:1513478/100‑1 (MQ=255)                           cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGTAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                     <  2:757510/100‑1 (MQ=255)                           cAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcg                                                                     <  2:500582/100‑1 (MQ=255)                            aGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAAcgc                                                                    <  2:92602/100‑1 (MQ=255)                              cACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACgcgc                                                                  >  2:2166681/1‑100 (MQ=255)                               aCTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGca                                                                 >  1:70299/1‑100 (MQ=255)                                  aaCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGcatca                                                              >  1:2160077/1‑100 (MQ=255)                                   aCAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATTAACGGCAGCAGCGCACGGCAAATGGAGCGGATACAGTGAACGCGcatcat                                                             <  1:693257/100‑1 (MQ=255)                                    cagcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCAtt                                                            <  2:2476522/100‑1 (MQ=255)                                     agcagTTCACTACCCGGCAATTGATACAACTTTTCTTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                           <  1:2391261/100‑1 (MQ=255)                                     agcagTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTa                                                           >  2:1837302/1‑100 (MQ=255)                                        agTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGc                                                        <  2:2140655/100‑1 (MQ=255)                                         gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAATTTGAGCGCACACAGTGAACACGCAACATTAAGcc                                                       >  2:899890/1‑100 (MQ=255)                                         gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  2:2286879/1‑100 (MQ=255)                                         gTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGcc                                                       >  1:1361836/1‑100 (MQ=255)                                          ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:948955/100‑1 (MQ=255)                                          ttCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCg                                                      <  2:1753634/100‑1 (MQ=255)                                            cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    >  2:2419435/1‑100 (MQ=255)                                            cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    >  1:2195068/1‑100 (MQ=255)                                            cACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTa                                                    >  1:882222/1‑100 (MQ=255)                                             aCTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTaa                                                   <  1:1842421/100‑1 (MQ=255)                                              cTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAAt                                                  >  1:182076/1‑100 (MQ=255)                                               tACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATg                                                 <  2:744387/100‑1 (MQ=255)                                               tACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATg                                                 <  2:1764240/100‑1 (MQ=255)                                                aCCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATga                                                >  1:1838225/1‑100 (MQ=255)                                                 cccGGCAATTGATACAACTTTTCGTTTTCCTGCATTAACGGCAGCAACGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               <  1:1704378/100‑1 (MQ=255)                                                 cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               <  1:304937/100‑1 (MQ=255)                                                 cccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgac                                               >  2:1881631/1‑100 (MQ=255)                                                  ccGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacg                                              <  2:576606/100‑1 (MQ=255)                                                    ggCAATTGATACGACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  1:2141679/100‑1 (MQ=255)                                                    ggCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATgacgac                                            <  1:2313774/100‑1 (MQ=255)                                                     gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  1:63484/100‑1 (MQ=255)                                                     gCAATTGATACAACTTTTCGTTTTCCTGCATCAACGACAGCAGCTCACGGCAAATTGAGCTGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACt                                           <  2:1575387/100‑1 (MQ=255)                                                      cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                          >  2:2053240/1‑100 (MQ=255)                                                      cAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                          >  2:2095798/1‑100 (MQ=255)                                                       aaTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTc                                          <  1:2619319/99‑1 (MQ=255)                                                        aTTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCAt                                        >  2:1451855/1‑100 (MQ=255)                                                         gttgtACAACTCTCCGTTTTCCTGGAGCAACGGCAGCAGCTCACGGCAAATTGAGAGGATACACTGAACGCGACACATTACGCCGTAATGACGACTCaaa                                       <  2:485444/96‑3 (MQ=255)                                                           gATACAACTTTTCGTTTTCCTGCATCAACGGCAGCATCTCACGGTAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaa                                     <  1:747504/100‑1 (MQ=255)                                                            aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  2:735455/100‑1 (MQ=255)                                                            aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  1:699053/100‑1 (MQ=255)                                                            aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  1:1983631/100‑1 (MQ=255)                                                            aTACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAAc                                    <  2:567493/100‑1 (MQ=255)                                                            aTACAACTTTTCGTTTTCCTGCATCAACGGCAGAAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGAAGACTCATAAAc                                    >  2:2250171/1‑100 (MQ=255)                                                             tACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATaaaat                                   >  1:2214167/1‑98 (MQ=255)                                                                aaCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCaa                                <  1:1873236/100‑1 (MQ=255)                                                                 aCTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAACGGATACAGTGAACGCGCATCATAAAGCCGTAATGACGACTCAAAAACTCAAg                               >  2:1464730/1‑100 (MQ=255)                                                                  cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  2:2274114/100‑1 (MQ=255)                                                                  cTTTTCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGa                              <  2:1389799/100‑1 (MQ=255)                                                                     ttCGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCGTAAACTCAAGATTa                           <  1:2008287/100‑1 (MQ=255)                                                                       cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCTTTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         <  2:1416030/100‑1 (MQ=255)                                                                       cGTTTTCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTAt                          <  1:67561/99‑1 (MQ=255)                                                                       cGTTTTCCTGCATCAACGGCAGAAGCGCAAGGCAAATTGAGCGGATACAGTGAACGCGAATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATc                         >  2:2150115/1‑100 (MQ=255)                                                                        gtcttCCTGCATCAACGGCAGCAGCGCACTGAAAATTGAGCGGATACAGTGAACGTTCATCATTATGCCGTAATGACGACTCATAAACTCAAGATTATCa                        <  2:2385011/97‑1 (MQ=255)                                                                           ttCCTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATg                     <  1:1560766/100‑1 (MQ=255)                                                                             ccTGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATgcg                   >  2:2397837/1‑100 (MQ=255)                                                                               tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 <  2:723006/100‑1 (MQ=255)                                                                               tGCATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgca                 >  2:2076542/1‑100 (MQ=255)                                                                                 cATCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcagg               >  2:779104/1‑100 (MQ=255)                                                                                  aTCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggc              >  1:2066075/1‑100 (MQ=255)                                                                                   tCAACGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCgcaggca             >  2:1522344/1‑100 (MQ=255)                                                                                     aaCGGCAGCAGCGCACGGCAAATTGAGCGGATACAGTTAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagc           >  2:1413182/1‑100 (MQ=255)                                                                                        ggCAGCAGCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGgcagca          <  2:2658119/98‑1 (MQ=255)                                                                                          cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      >  2:2202559/1‑100 (MQ=255)                                                                                          cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      >  2:2077452/1‑100 (MQ=255)                                                                                          cagcagCGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCaaaaa      <  2:2058548/100‑1 (MQ=255)                                                                                             cagcGCACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActc   <  2:305882/100‑1 (MQ=255)                                                                                              agcGCACGGCAAATTCAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAActct  <  1:2620312/100‑1 (MQ=255)                                                                                                  |                                                                                                AGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTTTTCCTGCATCAACGGCAGCAGCGTACGGCAAATTGAGCGGATACAGTGAACGCGCATCATTAAGCCGTAATGACGACTCATAAACTCAAGATTATCAATGCGCAGGCAGCAAAAACTCT  >  NC_000913/2516592‑2516784

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: ‑