BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,886,722	G→A	G122S (GGC→AGC)	b3708 →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,886,722	0	G	A	99.6%	831.2 / ‑5.4	235	G122S (GGC→AGC)	b3708	tryptophanase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (123/110); total (124/111)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 244 total aligned reads displayed.

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811

ggTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGc                                                                                         >  2:4582060/1‑93 (MQ=255)

 gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGcc                                                                                        <  2:1116985/93‑1 (MQ=255)

 gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGcc                                                                                        <  2:4400798/93‑1 (MQ=255)

  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCACATCTATATTACGGTACTGCTTAAAAAACGCGAGCAGGAACAAAGCCt                                                                                       >  1:574131/1‑93 (MQ=255)

  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                       <  1:3977172/93‑1 (MQ=255)

  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                       <  2:1682564/93‑1 (MQ=255)

   tATCAATTCACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      <  1:2097348/93‑1 (MQ=255)

   tATCAATACACCATTCCGACTCACCAGGGCCGTGGAGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAATCCTg                                                                                      >  1:2533079/1‑93 (MQ=255)

    aTCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTgg                                                                                     >  1:4312260/1‑93 (MQ=255)

      cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAt                                                                                   >  2:1234241/1‑93 (MQ=255)

       aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATc                                                                                  <  2:3698142/93‑1 (MQ=255)

       aaTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATc                                                                                  >  1:2196419/1‑93 (MQ=255)

        aTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCg                                                                                 >  2:852594/1‑93 (MQ=255)

         tACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgc                                                                                >  1:2184175/1‑93 (MQ=255)

          tcacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgca                                                                               <  1:977542/92‑1 (MQ=255)

          acacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgca                                                                               >  2:582579/1‑93 (MQ=255)

          acacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgca                                                                               >  2:3656993/1‑93 (MQ=255)

          acacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgca                                                                               >  2:2386012/1‑93 (MQ=255)

           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                              <  1:1207833/93‑1 (MQ=255)

           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                              >  2:1937203/1‑93 (MQ=255)

           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                              >  1:3055255/1‑93 (MQ=255)

            cccATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                             <  2:2675642/92‑1 (MQ=255)

             ccATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagca                                                                            >  2:2505350/1‑93 (MQ=255)

                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         >  1:1011744/1‑93 (MQ=255)

                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         <  1:2901111/93‑1 (MQ=255)

                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         <  1:1430457/93‑1 (MQ=255)

                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         <  1:3682022/93‑1 (MQ=255)

                   cGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtgg                                                                      <  1:4475377/93‑1 (MQ=255)

                   cGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtgg                                                                      >  2:2783017/1‑93 (MQ=255)

                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAAtggt                                                                     <  2:744151/93‑1 (MQ=255)

                     gcTCAGCAGGGCCGTGGCGCATAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    <  1:2595007/92‑1 (MQ=255)

                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    >  2:1292617/1‑93 (MQ=255)

                      cTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtgg                                                                   <  2:3327190/93‑1 (MQ=255)

                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                  >  2:3298235/1‑93 (MQ=255)

                         aCCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGGACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGt                                                                >  1:2686969/1‑93 (MQ=255)

                            aGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtct                                                             >  2:4504007/1‑93 (MQ=255)

                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                           <  1:3714845/93‑1 (MQ=255)

                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                           <  2:3432264/93‑1 (MQ=255)

                                ccGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTaa                                                         <  2:603619/93‑1 (MQ=255)

                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                        <  1:2869541/93‑1 (MQ=255)

                                   tGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTa                                                      >  1:1159193/1‑93 (MQ=255)

                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                    <  2:2444739/93‑1 (MQ=255)

                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                    <  1:3582199/93‑1 (MQ=255)

                                        cAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTcc                                                 <  2:2880171/93‑2 (MQ=255)

                                        cAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                                                 <  1:3433945/93‑1 (MQ=255)

                                        cAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                                                 <  1:613392/93‑1 (MQ=255)

                                         agagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                                                <  1:1406362/93‑1 (MQ=255)

                                         agagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                                                >  1:54830/1‑93 (MQ=255)

                                         agagCAAATCTATATACCGGTACTGATTAAAAAAAGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAACTGCTGGCGTTCTCAAACTCTCTCtt                                                >  2:3680057/1‑93 (MQ=255)

                                          gagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                                               <  1:4414749/93‑1 (MQ=255)

                                           agCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTg                                              >  1:3348747/1‑93 (MQ=255)

                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             >  2:4087486/1‑93 (MQ=255)

                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             >  1:3318369/1‑93 (MQ=255)

                                             cAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGAt                                            >  2:839115/1‑93 (MQ=255)

                                               aaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                                          <  2:2350606/93‑1 (MQ=255)

                                               aaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                                          >  2:4655227/1‑93 (MQ=255)

                                                aTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                                         <  1:1562801/93‑1 (MQ=255)

                                                aTCTATATTCCGGTACTGATTAAAAAAAGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                                         >  2:4023441/1‑93 (MQ=255)

                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        <  1:2272848/93‑1 (MQ=255)

                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCCCg                                      >  1:3161099/1‑93 (MQ=255)

                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                                      >  1:4714623/1‑93 (MQ=255)

                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                                      >  1:415001/1‑93 (MQ=255)

                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                                      >  2:2064711/1‑93 (MQ=255)

                                                   tataTTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACg                                      <  2:1474296/93‑1 (MQ=255)

                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    >  2:2157646/1‑93 (MQ=255)

                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  1:1936212/93‑1 (MQ=255)

                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  1:4292773/93‑1 (MQ=255)

                                                      atTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAg                                   >  1:174001/1‑93 (MQ=255)

                                                      atTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAg                                   >  1:82762/1‑93 (MQ=255)

                                                        tCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg                                 >  2:3864689/1‑93 (MQ=255)

                                                        tCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg                                 <  2:2760020/93‑1 (MQ=255)

                                                             tACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa                            <  1:708869/93‑1 (MQ=255)

                                                              aCTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAg                           >  1:1632677/1‑93 (MQ=255)

                                                                tGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc                         >  1:1677868/1‑93 (MQ=255)

                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                        >  1:1997571/1‑93 (MQ=255)

                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                        <  1:3972182/93‑1 (MQ=255)

                                                                  aTTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAg                       >  1:3115708/1‑93 (MQ=255)

                                                                  aTTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAg                       >  1:2514189/1‑93 (MQ=255)

                                                                     aaaaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc                    <  1:4332355/93‑1 (MQ=255)

                                                                     aaaaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc                    >  1:408361/1‑93 (MQ=255)

                                                                       aaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCaa                  <  1:2730730/93‑1 (MQ=255)

                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                 <  2:2216968/93‑1 (MQ=255)

                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                 >  2:1874267/1‑93 (MQ=255)

                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                 <  2:499723/93‑1 (MQ=255)

                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGACCAAc                 >  1:1329946/1‑93 (MQ=255)

                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg                >  1:4557131/1‑93 (MQ=255)

                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg                >  1:38902/1‑93 (MQ=255)

                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATAAACg                >  1:13484/1‑93 (MQ=255)

                                                                          aCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACgg               >  1:2420219/1‑93 (MQ=255)

                                                                           cgcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCTGTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGc              <  1:46373/93‑1 (MQ=255)

                                                                           cgcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGc              <  1:4705590/93‑1 (MQ=255)

                                                                              gAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGCTACCACGCAGGGCCATAGCCAGATCAACGGCTGt           <  2:2978240/93‑1 (MQ=255)

                                                                                  aGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCg       <  2:702612/93‑1 (MQ=255)

                                                                                  aGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCg       <  2:3110253/93‑1 (MQ=255)

                                                                                   ggAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGt      >  1:2626967/1‑93 (MQ=255)

                                                                                     aaaaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGc    >  1:2374951/1‑93 (MQ=255)

                                                                                      aaaaAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCg   >  1:786807/1‑93 (MQ=255)

                                                                                      aaaaAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCg   <  2:1606533/93‑1 (MQ=255)

                                                                                      aaaaAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCg   >  2:2731391/1‑93 (MQ=255)

                                                                                       aaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGt  >  2:1797582/1‑93 (MQ=255)

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811

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