BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,460,914	T→C	Y457C (TAC→TGC)	treC ←	trehalase 6‑P hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,460,914	0	T	C	99.7%	924.7 / ‑5.0	283	Y457C (TAC→TGC)	treC	trehalase 6‑P hydrolase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (155/127); minor base A (1/0); total (156/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 287 total aligned reads displayed.

ATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGT  >  NC_000913/4460825‑4461002

  aCAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAac                                                                                     >  2:6015058/1‑93 (MQ=255)

   cAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAaca                                                                                    >  2:4807839/1‑93 (MQ=255)

    aGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAacac                                                                                   >  2:4988619/1‑93 (MQ=255)

     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  >  1:2950533/1‑93 (MQ=255)

     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  2:3751052/93‑1 (MQ=255)

     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  1:2476915/93‑1 (MQ=255)

     gggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACAcc                                                                                  <  2:39721/93‑1 (MQ=255)

      ggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCg                                                                                 <  1:5373402/93‑1 (MQ=255)

      ggCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCg                                                                                 <  1:3193440/93‑1 (MQ=255)

       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                >  2:1515335/1‑93 (MQ=255)

       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                >  2:5529065/1‑93 (MQ=255)

       gCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGa                                                                                <  1:2714137/93‑1 (MQ=255)

        cTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAg                                                                               <  1:4345196/93‑1 (MQ=255)

         tGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAgg                                                                              <  2:3676240/93‑1 (MQ=255)

          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             >  1:1834515/1‑93 (MQ=255)

          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             <  1:4840938/93‑1 (MQ=255)

          gTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGa                                                                             >  1:4032184/1‑93 (MQ=255)

           tttGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGaa                                                                            >  1:946013/1‑93 (MQ=255)

           tttGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCGAGGaa                                                                            >  1:2098725/1‑93 (MQ=255)

                 agcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCACAACACCGAGGAATCATCg                                                                      >  2:3702465/1‑93 (MQ=255)

                 agcagATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCg                                                                      >  2:4435334/1‑93 (MQ=255)

                     gATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCa                                                                  >  2:5291917/1‑93 (MQ=255)

                     gATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCa                                                                  >  2:1325528/1‑93 (MQ=255)

                       tCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgc                                                                <  1:1397236/93‑1 (MQ=255)

                         cTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgcgc                                                              >  2:4754988/1‑93 (MQ=255)

                         cTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAgcgc                                                              >  1:1543900/1‑93 (MQ=255)

                          tGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGcgcg                                                             >  2:3362482/1‑93 (MQ=255)

                           ggTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCgg                                                            >  2:1021498/1‑93 (MQ=255)

                           ggTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCgg                                                            >  2:6172717/1‑93 (MQ=255)

                            gTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGc                                                           <  1:6087733/93‑1 (MQ=255)

                            gTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGc                                                           <  1:4194619/93‑1 (MQ=255)

                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         <  1:1152032/93‑1 (MQ=255)

                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         <  1:425405/93‑1 (MQ=255)

                              aaTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCtt                                                         >  1:4162005/1‑93 (MQ=255)

                               aTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTc                                                        >  2:2394044/1‑93 (MQ=255)

                               aTTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTc                                                        <  1:2852191/93‑1 (MQ=255)

                                 tGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACTGAGGAATCATCGGCCAGCGCGGCTTCTa                                                      <  1:3527911/93‑1 (MQ=255)

                                    cccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGt                                                   >  1:221800/1‑93 (MQ=255)

                                     ccATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGtt                                                  >  2:3768301/1‑93 (MQ=255)

                                      cATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTg                                                 >  1:4243533/1‑93 (MQ=255)

                                      cATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTg                                                 >  2:825312/1‑93 (MQ=255)

                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                <  1:1597787/93‑1 (MQ=255)

                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  2:5207747/1‑93 (MQ=255)

                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  1:1999197/1‑93 (MQ=255)

                                       aTGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGa                                                >  2:231689/1‑93 (MQ=255)

                                         gTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAtt                                              <  1:5462728/93‑1 (MQ=255)

                                         gTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAtt                                              >  2:3408603/1‑93 (MQ=255)

                                          tCAGGATGGCTTCCTTCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAttt                                             <  1:1757924/93‑1 (MQ=255)

                                          tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGAttt                                             <  2:5903307/93‑1 (MQ=255)

                                          tCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGATGAttt                                             >  1:4515427/1‑93 (MQ=255)

                                           cAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATttg                                            <  1:4898775/93‑1 (MQ=255)

                                             ggATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATttgtt                                          <  2:1282792/93‑1 (MQ=255)

                                             ggATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCAGCCAGCGCGGCTTCTACGTTGATttgtt                                          >  1:101018/1‑93 (MQ=255)

                                                 ggCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATa                                      <  2:1864724/93‑1 (MQ=255)

                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     <  1:5250961/93‑1 (MQ=255)

                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     >  1:578219/1‑93 (MQ=255)

                                                  gCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAg                                     >  1:636248/1‑93 (MQ=255)

                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  1:3636535/93‑1 (MQ=255)

                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  1:3784976/93‑1 (MQ=255)

                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGt                                    <  2:3037944/93‑1 (MQ=255)

                                                   cTTCCTGCTTACGCAGTGCGATTAACTTTTGGGAGGTGCAAAACACCCAGGAATCATCGGCCAGCGCGGCTTCTACGTTGCTTTGTTGATAGt                                    >  1:477787/1‑93 (MQ=255)

                                                    ttCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGtt                                   >  2:5870098/1‑93 (MQ=255)

                                                    tgccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGtt                                   <  2:3909628/91‑1 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  2:2726063/1‑93 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  1:2184692/1‑93 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  <  1:5718560/93‑1 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  2:3853017/1‑93 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  <  1:3278366/93‑1 (MQ=255)

                                                     tCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTa                                  >  1:2902277/1‑93 (MQ=255)

                                                      ccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAt                                 <  2:4070595/93‑1 (MQ=255)

                                                      ccTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTAt                                 >  2:3787310/1‑93 (MQ=255)

                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTTGTGATAGTTATc                                >  1:5623926/1‑93 (MQ=255)

                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:1477661/1‑93 (MQ=255)

                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:5200071/1‑93 (MQ=255)

                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  1:5118564/1‑93 (MQ=255)

                                                       cTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATc                                >  2:1759677/1‑93 (MQ=255)

                                                         gCTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGc                              <  2:4681820/93‑1 (MQ=255)

                                                          cTTACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGcc                             <  2:5673944/93‑1 (MQ=255)

                                                            tACGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCa                           >  2:3600752/1‑93 (MQ=255)

                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:2925888/93‑1 (MQ=255)

                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:748201/93‑1 (MQ=255)

                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          <  2:455048/93‑1 (MQ=255)

                                                             aCGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAg                          >  2:2821751/1‑93 (MQ=255)

                                                              cGCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAgg                         >  2:3252574/1‑93 (MQ=255)

                                                               gCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  1:2245177/1‑93 (MQ=255)

                                                               gCAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  2:6263732/1‑93 (MQ=255)

                                                               gCAATGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGc                        >  1:772177/1‑93 (MQ=255)

                                                                cAGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGcc                       >  1:622074/1‑93 (MQ=255)

                                                                 aGTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCa                      >  2:4648801/1‑93 (MQ=255)

                                                                  gTGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCaa                     >  1:3943662/1‑93 (MQ=255)

                                                                   tGCGATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAAt                    >  1:3993488/1‑93 (MQ=255)

                                                                      gATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCa                 <  1:1185556/93‑1 (MQ=255)

                                                                      gATTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAAGCCAATCCa                 >  2:3553560/1‑93 (MQ=255)

                                                                       aTTAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCAc                >  2:2582904/1‑93 (MQ=255)

                                                                        ttAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACg               >  2:156295/1‑93 (MQ=255)

                                                                        ttAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACg               >  2:24391/1‑93 (MQ=255)

                                                                         tAACTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACgg              >  2:2869978/1‑93 (MQ=255)

                                                                          aaCTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGt             <  1:3078043/93‑1 (MQ=255)

                                                                           aCTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGtt            <  2:3192721/93‑1 (MQ=255)

                                                                            cTTTTGGTAGGTGCAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTc           <  2:1853261/93‑1 (MQ=255)

ATACAGGGCTGTTTGGCAGCAGATCCTGGTAATTGCCCCATGTCAGGATGGCTTCCTGCTTACGCAGTGCGATTAACTTTTGGTAGGTGTAAAACACCGAGGAATCATCGGCCAGCGCGGCTTCTACGTTGATTTGTTGATAGTTATCGCCCAGGCCAATCCACGGTTCGCCAGCCGT  >  NC_000913/4460825‑4461002

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