BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,313	A→C	17.0%	E604A (GAG→GCG) ‡	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,313	0	A	C	17.0%	323.5 / 69.1	147	E604A (GAG→GCG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand): ref base A (62/59); new base C (8/17); total (70/77)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-03

Only 100 of 150 total aligned reads displayed.

TGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACC  >  NC_000913/805214‑805412

tGGCGGAGTTTACCCTGTCTCGCCGCGATCCAGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGa                                                                                                     >  1:2849204/1‑100 (MQ=255)

 ggCGGAGTTTACCCTGTCTCGCCGCGATTCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    <  1:3518421/100‑1 (MQ=255)

 ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    >  1:192669/1‑100 (MQ=255)

 ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    <  2:2711940/100‑1 (MQ=255)

 ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAg                                                                                                    >  2:2224150/1‑100 (MQ=255)

  gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGc                                                                                                   <  2:2169346/100‑1 (MQ=255)

      aGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAc                                                                                               <  2:463904/100‑1 (MQ=255)

       gTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACa                                                                                              >  1:1632172/1‑100 (MQ=255)

        tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACat                                                                                             <  1:2136210/100‑2 (MQ=255)

         ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTTCTAAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACaga                                                                                            >  1:1199711/1‑100 (MQ=255)

         ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACatt                                                                                            <  2:600047/100‑3 (MQ=255)

         ttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACaga                                                                                            <  2:2839956/100‑1 (MQ=255)

          tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAAaga                                                                                            >  1:1957040/1‑99 (MQ=255)

           aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcc                                                                                          <  1:3839949/100‑5 (MQ=255)

           aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcc                                                                                          <  1:2250302/100‑5 (MQ=255)

           aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAgg                                                                                          >  2:2968815/1‑100 (MQ=255)

            nccggTCTCGACTCGATCCCTGTCATGTTACCAGAAGCAAATCGACTGCTGAGCCCGACAATCAGCGTCTGGCGGGGAATGTCAGCGCGCTGACagagc                                                                                          >  2:3233209/5‑98 (MQ=255)

            cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattccc                                                                                         <  1:3030672/100‑6 (MQ=255)

            cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCACGCTGACattccc                                                                                         <  1:2629926/100‑6 (MQ=255)

            cccTGTCTCGCCGCGATCCCGCTTATTTTAGCAGATGTAAATCTACTGCTCATCTGTAAAATCAGCGTCTATCGCTCAATGTTATCGAGCTGACaga                                                                                            <  2:1867742/97‑1 (MQ=255)

             ccTGTCTCGCCGCGATCCCGGTTATGTTAGTAGAAGTAAAGCTTCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCTCGCTGACattcccc                                                                                        <  1:3275026/100‑7 (MQ=255)

             ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                        <  2:1681635/100‑1 (MQ=255)

             ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                        <  2:3838922/100‑1 (MQ=255)

              cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtgt                                                                                       >  2:1142333/1‑100 (MQ=255)

                 tctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTg                                                                                    <  2:3105439/100‑1 (MQ=255)

                  ctcGCCGCGATCCCGGGTATGTTAGCAGAAGTCAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTg                                                                                    >  1:3728332/1‑99 (MQ=255)

                   tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcg                                                                                  <  1:1716565/100‑1 (MQ=255)

                   tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcg                                                                                  <  1:782719/100‑1 (MQ=255)

                    cgccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcgc                                                                                 >  1:993903/1‑100 (MQ=255)

                     gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTTACAGAGGTGTTTgcgcg                                                                                <  2:2559694/100‑1 (MQ=255)

                     gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcgcg                                                                                >  1:3371067/1‑100 (MQ=255)

                            ttccgatcttGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAg                                                                         <  1:451669/91‑1 (MQ=255)

                             cccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGc                                                                        <  2:564818/100‑1 (MQ=255)

                              ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCa                                                                       <  2:1382318/100‑1 (MQ=255)

                                    aTGTTAGCAGAAGTAAAGCGACTGCTGATCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGc                                                                 <  1:1726108/100‑1 (MQ=255)

                                    aTGTTAGCAGAAGTAAAGCGACTGCTGAGTTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTgn                                                                 <  2:2347399/100‑2 (MQ=255)

                                     tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                >  1:2986714/1‑100 (MQ=255)

                                     tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                >  1:1252576/1‑100 (MQ=255)

                                     tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                >  2:1966170/1‑100 (MQ=255)

                                     tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCAGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGGGCATTAAGCAGATTGc                                                                 >  2:3460297/1‑99 (MQ=255)

                                      gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                               >  2:1952083/1‑100 (MQ=255)

                                      gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTg                                                               >  2:2630175/1‑100 (MQ=255)

                                           cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                           >  2:3689711/1‑99 (MQ=255)

                                             gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGAAAGAGGTGTTTGCGAGCATTAAGCAGATTGCTGGTCAGGa                                                        >  1:3571040/1‑100 (MQ=255)

                                                 tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                    <  1:2500150/100‑1 (MQ=255)

                                                  aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCata                                                   >  1:2181714/1‑100 (MQ=255)

                                                    aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATAtt                                                 <  1:3785099/100‑1 (MQ=255)

                                                    aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATAtt                                                 <  1:1674177/100‑1 (MQ=255)

                                                        aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                             <  2:3879067/100‑1 (MQ=255)

                                                         ctgctgCGCTGGAAAATCAGCGTCTGCCGCGGAATGCCAGCCAGCCGACAGAGGTCTTTGCGCACATTAAGCAGATTGCTTGTTAGGAGGATAGTGATcc                                            <  2:1723528/100‑1 (MQ=255)

                                                          tgctgAGCTTGAATATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                    <  1:2681054/91‑1 (MQ=255)

                                                           gctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                          >  2:1629618/1‑100 (MQ=255)

                                                            ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                         <  1:3430128/100‑1 (MQ=255)

                                                             tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                        >  1:1216721/1‑100 (MQ=255)

                                                              gAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAAGAGCATATTGATCCGCTGc                                       >  1:1259433/1‑100 (MQ=255)

                                                                gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTAGTCAGGAGCATATTGATCCGCTGCaa                                     >  1:863430/1‑100 (MQ=255)

                                                                 cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaaa                                    >  2:668165/1‑100 (MQ=255)

                                                                  tGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                   >  2:1379168/1‑100 (MQ=255)

                                                                   ggAACATCAGCGTCTGGCGGAGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                  >  2:462640/1‑100 (MQ=255)

                                                                   ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                  >  1:1790336/1‑100 (MQ=255)

                                                                          cAGCGTCTGGCGGGGAATGTCAGCGTGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                           <  1:1542428/100‑1 (MQ=255)

                                                                          aaGCGTCTGGCGGGGAATGTCAGCGAGATGACAGAGGTGTTTGCGCTCATTAAGCAGATTGCTGGGCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                           >  2:1295736/2‑100 (MQ=255)

                                                                           aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                          <  2:1059125/100‑1 (MQ=255)

                                                                            gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                         <  2:231445/100‑1 (MQ=255)

                                                                            gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                         >  2:962909/1‑100 (MQ=255)

                                                                            gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCAGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                         >  1:446619/1‑100 (MQ=255)

                                                                             cGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAAGAGCATATTGATCCGCTGCAAACTGAAATTGgca                        >  1:1694305/1‑100 (MQ=255)

                                                                              gTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGGTGGACAAGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                       >  1:1063368/1‑100 (MQ=255)

                                                                              gTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                       <  1:3326836/100‑1 (MQ=255)

                                                                              gTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcag                       <  2:3291508/100‑1 (MQ=255)

                                                                               tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                      >  2:2230328/1‑100 (MQ=255)

                                                                               tCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcaga                      >  1:2352934/1‑99 (MQ=255)

                                                                                cTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagca                     <  1:3685678/100‑1 (MQ=255)

                                                                                 tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                    >  2:442416/1‑100 (MQ=255)

                                                                                 tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                    <  1:896851/100‑1 (MQ=255)

                                                                                  ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGTATTTAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                   <  1:357561/100‑1 (MQ=255)

                                                                                  ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                   >  1:407178/1‑100 (MQ=255)

                                                                                  ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                   >  2:1129952/1‑100 (MQ=255)

                                                                                  ggCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                   >  1:273396/1‑100 (MQ=255)

                                                                                   gCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                  <  2:3597463/100‑1 (MQ=255)

                                                                                    cGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGt                 <  2:443938/100‑1 (MQ=255)

                                                                                     ggggAATGTCAGCGAGCTGACAGAGGAGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGCGCATATTGATCCGCTGCAAACTGAAATTGGCAGAATGGTc                >  2:3692485/1‑100 (MQ=255)

                                                                                     cacctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                >  1:3504710/7‑100 (MQ=255)

                                                                                       ggAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              <  2:2986687/100‑1 (MQ=255)

                                                                                       cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCCTATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCt               >  2:2454240/5‑99 (MQ=255)

                                                                                       cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              >  2:2920091/5‑100 (MQ=255)

                                                                                       cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              >  1:2794738/5‑100 (MQ=255)

                                                                                       cctctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa              >  1:1407614/5‑100 (MQ=255)

                                                                                         tctGTCAGCTCGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGCGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATg            >  1:944383/3‑100 (MQ=255)

                                                                                          aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTTATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc           <  2:1420905/100‑1 (MQ=255)

                                                                                          aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc           <  2:3264279/100‑1 (MQ=255)

                                                                                          aTGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc           <  1:3684774/100‑1 (MQ=255)

                                                                                           tGTCAGCGAGCTGACAGAGGTGATTGCGCGAATTAAGCAGATTGATGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTAAAAGCAAGGTCTATGCg          >  2:2114728/1‑100 (MQ=255)

                                                                                            gTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg         <  2:1390844/100‑1 (MQ=255)

                                                                                              cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       <  1:653694/100‑1 (MQ=255)

                                                                                              cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       <  1:1379183/100‑1 (MQ=255)

                                                                                              cAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg       >  2:632581/1‑100 (MQ=255)

                                                                                               aGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa      >  1:1958682/1‑100 (MQ=255)

                                                                                                gCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGTTGGTTATGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa     <  1:1219889/100‑1 (MQ=255)

                                                                                                   aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAGCTGAAATTAGCAGCATGGTCTATGCGGTGAAAc   >  2:886099/1‑99 (MQ=255)

TGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACC  >  NC_000913/805214‑805412

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