BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,409,046	C→T	A56V (GCC→GTC)	yjfK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,409,046	0	C	T	99.7%	1051.1 / ‑6.0	301	A56V (GCC→GTC)	yjfK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (146/154); total (146/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 332 total aligned reads displayed.

TCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGC  >  NC_000913/4408956‑4409137

tcAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAgtcg                                                                                           >  1:1755804/1‑93 (MQ=255)

  aaTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAgtcgcc                                                                                         >  2:3980960/1‑93 (MQ=255)

  aaTAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAgtcgcc                                                                                         >  2:6919083/1‑93 (MQ=255)

     aGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGTAATTTACGGTAGTCGCCGTc                                                                                      <  1:3765642/93‑1 (MQ=255)

     aGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTc                                                                                      <  1:1033003/93‑1 (MQ=255)

      gTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGATGCCGGGTGAGGAATTTACGGTAGTCGCCGTCa                                                                                     >  1:7615805/1‑93 (MQ=255)

        ggCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGc                                                                                   >  1:6166067/1‑93 (MQ=255)

           tttACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCcac                                                                                >  1:4963349/1‑93 (MQ=255)

           tttACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCcac                                                                                >  1:1251591/1‑93 (MQ=255)

            ttACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCcaca                                                                               <  1:4219557/93‑1 (MQ=255)

             tACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACAt                                                                              <  2:2009228/93‑1 (MQ=255)

             tACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACAt                                                                              >  1:5510090/1‑93 (MQ=255)

                gCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGa                                                                           <  1:2221479/93‑1 (MQ=255)

                 cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCgct                                                                          >  2:5876744/1‑91 (MQ=255)

                 cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGAt                                                                          <  1:902166/93‑1 (MQ=255)

                 cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGAt                                                                          <  2:5379516/93‑1 (MQ=255)

                 cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGAt                                                                          >  1:5701504/1‑93 (MQ=255)

                 cTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGAt                                                                          <  1:5085918/93‑1 (MQ=255)

                  tCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATc                                                                         <  1:7778128/93‑1 (MQ=255)

                    gACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTg                                                                       <  1:879606/93‑1 (MQ=255)

                    gACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTg                                                                       <  1:6224027/93‑1 (MQ=255)

                    gACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTg                                                                       <  2:300429/93‑1 (MQ=255)

                      cacGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTggg                                                                     >  1:11667/1‑93 (MQ=255)

                      cacGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTggg                                                                     >  2:7135198/1‑93 (MQ=255)

                       acGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggc                                                                    >  2:6148785/1‑93 (MQ=255)

                       acGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggc                                                                    <  1:7464085/93‑1 (MQ=255)

                       acGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggc                                                                    >  1:2724462/1‑93 (MQ=255)

                          ttAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCCGCCACATCGATCTGggcggc                                                                 >  1:4066724/1‑93 (MQ=255)

                          ttAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggcggc                                                                 >  1:6976058/1‑93 (MQ=255)

                          ttAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggcggc                                                                 <  1:5202209/93‑1 (MQ=255)

                          ttAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGggcggc                                                                 <  1:4669257/93‑1 (MQ=255)

                           tAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGgcggcg                                                                >  2:1653874/1‑93 (MQ=255)

                           tAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGgcggcg                                                                <  1:6302587/93‑1 (MQ=255)

                           tAGCGTTTCGTTTGCTGGAAGATCAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGgcggcg                                                                <  2:4282294/93‑1 (MQ=255)

                            aGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGACCTCGTcggcgg                                                               >  2:1151961/1‑93 (MQ=255)

                              cgtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTa                                                             <  2:367839/93‑1 (MQ=255)

                              cgtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTa                                                             <  1:7664880/93‑1 (MQ=255)

                               gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCTGTCAGCCACATCGATCTGGGCGGCGGTAg                                                            <  2:583005/93‑1 (MQ=255)

                               gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAg                                                            >  2:3499796/1‑93 (MQ=255)

                               gtttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAg                                                            <  1:6327588/93‑1 (MQ=255)

                                tttcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGt                                                           >  1:5135313/1‑93 (MQ=255)

                                  tcgtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCa                                                         <  1:3970094/93‑1 (MQ=255)

                                    gtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCCGCCACATCGATCTGGGCGGCGGTAGTCAGa                                                       >  2:4057692/1‑93 (MQ=255)

                                    gtttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGa                                                       >  1:2562425/1‑93 (MQ=255)

                                     tttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGAt                                                      >  1:6880120/1‑93 (MQ=255)

                                      ttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGAtt                                                     >  1:2304762/1‑93 (MQ=255)

                                      ttGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGCCAGAtt                                                     >  2:2590914/1‑93 (MQ=255)

                                      ttGCTGGAAGATGAATTGCTGATGGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGAtt                                                     <  2:4356081/93‑1 (MQ=255)

                                       tGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGAttt                                                    >  2:4997969/1‑93 (MQ=255)

                                          tGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTcc                                                 <  2:7756054/93‑1 (MQ=255)

                                          tGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTcc                                                 >  2:7529377/1‑93 (MQ=255)

                                            gAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGc                                               >  1:6749163/1‑93 (MQ=255)

                                            gAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGc                                               >  2:6632959/1‑93 (MQ=255)

                                              aGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGcta                                             <  2:2249091/93‑1 (MQ=255)

                                              aGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGcta                                             >  1:2891971/1‑93 (MQ=255)

                                               gATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGctac                                            <  2:4650114/93‑1 (MQ=255)

                                               gATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGctac                                            >  2:5082655/1‑93 (MQ=255)

                                                aTGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGctact                                           >  2:2100634/1‑93 (MQ=255)

                                                 tGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGctacta                                          >  2:3805934/1‑93 (MQ=255)

                                                  gAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCtactac                                         <  1:3268437/93‑1 (MQ=255)

                                                      tGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACtt                                     <  2:1420060/93‑1 (MQ=255)

                                                      tGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACtt                                     >  2:2526902/1‑93 (MQ=255)

                                                       gCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTc                                    <  2:6417276/93‑1 (MQ=255)

                                                       gCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTc                                    <  2:3946088/93‑1 (MQ=255)

                                                       gCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTc                                    >  1:205931/1‑93 (MQ=255)

                                                         ttaTTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAg                                  <  1:7859617/91‑1 (MQ=255)

                                                         tGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAg                                  >  2:2167485/1‑93 (MQ=255)

                                                         tGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAg                                  <  2:1482248/93‑1 (MQ=255)

                                                         tGATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAg                                  <  2:5250422/93‑1 (MQ=255)

                                                          gATTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAgg                                 <  1:7334266/93‑1 (MQ=255)

                                                           aTTGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGc                                >  2:7297012/1‑93 (MQ=255)

                                                            ttGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCg                               >  2:5502536/1‑93 (MQ=255)

                                                             tGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGa                              <  2:4074520/93‑1 (MQ=255)

                                                             tGCGCTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGa                              <  1:3544363/93‑1 (MQ=255)

                                                              gcgcTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGAt                             <  1:1895218/93‑1 (MQ=255)

                                                                gcTGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGa                           <  1:1067735/93‑1 (MQ=255)

                                                                  tGCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGt                         >  2:3408908/1‑93 (MQ=255)

                                                                   gCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGtt                        <  2:3226034/93‑1 (MQ=255)

                                                                   gCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGtt                        <  1:1493377/93‑1 (MQ=255)

                                                                   gCCGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGtt                        >  1:4527819/1‑93 (MQ=255)

                                                                    ccGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGttt                       <  1:243378/93‑1 (MQ=255)

                                                                    ccGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGttt                       <  2:682914/93‑1 (MQ=255)

                                                                     cGGGTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTc                      <  1:5261211/93‑1 (MQ=255)

                                                                      gggTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCt                     >  1:3181160/1‑93 (MQ=255)

                                                                      gggTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCt                     <  1:4840618/93‑1 (MQ=255)

                                                                       ggTGAGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTg                    <  2:2422899/93‑1 (MQ=255)

                                                                           agcaaTTTGCGGTAGGCCGCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCaaa                <  2:7709623/90‑1 (MQ=255)

                                                                           aGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCaaa                >  2:6885383/1‑93 (MQ=255)

                                                                           aGGAATTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCaaa                <  1:3452175/93‑1 (MQ=255)

                                                                               aTTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCaa            <  1:4171173/93‑1 (MQ=255)

                                                                               aTTTACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCaa            >  2:2300485/1‑93 (MQ=255)

                                                                                tttACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAAt           <  1:3835462/93‑1 (MQ=255)

                                                                                tttACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAAt           <  1:5994734/93‑1 (MQ=255)

                                                                                 ttACGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATa          <  2:879804/93‑1 (MQ=255)

                                                                                  tcgggTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATac         <  1:1631926/90‑1 (MQ=255)

                                                                                   aCGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacc        >  2:2959325/1‑93 (MQ=255)

                                                                                    cGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacca       >  2:1866531/1‑93 (MQ=255)

                                                                                    cGGTAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATacca       >  1:2143413/1‑93 (MQ=255)

                                                                                      gTAGTCCCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATaccacc     <  1:7152996/93‑1 (MQ=255)

                                                                                       tAGTCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACcg    <  2:4849396/93‑1 (MQ=255)

TCAATAGTGGCTTTACGCTCGACACGTTAGCGTTTCGTTTGCTGGAAGATGAATTGCTGATTGCGCTGCCGGGTGAGGAATTTACGGTAGCCGCCGTCAGCCACATCGATCTGGGCGGCGGTAGTCAGATTTTCCGCTACTACACTTCAGGCGATGAGTTTCTGCAAATCAATACCACCGGC  >  NC_000913/4408956‑4409137

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