BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,739,740	A→G	M141T (ATG→ACG)	yiaJ ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,739,740	0	A	G	99.4%	504.8 / ‑4.6	157	M141T (ATG→ACG)	yiaJ	putative regulator
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (64/92); total (64/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

Only 100 of 168 total aligned reads displayed.

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

 aGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTcc                                                                                        <  2:306051/93‑1 (MQ=255)

  gTCCGGGTGACTAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTccc                                                                                       <  2:143281/93‑1 (MQ=255)

   tCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCg                                                                                      <  1:2161465/93‑1 (MQ=255)

     cGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGt                                                                                    <  2:4188388/93‑1 (MQ=255)

      gggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGtt                                                                                   <  1:353383/93‑1 (MQ=255)

      gggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGtt                                                                                   <  2:2410752/93‑1 (MQ=255)

       ggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtg                                                                                  >  2:2062376/1‑93 (MQ=255)

       ggTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtg                                                                                  <  2:2187208/93‑1 (MQ=255)

        gTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTtgt                                                                                 <  2:4247560/93‑1 (MQ=255)

          tgccAAACGCCATGTAGATCTTGCGCATTGCGGAACAGTAGAGCGGCATATGCTGCCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTgg                                                                               <  2:1547847/91‑1 (MQ=255)

          gACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTgg                                                                               <  1:3100675/93‑1 (MQ=255)

          gACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTgg                                                                               <  1:4532259/93‑1 (MQ=255)

            ccAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGt                                                                             <  2:2258358/93‑1 (MQ=255)

               aaCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTcc                                                                          <  2:1414645/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  1:1732440/93‑1 (MQ=255)

                aCGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCa                                                                         <  1:1586164/93‑1 (MQ=255)

                  gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       >  1:260415/1‑93 (MQ=255)

                  gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       >  2:4628307/1‑93 (MQ=255)

                  gCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGc                                                                       <  1:1612672/93‑1 (MQ=255)

                   ccATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCt                                                                      <  1:655264/93‑1 (MQ=255)

                    cATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCtt                                                                     <  1:4126780/93‑1 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCGGCTta                                                                    >  2:2394952/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    >  1:991336/1‑93 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    <  1:772110/93‑1 (MQ=255)

                     aTGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTta                                                                    <  2:2805785/93‑1 (MQ=255)

                      tGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTtat                                                                   >  1:2772738/1‑93 (MQ=255)

                          gATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAAt                                                               >  2:3017970/1‑93 (MQ=255)

                           aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATc                                                              >  1:2572843/1‑93 (MQ=255)

                           aTCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATc                                                              <  2:675796/93‑1 (MQ=255)

                              ttGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCaaa                                                           >  1:998342/1‑93 (MQ=255)

                                gCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAAt                                                         <  2:4080748/93‑1 (MQ=255)

                                 gccATTGCGGAACAGTAGAGCTGGATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATa                                                        <  2:3907728/92‑1 (MQ=255)

                                 cccATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATa                                                        <  2:1165695/93‑1 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      >  1:4415119/1‑93 (MQ=255)

                                   cATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  2:1164755/93‑1 (MQ=255)

                                   cATTGCGGAACAGTAAAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGc                                                      <  2:2450270/93‑1 (MQ=255)

                                     ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCACAATAGCGt                                                    >  2:1843232/1‑93 (MQ=255)

                                     ttGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGt                                                    <  1:3029144/93‑1 (MQ=255)

                                       gCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGa                                                  <  1:3219357/93‑1 (MQ=255)

                                        cGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGAt                                                 <  2:4306400/93‑1 (MQ=255)

                                            aCAGTAGAGCGGTATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTc                                             <  1:3476859/93‑1 (MQ=255)

                                            aCAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTc                                             >  2:1074361/1‑93 (MQ=255)

                                             cAGTAGAGTGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCt                                            <  1:3845596/93‑1 (MQ=255)

                                             cAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCt                                            <  2:3040476/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  1:1877714/1‑93 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  2:2716759/1‑93 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          <  2:97744/93‑1 (MQ=255)

                                               gTAGAGCGGCATATGCTGGCCAATATAGGCACGGGCTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTc                                          >  2:2507482/1‑93 (MQ=255)

                                                tAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcg                                         <  2:2633091/93‑1 (MQ=255)

                                                 agagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgc                                        <  2:1584412/93‑1 (MQ=255)

                                                 agagCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTcgc                                        <  2:1388022/93‑1 (MQ=255)

                                                    gCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGc                                     <  2:3292935/93‑1 (MQ=255)

                                                    gCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGc                                     <  1:2650035/93‑1 (MQ=255)

                                                     gggaatatGCTGCCCAATATAGGCATGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCt                                    <  2:1195341/89‑1 (MQ=255)

                                                      ggCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTg                                   >  2:1613893/1‑93 (MQ=255)

                                                      ggCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTg                                   <  1:3927811/93‑1 (MQ=255)

                                                       gCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCTGCTgg                                  >  1:2888020/1‑93 (MQ=255)

                                                       gCATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTgg                                  >  1:1397794/1‑93 (MQ=255)

                                                        cATATGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGga                                 >  1:3420855/1‑93 (MQ=255)

                                                         atatGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGCTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTGGCGGCTGgag                                >  1:1102190/1‑93 (MQ=255)

                                                           atGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGgagca                              >  2:2152494/1‑91 (MQ=255)

                                                            tGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAg                             >  2:3329419/1‑93 (MQ=255)

                                                            tGCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAg                             >  2:1736735/1‑93 (MQ=255)

                                                             gCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGt                            <  1:1437500/93‑1 (MQ=255)

                                                             gCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGt                            >  1:3007075/1‑93 (MQ=255)

                                                             gCTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGt                            >  1:4669850/1‑93 (MQ=255)

                                                              cTGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGtt                           <  2:1605724/93‑1 (MQ=255)

                                                               tGGCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTa                          >  2:4090128/1‑93 (MQ=255)

                                                                ggCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTaa                         >  1:4734836/1‑93 (MQ=255)

                                                                 gCCAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAAt                        >  2:1667689/1‑93 (MQ=255)

                                                                  ccAATATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATg                       >  1:4401277/1‑93 (MQ=255)

                                                                    aaTATAGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATTGt                     >  2:2415989/1‑93 (MQ=255)

                                                                      tataGGCGCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGttt                   <  2:1403707/93‑1 (MQ=255)

                                                                      tataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGGGAAGTTAATAGttt                   >  1:1124922/1‑93 (MQ=255)

                                                                      tataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGttt                   >  2:212436/1‑93 (MQ=255)

                                                                      tataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGCTAATGGttt                   >  2:4545188/1‑93 (MQ=255)

                                                                       ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  <  2:758273/93‑1 (MQ=255)

                                                                       ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  <  2:730977/93‑1 (MQ=255)

                                                                       ataGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTc                  <  1:1265158/93‑1 (MQ=255)

                                                                        tgggCGCGGGTTCGCAGCGTCCCGGTTGTGGGATCCAGCTTATAGATCAAAGTAGCGTGATCGTCTTCGCGGCTGGTGAAGTTAATGGTTTcc                 <  2:3601821/91‑2 (MQ=255)

                                                                        taGGCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCa                 <  1:2752654/93‑1 (MQ=255)

                                                                          ggCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCAcc               <  1:974623/93‑1 (MQ=255)

                                                                          ggCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCAcc               <  1:2058386/93‑1 (MQ=255)

                                                                          ggCACGGGTTCGCAGCGTACCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCAcc               <  2:2277741/93‑1 (MQ=255)

                                                                           gCACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCa              <  2:3421667/93‑1 (MQ=255)

                                                                            gaCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAg             <  1:2846819/92‑1 (MQ=255)

                                                                            cACGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAg             >  1:632349/1‑93 (MQ=255)

                                                                             aCGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATTGGTTCACCAGt            >  2:2943203/1‑93 (MQ=255)

                                                                              cGGGTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTg           >  2:4611386/1‑93 (MQ=255)

                                                                               gggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTgg          <  2:3663490/93‑1 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGCAAAGTTAATGGGTTCACCAGTGGc         >  2:2581076/1‑93 (MQ=255)

                                                                                ggTTCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGc         >  1:3560993/1‑93 (MQ=255)

                                                                                  ttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGa       >  2:2687677/1‑93 (MQ=255)

                                                                                  ttCGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGa       <  1:955512/93‑1 (MQ=255)

                                                                                    cGCAGCGTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATg     >  2:1512873/1‑93 (MQ=255)

                                                                                     gcagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGt    <  2:826227/93‑1 (MQ=255)

                                                                                      cagcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGtt   >  1:2880188/1‑93 (MQ=255)

                                                                                       agcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTc  <  1:2454946/93‑1 (MQ=255)

                                                                                       agcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTc  >  1:4687417/1‑93 (MQ=255)

                                                                                       agcgTCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTc  >  1:894255/1‑93 (MQ=255)

TAGTCCGGGTGACCAAACGCCATGTAGATCTTGCCCATTGCGGAACAGTAGAGCGGCATATGCTGGCCAATATAGGCACGGGTTCGCAGCATCCCGGTTGTGGGTTCCAGCTTATAAATCAAAATAGCGTGATCGTCTTCGCGGCTGGAGAAGTTAATGGTTTCACCAGTGGCGATGTTC  >  NC_000913/3739650‑3739829

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