Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AC_000091	3,222,265	T→C	L326P (CTG→CCG)	ebgA →	cryptic beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AC_000091	3,222,265	0	T	C	100.0%	384.6 / NA	98	L326P (CTG→CCG)	ebgA	cryptic beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (49/49);  total (49/49)

CATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAAT  >  AC_000091/3222179‑3222349                                                                                       |                                                                                     cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGtt                                                                                   >  2:2130902/1‑90 (MQ=255) cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGtt                                                                                   <  1:1527116/90‑1 (MQ=255) cATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGtt                                                                                   >  1:2004131/1‑90 (MQ=255)  aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTc                                                                                  >  2:833950/1‑90 (MQ=255)  aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTc                                                                                  <  2:1870115/90‑1 (MQ=255)  aTGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGACATCAAAGTGCGCGACGGTCCGTTc                                                                                  <  2:2477352/90‑1 (MQ=255)   tGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCt                                                                                 <  1:2530311/90‑1 (MQ=255)   tGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCt                                                                                 <  2:1830559/90‑1 (MQ=255)    gACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTg                                                                                >  1:975763/1‑90 (MQ=255)      cGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCCCAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGa                                                                              <  2:1326436/90‑1 (MQ=255)      cGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGa                                                                              >  1:2088106/1‑90 (MQ=255)      cGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGa                                                                              <  1:2716500/90‑1 (MQ=255)       gCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGAt                                                                             >  2:3484390/1‑90 (MQ=255)        cTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATc                                                                            <  1:794065/90‑1 (MQ=255)        cTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATc                                                                            <  1:3212179/90‑1 (MQ=255)        cTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATc                                                                            >  2:964547/1‑90 (MQ=255)         tGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCa                                                                           >  1:1083773/1‑90 (MQ=255)           aaaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAAt                                                                         <  2:2462184/90‑1 (MQ=255)           aaaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAAt                                                                         <  2:2512102/90‑1 (MQ=255)           aaaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAAt                                                                         <  1:1767332/90‑1 (MQ=255)            aaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATa                                                                        >  1:3462115/1‑90 (MQ=255)            aaGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATa                                                                        <  1:3181046/90‑1 (MQ=255)             aGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATaa                                                                       <  1:2082921/90‑1 (MQ=255)             aGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATaa                                                                       >  1:3247164/1‑90 (MQ=255)             aGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATaa                                                                       >  1:2376660/1‑90 (MQ=255)              gACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAAc                                                                      <  1:2917610/90‑1 (MQ=255)              gACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAAc                                                                      >  2:1354911/1‑90 (MQ=255)              gACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAAc                                                                      <  2:1354105/90‑1 (MQ=255)                cGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCg                                                                    <  1:3033238/90‑1 (MQ=255)                cGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCg                                                                    <  1:1521328/90‑1 (MQ=255)                    aaCGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTAt                                                                >  2:331801/1‑90 (MQ=255)                      cGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTAtgt                                                              >  2:3130619/1‑90 (MQ=255)                      cGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTAtgt                                                              <  2:1706769/90‑1 (MQ=255)                       ggCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTAtgtg                                                             >  2:1936288/1‑90 (MQ=255)                       ggCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTAtgtg                                                             >  2:693678/1‑90 (MQ=255)                        gCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAACGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGa                                                            <  1:730201/90‑1 (MQ=255)                        gCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGa                                                            <  1:1882537/90‑1 (MQ=255)                        gCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGa                                                            <  2:2741462/90‑1 (MQ=255)                         cAACGTTCTGGAAGTGGTGCCACAACTCGTTTGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAt                                                           <  2:2658030/90‑1 (MQ=255)                         cAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAt                                                           >  1:732043/1‑90 (MQ=255)                          aaCGTTCTGGAAGTGGTGCCCCAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtg                                                          <  2:2472990/90‑1 (MQ=255)                          aaCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtg                                                          <  2:2167793/90‑1 (MQ=255)                          aaCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtg                                                          <  2:2206703/90‑1 (MQ=255)                           aCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgc                                                         >  1:1748305/1‑90 (MQ=255)                           aCGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgc                                                         <  2:1350195/90‑1 (MQ=255)                            cGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgct                                                        <  2:2541819/90‑1 (MQ=255)                             gTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgctg                                                       <  2:2847546/90‑1 (MQ=255)                             gTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgctg                                                       <  2:130930/90‑1 (MQ=255)                              ttCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGAtgctgc                                                      <  1:963955/90‑1 (MQ=255)                               tCTGGAAGTGGTGCGACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCa                                                     >  2:1439037/1‑90 (MQ=255)                               tCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCa                                                     <  2:528277/90‑1 (MQ=255)                               tCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCa                                                     <  2:3043254/90‑1 (MQ=255)                               tCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCa                                                     <  1:2908799/90‑1 (MQ=255)                               tCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCCGCa                                                     <  2:163610/90‑1 (MQ=255)                                cTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCAc                                                    >  1:3400640/1‑90 (MQ=255)                                cTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCAACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCAc                                                    <  2:255555/90‑1 (MQ=255)                                 tGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACg                                                   <  2:1508951/90‑1 (MQ=255)                                  ggAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACgg                                                  >  2:2631422/1‑90 (MQ=255)                                   gAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGc                                                 >  2:793954/1‑90 (MQ=255)                                    aaGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCg                                                >  2:2219966/1‑90 (MQ=255)                                     aGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGt                                               >  2:3350076/1‑90 (MQ=255)                                        ggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCaa                                            >  1:1438529/1‑90 (MQ=255)                                        ggtgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCaa                                            <  2:1310241/90‑1 (MQ=255)                                          tgCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAAcc                                          <  2:1463934/90‑1 (MQ=255)                                           gCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCg                                         <  2:348192/90‑1 (MQ=255)                                            ccATAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGt                                        <  2:1495521/90‑1 (MQ=255)                                            ccACAACGCGTTTGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGCTATGTGATGCTGCACGTCGTCAACCGt                                        <  2:1006726/90‑1 (MQ=255)                                               caACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCAc                                     >  2:1600710/1‑90 (MQ=255)                                               caACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCAc                                     <  1:2097171/90‑1 (MQ=255)                                                aaCGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACg                                    >  2:3014291/1‑90 (MQ=255)                                                   gcgTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACa                                 >  1:2339927/1‑90 (MQ=255)                                                    cgTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACaa                                >  1:2101638/1‑90 (MQ=255)                                                     gTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGAGAAc                               >  2:3302902/1‑90 (MQ=255)                                                     gTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAAc                               <  2:1838335/90‑1 (MQ=255)                                                       tGGCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGa                             >  2:244976/1‑90 (MQ=255)                                                        ggCTTCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATCCTGCACGGCGTCAACCGACACGACAACGat                            <  2:192522/90‑1 (MQ=255)                                                           ttCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGatcat                         >  2:927158/1‑90 (MQ=255)                                                           ttCCGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGatcat                         >  1:2398602/1‑90 (MQ=255)                                                             ccGTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCg                       >  2:364687/1‑90 (MQ=255)                                                               gTGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCa                     >  1:519435/1‑90 (MQ=255)                                                                tGATATCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCaa                    >  1:578395/1‑90 (MQ=255)                                                                  atatCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAg                  >  2:1530998/1‑90 (MQ=255)                                                                     tCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGcc               >  2:1437084/1‑90 (MQ=255)                                                                     tCAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGcc               >  1:2710708/1‑90 (MQ=255)                                                                      cAAAGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcg              >  1:2464340/1‑90 (MQ=255)                                                                       aaaGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcgc             <  1:2782319/90‑1 (MQ=255)                                                                        aaGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcgcg            <  2:3300444/90‑1 (MQ=255)                                                                        aaGTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCcgcg            >  1:3043123/1‑90 (MQ=255)                                                                          gTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc          >  2:3347866/1‑90 (MQ=255)                                                                          gTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc          >  1:2443488/1‑90 (MQ=255)                                                                          gTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc          >  2:852293/1‑90 (MQ=255)                                                                          gTGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGcc          >  2:3037797/1‑90 (MQ=255)                                                                           tGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCg         >  2:1206591/1‑90 (MQ=255)                                                                           tGCGCGACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCg         >  1:373485/1‑90 (MQ=255)                                                                             cgcgACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGtt       >  2:278610/1‑90 (MQ=255)                                                                             cgcgACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGtt       >  1:2723960/1‑90 (MQ=255)                                                                              gcgACGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTg      >  2:2486518/1‑90 (MQ=255)                                                                                  cGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAAt  <  2:3321319/90‑1 (MQ=255)                                                                                  cGGTCCGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAAt  <  1:1441957/90‑1 (MQ=255)                                                                                       |                                                                                     CATGACGCTGAAAGACGCCAACGGCAACGTTCTGGAAGTGGTGCCACAACGCGTTGGCTTCCGTGATATCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAAT  >  AC_000091/3222179‑3222349

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