BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,205,190	T→C	G492G (GGT→GGC)	yehM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,205,190	0	T	C	96.2%	648.8 / NA	213	G492G (GGT→GGC)	yehM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (89/116); minor base A (8/0); total (97/116)
Fisher's exact test for biased strand distribution p-value = 1.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

AGCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGTCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCTGG  >  NC_000913/2205100‑2205278

aGCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGGCt                                                                                        >  2:2361748/1‑93 (MQ=255)

 gCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTg                                                                                       <  2:5623304/93‑1 (MQ=255)

 gCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTg                                                                                       >  2:5508705/1‑93 (MQ=255)

 gCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTg                                                                                       >  2:3482371/1‑93 (MQ=255)

  ccGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTgg                                                                                      <  1:5841290/93‑1 (MQ=255)

  ccGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTgg                                                                                      >  2:5813422/1‑93 (MQ=255)

  ccGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTgg                                                                                      <  1:2266552/93‑1 (MQ=255)

    gcacgcaATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCg                                                                                    >  1:6202653/1‑93 (MQ=255)

      acgcaATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGTCGTTAGCGGGCCTGGCGTc                                                                                  <  1:6238711/93‑1 (MQ=255)

      acgcaATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTc                                                                                  >  1:2457433/1‑93 (MQ=255)

        gcaATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCt                                                                                <  1:2070586/93‑1 (MQ=255)

          aaTAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTc                                                                              >  2:3376120/1‑93 (MQ=255)

          aaTAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTc                                                                              <  1:1027339/93‑1 (MQ=255)

             aaaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGt                                                                           <  2:358400/93‑1 (MQ=255)

             aaaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGt                                                                           >  2:2272387/1‑93 (MQ=255)

              aaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTc                                                                          <  2:3303489/93‑1 (MQ=255)

              aaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTc                                                                          >  1:5044069/1‑93 (MQ=255)

              aaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTc                                                                          <  2:3779480/93‑1 (MQ=255)

              aaTTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTc                                                                          >  1:1701459/1‑93 (MQ=255)

                  aGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCTAGGCGGCGTTAGCGGGCCTGGCGTGCTTCAGTcaaca                                                                      >  2:2640530/1‑93 (MQ=255)

                  aGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTcaaca                                                                      >  1:344510/1‑93 (MQ=255)

                    aaGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAAt                                                                    <  2:3075029/93‑1 (MQ=255)

                    aaGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAAt                                                                    <  1:1875038/93‑1 (MQ=255)

                      gCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTCAACAATTa                                                                  >  1:1691710/1‑93 (MQ=255)

                       ccGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTAc                                                                 <  2:4011534/93‑1 (MQ=255)

                       ccGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTAc                                                                 <  1:3893625/93‑1 (MQ=255)

                        cGATATGCTGGACGCAGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACt                                                                <  1:2574246/93‑1 (MQ=255)

                          atatGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTgg                                                              >  1:3111753/1‑93 (MQ=255)

                          atatGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTgg                                                              <  1:3521880/93‑1 (MQ=255)

                           tatGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGa                                                             <  2:1132451/93‑1 (MQ=255)

                           tatGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGa                                                             <  2:1948043/93‑1 (MQ=255)

                              tctagACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCa                                                          <  1:4905551/89‑1 (MQ=255)

                              gCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCa                                                          >  2:4783231/1‑93 (MQ=255)

                               gtGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCaa                                                         <  2:4818973/92‑1 (MQ=255)

                                tGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAAc                                                        >  1:432692/1‑93 (MQ=255)

                                tGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAAc                                                        <  1:1384924/93‑1 (MQ=255)

                                tGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTCAACAATTACTGGAGCAAc                                                        >  1:1821897/1‑93 (MQ=255)

                                 ggACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACt                                                       >  2:2864788/1‑93 (MQ=255)

                                 ggACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTCAACAATTACTGGAGCAACt                                                       >  1:1624972/1‑93 (MQ=255)

                                  gACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTc                                                      >  2:266719/1‑93 (MQ=255)

                                    cgcgGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTcac                                                    <  2:2716919/93‑1 (MQ=255)

                                     gcgggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTcaca                                                   >  2:1061529/1‑93 (MQ=255)

                                      cgggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCacac                                                  <  2:772894/93‑1 (MQ=255)

                                      cgggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCacac                                                  <  2:4199238/93‑1 (MQ=255)

                                      cgggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCacac                                                  <  2:173872/93‑1 (MQ=255)

                                       gggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACt                                                 <  1:4735804/93‑1 (MQ=255)

                                        ggcggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTa                                                >  2:3308695/1‑93 (MQ=255)

                                           ggAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTa                                             >  1:1280259/1‑93 (MQ=255)

                                            gAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTaa                                            <  2:3279485/93‑1 (MQ=255)

                                             aaTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAAt                                           >  2:524556/1‑93 (MQ=255)

                                              aTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATc                                          <  2:5668363/93‑1 (MQ=255)

                                              aTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATc                                          <  1:5483067/93‑1 (MQ=255)

                                              aTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATc                                          >  1:6466647/1‑93 (MQ=255)

                                              aTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGAAACTCACACTATTAATc                                          >  1:4036065/1‑93 (MQ=255)

                                              aTCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGCAGCAACTCACACTATTAATc                                          <  1:2679540/93‑1 (MQ=255)

                                               tCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCg                                         >  1:1017885/1‑93 (MQ=255)

                                               tCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCg                                         >  1:3327757/1‑93 (MQ=255)

                                                cGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGc                                        <  1:6036910/93‑1 (MQ=255)

                                                cGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGc                                        >  2:1356263/1‑93 (MQ=255)

                                                 ggCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGcc                                       <  1:2753356/93‑1 (MQ=255)

                                                   cAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCa                                     <  1:540863/93‑1 (MQ=255)

                                                    aGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAg                                    <  2:6375685/93‑1 (MQ=255)

                                                     gTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAgg                                   <  2:2576382/93‑1 (MQ=255)

                                                       aTCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGaa                                 <  2:2298476/93‑1 (MQ=255)

                                                       aTCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGaa                                 >  2:2611519/1‑93 (MQ=255)

                                                        tCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGACTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGaaa                                >  2:2473148/1‑93 (MQ=255)

                                                          accacATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAAt                              <  2:1406757/93‑1 (MQ=255)

                                                           ccacATGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAAACGCCCAGGAAAATc                             >  2:395323/1‑93 (MQ=255)

                                                            cacaTGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCa                            >  1:2984059/1‑93 (MQ=255)

                                                             acaTGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCaa                           <  1:5562151/93‑1 (MQ=255)

                                                               aTGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAAtt                         >  2:655412/1‑93 (MQ=255)

                                                               aTGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAAtt                         >  2:4904371/1‑93 (MQ=255)

                                                                tGTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAAttt                        <  1:1474795/93‑1 (MQ=255)

                                                                 gTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTg                       <  2:3264019/93‑1 (MQ=255)

                                                                 gTCTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTg                       >  2:1011101/1‑93 (MQ=255)

                                                                   cTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGcc                     <  2:6041576/93‑1 (MQ=255)

                                                                   cTTAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGcc                     >  1:2210030/1‑93 (MQ=255)

                                                                   cTTAGCCAGGCGGCGTTAGCGGGCCAGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGcc                     >  1:3627072/1‑93 (MQ=255)

                                                                     tAGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGa                   >  2:6188278/1‑93 (MQ=255)

                                                                      aGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGaa                  >  1:6273944/1‑93 (MQ=255)

                                                                      aGCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGaa                  >  1:1140159/1‑93 (MQ=255)

                                                                       gCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGaaa                 <  1:2828823/93‑1 (MQ=255)

                                                                       gCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGaaa                 >  2:3709936/1‑93 (MQ=255)

                                                                       gCCAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGaaa                 >  1:264573/1‑93 (MQ=255)

                                                                        ccAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAAt                <  2:4622851/93‑1 (MQ=255)

                                                                         cAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATg               <  1:1351741/93‑1 (MQ=255)

                                                                         cAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATg               <  1:4525790/93‑1 (MQ=255)

                                                                         cAGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATg               <  1:2772382/93‑1 (MQ=255)

                                                                          aGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATgg              >  1:3729712/1‑93 (MQ=255)

                                                                          aGGCGGCGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATgg              <  2:1895879/93‑1 (MQ=255)

                                                                           ggcggcGTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATggg             >  1:6178629/1‑93 (MQ=255)

                                                                               gcgTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAg         <  2:5971191/93‑1 (MQ=255)

                                                                                cgTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAgg        <  2:6352550/93‑1 (MQ=255)

                                                                                cgTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAgg        <  2:2062232/93‑1 (MQ=255)

                                                                                 gTTAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGgc       <  1:4121015/93‑1 (MQ=255)

                                                                                  ttAGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGgcg      >  2:1435647/1‑93 (MQ=255)

                                                                                    aGCGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCt    >  1:5896497/1‑93 (MQ=255)

                                                                                     gCGGGCCTGGCGTCCTTCAGTNAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCTg   <  2:3184783/93‑1 (MQ=255)

                                                                                      cGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCTgg  <  2:4512228/93‑1 (MQ=255)

                                                                                      cGGGCCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCTgg  <  1:379968/93‑1 (MQ=255)

AGCCGCACGCAATAAATTAGAAGCCGATATGCTGGACGCGGGCGGAATCGGCAGTATCACCACATGTCTTAGCCAGGCGGCGTTAGCGGGTCTGGCGTCCTTCAGTCAACAATTACTGGAGCAACTCACACTATTAATCGCCCAGGAAAATCAATTTGCCGAAATGGGCCAGGCGCTGG  >  NC_000913/2205100‑2205278

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