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.8%	1633.1 / ‑6.3	448	G122S (GGC→AGC)	b3708	tryptophanase
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (220/226); total (221/227)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 468 total aligned reads displayed.

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811

 gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAACAGGcc                                                                                        >  1:7195748/1‑93 (MQ=255)

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

  ttATCAATACACCATGCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                       <  2:5962573/93‑1 (MQ=255)

   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      >  2:6434852/1‑93 (MQ=255)

   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      <  1:7743252/93‑1 (MQ=255)

    aTCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAACGCCTgg                                                                                     >  2:5306893/1‑93 (MQ=255)

     tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTggc                                                                                    >  1:2430998/1‑92 (MQ=255)

     tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGa                                                                                    >  1:7464702/1‑93 (MQ=255)

      cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAt                                                                                   <  2:167223/93‑1 (MQ=255)

      cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAt                                                                                   <  1:3688888/93‑1 (MQ=255)

      cAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGAt                                                                                   <  2:7218024/93‑1 (MQ=255)

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

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

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

            acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                             >  2:1805700/1‑93 (MQ=255)

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

              cATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaa                                                                           <  1:6173560/93‑1 (MQ=255)

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

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

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

                  ccGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtg                                                                       <  1:963451/93‑1 (MQ=255)

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

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

                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                     <  1:3780532/93‑1 (MQ=255)

                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                     <  1:4367381/93‑1 (MQ=255)

                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggt                                                                     >  2:7539855/1‑93 (MQ=255)

                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCCCAGCAACAtggg                                                                     >  1:6747708/1‑92 (MQ=255)

                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGACAAAAGCCTGGATCGCAGCAACAtgggg                                                                    >  2:798183/1‑91 (MQ=255)

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

                      cTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtgg                                                                   >  2:4920120/1‑93 (MQ=255)

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

                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                  <  1:6385041/93‑1 (MQ=255)

                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                  <  2:6964253/93‑1 (MQ=255)

                        cACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCg                                                                 <  2:5520333/93‑1 (MQ=255)

                         aCCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGt                                                                <  2:7786306/93‑1 (MQ=255)

                          ccAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGtt                                                               <  2:4405405/93‑1 (MQ=255)

                           cAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAACATGGTGGCGTtc                                                              >  1:3970223/1‑93 (MQ=255)

                           cAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtc                                                              >  2:4520265/1‑93 (MQ=255)

                           cAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtc                                                              >  2:6727417/1‑93 (MQ=255)

                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                                            >  1:4240718/1‑93 (MQ=255)

                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                                            <  1:3110471/93‑1 (MQ=255)

                                ccGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTaa                                                         >  1:6070257/1‑93 (MQ=255)

                                  gTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACt                                                       >  1:4128471/1‑93 (MQ=255)

                                   tGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAGAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTa                                                      <  2:6307398/93‑1 (MQ=255)

                                    ggCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAt                                                     <  1:272280/93‑1 (MQ=255)

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

                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                    >  1:767351/1‑93 (MQ=255)

                                        cAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCt                                                 >  1:2649734/1‑93 (MQ=255)

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

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

                                           agCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTg                                              <  1:5290795/93‑1 (MQ=255)

                                           agCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAATTCTTTg                                              <  2:2766693/93‑1 (MQ=255)

                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             <  1:5701414/93‑1 (MQ=255)

                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             <  2:2374687/93‑1 (MQ=255)

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

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

                                              aaaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                                           >  2:506370/1‑93 (MQ=255)

                                               aaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGCTCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAc                                          >  1:2813143/1‑93 (MQ=255)

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

                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        >  2:306498/1‑93 (MQ=255)

                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        >  2:3785898/1‑93 (MQ=255)

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

                                                  cTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                                       <  1:7585170/93‑1 (MQ=255)

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

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

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

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

                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATAGAAGAAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  2:6000276/93‑1 (MQ=255)

                                                      atTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAg                                   <  2:7340693/93‑1 (MQ=255)

                                                      atTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAg                                   <  1:1304704/93‑1 (MQ=255)

                                                       ttCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg                                  <  2:936607/93‑1 (MQ=255)

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

                                                        tCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggc                                 >  1:1229493/1‑92 (MQ=255)

                                                         ccGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGc                                >  1:161828/1‑93 (MQ=255)

                                                           ggTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCa                              >  1:3299418/1‑93 (MQ=255)

                                                            gTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt                             >  2:3453861/1‑93 (MQ=255)

                                                             tACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACCATTTCTTTGATACCACGCAGGGCCATa                            >  1:5614259/1‑93 (MQ=255)

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

                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc                          <  2:6768210/93‑1 (MQ=255)

                                                                tGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGTTACCACGCAGGGCCATAgac                         >  2:3339107/1‑91 (MQ=255)

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

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

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

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

                                                                   ttAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGa                      >  1:6894421/1‑93 (MQ=255)

                                                                    tAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGAt                     <  2:1016151/93‑1 (MQ=255)

                                                                    tAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGAt                     <  1:5508590/93‑1 (MQ=255)

                                                                     aaaaaaCGTGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATc                    >  2:2930520/1‑93 (MQ=255)

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

                                                                            gcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCt             >  1:20566/1‑93 (MQ=255)

                                                                            gcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCt             <  2:6458640/93‑1 (MQ=255)

                                                                            gcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCt             <  2:5827889/93‑1 (MQ=255)

                                                                                gCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTAc         <  2:4813300/93‑1 (MQ=255)

                                                                                 cAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTAcc        <  1:6446366/93‑1 (MQ=255)

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

                                                                                     aaaaaaGCCTGGATGGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGc    <  1:4073130/93‑1 (MQ=255)

                                                                                     aaaaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGc    <  2:6569053/93‑1 (MQ=255)

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

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

                                                                                       aaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGt  <  2:2551734/93‑1 (MQ=255)

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811

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