BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,115,180	T→C	Q416R (CAG→CGG)	wcaK ←	putative galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,115,180	0	T	C	99.6%	819.3 / NA	251	Q416R (CAG→CGG)	wcaK	putative galactokinase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (157/93); minor base A (1/0); total (158/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

GTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCA  >  NC_000913/2115090‑2115269

gTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGc                                                                                         >  2:2994388/1‑93 (MQ=255)

 tcaaGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCa                                                                                        >  1:6618385/3‑93 (MQ=255)

  tAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCAc                                                                                       <  2:3770348/93‑1 (MQ=255)

   aaGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCAcc                                                                                      <  1:5888679/93‑1 (MQ=255)

   aaGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCAcc                                                                                      >  1:289829/1‑93 (MQ=255)

   aaGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCAcc                                                                                      <  1:2284621/93‑1 (MQ=255)

     ggACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCAt                                                                                    <  1:3424048/93‑1 (MQ=255)

     ggACGAAGGTTTCTGACGACAGCGGGAATTTCAGAAAAAAGAACCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACACACCGCACCAt                                                                                    >  1:2942599/1‑93 (MQ=255)

       aCGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCt                                                                                  <  1:5658934/93‑1 (MQ=255)

       aCGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCt                                                                                  >  1:1462022/1‑93 (MQ=255)

       aCGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCt                                                                                  >  1:2512251/1‑93 (MQ=255)

          aaGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCa                                                                               >  1:3334319/1‑93 (MQ=255)

           aGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCAt                                                                              >  2:1124618/1‑93 (MQ=255)

           aGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCAt                                                                              >  1:138201/1‑93 (MQ=255)

           aGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCAt                                                                              >  1:4287987/1‑93 (MQ=255)

            ggTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCAtt                                                                             >  2:770883/1‑93 (MQ=255)

              tttCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTcc                                                                           >  1:6080839/1‑93 (MQ=255)

                tCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTg                                                                         <  1:2171829/93‑1 (MQ=255)

                    acgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTg                                                                     <  2:5745397/93‑1 (MQ=255)

                     cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGa                                                                    >  1:2589381/1‑93 (MQ=255)

                     cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGa                                                                    <  2:171008/93‑1 (MQ=255)

                     cgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGa                                                                    >  2:4883888/1‑93 (MQ=255)

                      gacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGAc                                                                   <  1:4073675/93‑1 (MQ=255)

                        cAGCGGTAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGc                                                                 <  1:5465098/93‑1 (MQ=255)

                         aGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCt                                                                <  1:3045132/93‑1 (MQ=255)

                         aGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCt                                                                <  2:662592/93‑1 (MQ=255)

                          gCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTc                                                               >  1:6236032/1‑93 (MQ=255)

                          gCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTc                                                               <  1:1069685/93‑1 (MQ=255)

                           cGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCa                                                              >  2:5148229/1‑93 (MQ=255)

                            gggAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCAc                                                             >  1:1830027/1‑93 (MQ=255)

                            gggAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCAc                                                             >  1:3370490/1‑93 (MQ=255)

                             ggAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACg                                                            >  2:3908926/1‑93 (MQ=255)

                               aaTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGAc                                                          <  1:2751166/93‑1 (MQ=255)

                                aTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACt                                                         >  2:551385/1‑93 (MQ=255)

                                aTTTCAGTAAAAAGAAGACGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACt                                                         >  2:4426547/1‑93 (MQ=255)

                                 tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTg                                                        >  1:6702368/1‑93 (MQ=255)

                                 tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTg                                                        <  2:2388548/93‑1 (MQ=255)

                                 tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTg                                                        >  1:2456435/1‑93 (MQ=255)

                                   tCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                                                      >  1:3749545/1‑93 (MQ=255)

                                      gTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGc                                                   >  1:4963814/1‑93 (MQ=255)

                                      gTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGc                                                   <  1:2078235/93‑1 (MQ=255)

                                       tAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCt                                                  <  2:5575021/93‑1 (MQ=255)

                                        aaaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCtt                                                 <  1:4710594/93‑1 (MQ=255)

                                         aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTc                                                >  2:8410/1‑93 (MQ=255)

                                           aagaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc                                              >  1:5073254/1‑93 (MQ=255)

                                            agaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt                                             >  2:6417369/1‑93 (MQ=255)

                                             gaagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTa                                            >  2:492805/1‑93 (MQ=255)

                                              aagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           >  1:5311340/1‑93 (MQ=255)

                                              aagCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTaa                                           <  2:1528070/93‑1 (MQ=255)

                                               agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAg                                          <  1:361455/93‑1 (MQ=255)

                                               agCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAg                                          >  2:3304933/1‑93 (MQ=255)

                                               agCCGACCTTCATTTCACCTCCCCGATGCGATCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGCCGGCTTCACTAAg                                          >  2:4230264/1‑93 (MQ=255)

                                                gCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGt                                         >  2:3601586/1‑93 (MQ=255)

                                                gCCGACCGTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGt                                         >  1:5402255/1‑93 (MQ=255)

                                                  cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCCTTCCTGTCTGACGCTCACGACCGACGGCGTCACCAAGttg                                       >  1:1160577/1‑91 (MQ=255)

                                                  cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcg                                       >  1:3848063/1‑93 (MQ=255)

                                                   gACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgc                                      >  2:2263278/1‑93 (MQ=255)

                                                   gACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgc                                      >  2:1963776/1‑93 (MQ=255)

                                                    aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg                                     <  1:5808785/93‑1 (MQ=255)

                                                    aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg                                     >  1:1605167/1‑93 (MQ=255)

                                                    aCCTTCATATCACCTCCCCGGTGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg                                     >  1:6570079/1‑93 (MQ=255)

                                                     ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCgcgc                                    >  1:976224/1‑93 (MQ=255)

                                                     ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCgcgc                                    >  1:5517099/1‑93 (MQ=255)

                                                     ccTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCgcgc                                    >  2:5068277/1‑93 (MQ=255)

                                                      cTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                   >  2:3135636/1‑93 (MQ=255)

                                                       ttCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCAt                                  >  1:6367223/1‑93 (MQ=255)

                                                       ttCATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCAt                                  >  1:513878/1‑93 (MQ=255)

                                                        tCATTTCCCCTCCCCGATGCGCTCAAGCCCAGGCCGCAACGTCTGCATTTCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCAtt                                 <  2:3517304/93‑1 (MQ=255)

                                                         cATTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                                >  1:2822389/1‑93 (MQ=255)

                                                          aTTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCa                               <  2:2593578/93‑1 (MQ=255)

                                                          aTTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCa                               <  1:5216127/93‑1 (MQ=255)

                                                          aTTTCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCa                               >  1:2634888/1‑93 (MQ=255)

                                                            ttCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgc                             <  2:2462431/93‑1 (MQ=255)

                                                            ttCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgc                             >  2:2183662/1‑93 (MQ=255)

                                                            ttCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgc                             <  1:6369792/93‑1 (MQ=255)

                                                             tCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcg                            >  2:2081807/1‑93 (MQ=255)

                                                             tCACCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcg                            >  2:1761114/1‑93 (MQ=255)

                                                              cACCCCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcgc                           >  1:697665/1‑93 (MQ=255)

                                                               cccTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGcc                          <  1:4601085/92‑1 (MQ=255)

                                                               aCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGcc                          >  2:6065565/1‑93 (MQ=255)

                                                               aCCTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGcc                          >  2:5961432/1‑93 (MQ=255)

                                                                 cTCCCCGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCgg                        >  2:2768660/1‑93 (MQ=255)

                                                                  tcaccGATGCGTTCAAGCACAAACCGCACCATCTGCATTCCTGTCTTACGCTCACGACTGGCGGGTTCACTAAGTCGCGGATTCAGGGCCGGa                       >  1:5852072/4‑93 (MQ=255)

                                                                   ccccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGaa                      >  2:2526159/1‑93 (MQ=255)

                                                                     ccGATGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGc                    >  1:4454585/1‑93 (MQ=255)

                                                                        aTGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTgg                 >  1:632041/1‑93 (MQ=255)

                                                                         tGCGCTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGc                >  1:3089227/1‑93 (MQ=255)

                                                                           cgcTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGGCt              >  1:3307458/1‑93 (MQ=255)

                                                                           cgcTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTCGCCt              >  2:1196689/1‑93 (MQ=255)

                                                                             cTCAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTaa            <  2:4292881/93‑1 (MQ=255)

                                                                               cAAGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAgg          >  1:889515/1‑93 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         >  2:4729705/1‑93 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         >  2:4592439/1‑93 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         <  1:4740819/93‑1 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         >  2:6258671/1‑93 (MQ=255)

                                                                                aaGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt         >  2:1915224/1‑93 (MQ=255)

                                                                                 aGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        >  2:4009036/1‑93 (MQ=255)

                                                                                 aGCACAGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        >  1:17457/1‑93 (MQ=255)

                                                                                 aGCACAGACCGCACCACCCGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa        <  1:1796050/93‑1 (MQ=255)

                                                                                      aGACCGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGc   <  1:6109538/93‑1 (MQ=255)

GTTAAGGACGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCCGCA  >  NC_000913/2115090‑2115269

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