BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,563,909	A→G	W303R (TGG→CGG)	yjiN ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,563,909	0	A	G	99.7%	1169.5 / ‑5.2	366	W303R (TGG→CGG)	yjiN	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (181/184); total (181/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 374 total aligned reads displayed.

CAAACCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTT  >  NC_000913/4563819‑4563998

cAAACCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCcgca                                                                                         >  2:6502548/1‑93 (MQ=255)

cAAACCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCcgca                                                                                         >  1:1828170/1‑93 (MQ=255)

  aaCCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAAt                                                                                       <  1:3888885/93‑1 (MQ=255)

  aaCCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAAt                                                                                       <  1:3742283/93‑1 (MQ=255)

  aaCCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAAt                                                                                       >  2:1130260/1‑93 (MQ=255)

   aCCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAAtt                                                                                      <  1:1666906/93‑1 (MQ=255)

     cATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCa                                                                                    <  2:2065477/93‑1 (MQ=255)

     cATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCa                                                                                    >  2:3766583/1‑93 (MQ=255)

      aTTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCAc                                                                                   >  2:3402671/1‑93 (MQ=255)

      aTTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCAc                                                                                   <  1:3416345/93‑1 (MQ=255)

       ttGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACt                                                                                  >  1:198467/1‑93 (MQ=255)

       ttGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACt                                                                                  <  1:2546312/93‑1 (MQ=255)

         gtcccGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTga                                                                                <  1:1162666/91‑1 (MQ=255)

           cccGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTgaga                                                                              <  1:4371401/93‑1 (MQ=255)

            ccGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGAt                                                                             >  1:7354911/1‑93 (MQ=255)

            ccGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGAt                                                                             >  1:1080656/1‑93 (MQ=255)

              gCTCGTGCGATGCGTTCTTTCACACGAGCATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAg                                                                           >  1:1021719/1‑93 (MQ=255)

              gCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAg                                                                           >  1:2454688/1‑93 (MQ=255)

                tCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCg                                                                         >  1:1382984/1‑93 (MQ=255)

                 cGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGa                                                                        <  1:6432066/93‑1 (MQ=255)

                 cGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGa                                                                        <  1:5134122/93‑1 (MQ=255)

                  gTGCTATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGAt                                                                       <  1:7269870/93‑1 (MQ=255)

                  gTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGAt                                                                       >  2:1578265/1‑93 (MQ=255)

                  gTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGAt                                                                       <  2:4222397/93‑1 (MQ=255)

                      gATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTaaa                                                                   <  1:4525201/93‑1 (MQ=255)

                      gATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTaaa                                                                   >  2:322612/1‑93 (MQ=255)

                      gATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTaaa                                                                   <  2:7616701/93‑1 (MQ=255)

                       aTGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTaaaa                                                                  <  1:122942/93‑1 (MQ=255)

                       aTGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTaaaa                                                                  >  1:7731762/1‑93 (MQ=255)

                        tGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAg                                                                 >  2:3554756/1‑93 (MQ=255)

                        tGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAg                                                                 >  2:1380499/1‑93 (MQ=255)

                        tGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAg                                                                 >  1:7139255/1‑93 (MQ=255)

                           gTTCTTTCACACGAGAATCTTGACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCtt                                                              <  1:6302141/93‑1 (MQ=255)

                           gTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCtt                                                              >  2:3445857/1‑93 (MQ=255)

                           gTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCtt                                                              >  2:1801850/1‑93 (MQ=255)

                           gTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCtt                                                              >  1:3227958/1‑93 (MQ=255)

                           gTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCtt                                                              <  1:3607116/93‑1 (MQ=255)

                            ttCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTc                                                             >  1:2983311/1‑93 (MQ=255)

                             tCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCa                                                            >  2:5893208/1‑93 (MQ=255)

                               tttCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCAtc                                                          <  2:773182/93‑1 (MQ=255)

                               tttCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCAtc                                                          <  2:5667329/93‑1 (MQ=255)

                               tttCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCAtc                                                          <  2:3312392/93‑1 (MQ=255)

                                 tCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCAtctt                                                        >  2:255041/1‑93 (MQ=255)

                                 tCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCAtctt                                                        <  1:7810824/93‑1 (MQ=255)

                                    cacGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATcttctt                                                     <  2:3952774/93‑1 (MQ=255)

                                    cacGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATcttctt                                                     >  1:5330327/1‑93 (MQ=255)

                                    cacGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATcttctt                                                     <  2:719885/93‑1 (MQ=255)

                                     acGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCttt                                                    <  2:5775821/93‑1 (MQ=255)

                                     acGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCttt                                                    <  2:5430554/93‑1 (MQ=255)

                                       gagaATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCa                                                  <  1:2768230/93‑1 (MQ=255)

                                       gagaATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCa                                                  >  1:3646784/1‑93 (MQ=255)

                                          aaTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGAt                                               <  1:5919302/93‑1 (MQ=255)

                                          aaTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGAt                                               <  2:743619/93‑1 (MQ=255)

                                          aaTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGAt                                               >  1:4950691/1‑93 (MQ=255)

                                           aTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATa                                              <  1:2337782/93‑1 (MQ=255)

                                           aTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATa                                              <  2:2222931/93‑1 (MQ=255)

                                           aTCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATa                                              <  2:4490208/93‑1 (MQ=255)

                                            tCTTCGCTGTTGATATCCTCTTTCAGCCATTCCCGTAAAGCGCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAg                                             <  2:2324229/93‑1 (MQ=255)

                                             cTTCACTGTTGATATCCACTTTCAGCCATTCCTGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGc                                            <  2:5092487/93‑1 (MQ=255)

                                             cTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGc                                            <  1:5197772/93‑1 (MQ=255)

                                             cTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGc                                            <  2:3560904/93‑1 (MQ=255)

                                               tCACTGTTGATATCCACTTTCAGCCATTCACGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCtt                                          >  2:4839126/1‑93 (MQ=255)

                                                 aCTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCtttt                                        <  2:917010/93‑1 (MQ=255)

                                                 aCTGTTGATATCCACTATCCGCCCTTCCCGTAATTCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCGGATAGCtttt                                        <  1:7015312/93‑1 (MQ=255)

                                                   tGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTa                                      <  2:2813840/93‑1 (MQ=255)

                                                   tGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTa                                      <  1:6709163/93‑1 (MQ=255)

                                                   tGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTa                                      >  1:930911/1‑93 (MQ=255)

                                                   tGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTa                                      <  2:808649/93‑1 (MQ=255)

                                                     ttGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACg                                    >  1:2070209/1‑93 (MQ=255)

                                                       gATATCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGc                                  >  2:133736/1‑93 (MQ=255)

                                                         tatCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGTTAGCTTTTTACGGCAt                                <  2:7044025/93‑1 (MQ=255)

                                                         tatCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCAt                                >  1:6237399/1‑93 (MQ=255)

                                                         tatCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCAt                                <  2:1057273/93‑1 (MQ=255)

                                                          atCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATc                               >  2:6754773/1‑93 (MQ=255)

                                                          atCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATc                               <  1:3083395/93‑1 (MQ=255)

                                                          atCCACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATc                               >  2:2405696/1‑93 (MQ=255)

                                                            ccACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCgg                             >  2:1905171/1‑93 (MQ=255)

                                                             cACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGc                            <  1:5453809/93‑1 (MQ=255)

                                                             cACTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGc                            <  1:5037538/93‑1 (MQ=255)

                                                              aCTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCt                           >  2:2615274/1‑93 (MQ=255)

                                                               cTTTCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTc                          >  2:4587959/1‑93 (MQ=255)

                                                                tttCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGTATCGGCTcc                         >  1:6080117/1‑92 (MQ=255)

                                                                tttCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTcc                         >  2:5011019/1‑92 (MQ=255)

                                                                  tCAGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTcgcg                       >  2:6972414/1‑93 (MQ=255)

                                                                    aGCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCt                     <  1:5606199/93‑1 (MQ=255)

                                                                     gCCATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTg                    >  2:2782577/1‑93 (MQ=255)

                                                                      ccATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGc                   <  2:3100798/93‑1 (MQ=255)

                                                                       cATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGcc                  >  2:6217841/1‑93 (MQ=255)

                                                                       cATTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCCCGCTGcc                  >  2:633512/1‑93 (MQ=255)

                                                                        aTTCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCCCGCTGccc                 >  2:7358554/1‑92 (MQ=255)

                                                                         ttCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCCCGCTGccct                >  1:1438117/1‑91 (MQ=255)

                                                                          tCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCAtt               <  2:6131074/93‑1 (MQ=255)

                                                                          tCCCGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCAtt               <  2:5322316/93‑1 (MQ=255)

                                                                           cccGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCCttt              >  1:3521751/1‑93 (MQ=255)

                                                                            ccGTAAATCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCTGCTCGCGCCGCCCttgc             >  2:2973494/1‑91 (MQ=255)

                                                                                aaaTCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGa         <  2:4624525/93‑1 (MQ=255)

                                                                                 aaTCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGAt        <  1:4173034/93‑1 (MQ=255)

                                                                                 aaTCCCCCCGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGAt        <  1:4122450/93‑1 (MQ=255)

                                                                                     cccccGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGGTCggt    >  1:4625189/1‑91 (MQ=255)

                                                                                     cccccGCAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGtt    <  2:402807/93‑1 (MQ=255)

CAAACCATTGACCCGCTCGTGCGATGCGTTCTTTCACACGAGAATCTTCACTGTTGATATCCACTTTCAGCCATTCCCGTAAATCCCCCCACAATTCACTGAGATAGCGATTAAAAGCTTCATCTTCTTTCAGATAGCTTTTTACGGCATCGGCTCGCGCTGCCATTTCCGGATCGTTTT  >  NC_000913/4563819‑4563998

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