BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,542,476	C→T	6.6%	A147V (GCG→GTG)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,476	0	C	T	6.6%	656.0 / 29.0	182	A147V (GCG→GTG)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand): ref base C (92/78); new base T (5/7); total (97/85)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 183 total aligned reads displayed.

GCAGGCCTATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTT  >  NC_000913/4542378‑4542573

gcagGCCTATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGc                                                                                                   >  1:4020656/1‑99 (MQ=255)

   ggCCTATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGa                                                                                                >  1:3573011/1‑99 (MQ=255)

          cGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggt                                                                                         >  2:872092/1‑99 (MQ=255)

          cGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggt                                                                                         <  1:2803551/99‑1 (MQ=255)

           gCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                                                                                        >  2:2664208/1‑99 (MQ=255)

           gCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCgtggtg                                                                                        >  1:479233/1‑99 (MQ=255)

               tattCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATCTGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGa                                                                                    <  2:2853339/99‑1 (MQ=255)

               tattCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGTGGAGCGTGGTGCAGa                                                                                    >  2:1208518/1‑99 (MQ=255)

               tattCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGTGGAGCGTGGTGCAGa                                                                                    <  1:3955585/99‑1 (MQ=255)

               tattCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGa                                                                                    <  2:1109943/99‑1 (MQ=255)

                attCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGAt                                                                                   >  2:728035/1‑99 (MQ=255)

                 ttCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATa                                                                                  <  1:1369695/99‑1 (MQ=255)

                   cgcgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACt                                                                                <  2:3148519/99‑1 (MQ=255)

                    gcgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGTGGAGCGTGGTGCAGATACTc                                                                               >  2:2813261/1‑99 (MQ=255)

                     cgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCg                                                                              <  1:426408/99‑1 (MQ=255)

                     cgATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCg                                                                              <  1:4245215/99‑1 (MQ=255)

                      gATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGTGGAGCGTGGTGCAGATACTCGt                                                                             <  2:1830504/99‑1 (MQ=255)

                       aTGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGtt                                                                            <  1:1107153/99‑1 (MQ=255)

                       aTGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGtt                                                                            >  1:1972365/1‑99 (MQ=255)

                         gCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTa                                                                          <  1:2276027/99‑1 (MQ=255)

                           cGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAAt                                                                        <  1:457487/99‑1 (MQ=255)

                           cGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAAt                                                                        <  2:3532029/99‑1 (MQ=255)

                            ggTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCTGAGCGTGGTGCAGATACTCGTTTAAtt                                                                       <  2:2933973/99‑1 (MQ=255)

                             gTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTc                                                                      <  1:1786207/99‑1 (MQ=255)

                             gTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTc                                                                      >  2:3741471/1‑99 (MQ=255)

                              tATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCa                                                                     <  1:4104159/99‑1 (MQ=255)

                               aTTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAg                                                                    <  1:1414594/99‑1 (MQ=255)

                                ttGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGTGGAGCCTGGTGCAGATACTCGTTTAATTCAgg                                                                   <  1:1380176/99‑1 (MQ=255)

                                    aGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTa                                                               >  1:3154/1‑99 (MQ=255)

                                     gCTTGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAt                                                              >  1:3147907/1‑99 (MQ=255)

                                      cTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtc                                                             <  1:622203/99‑1 (MQ=255)

                                      cTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtc                                                             >  1:2621587/1‑99 (MQ=255)

                                      cTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtc                                                             <  1:3205640/99‑1 (MQ=255)

                                      cTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtc                                                             <  2:4292714/99‑1 (MQ=255)

                                       tGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtct                                                            >  2:1767404/1‑99 (MQ=255)

                                        ggAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTAtctc                                                           <  1:3821090/99‑1 (MQ=255)

                                         gAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCg                                                          >  2:4064828/1‑99 (MQ=255)

                                         gAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCg                                                          >  2:2489022/1‑99 (MQ=255)

                                         gAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCg                                                          >  2:681047/1‑99 (MQ=255)

                                         gAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCg                                                          <  2:3319787/99‑1 (MQ=255)

                                          aaCCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCgg                                                         >  1:1835253/1‑99 (MQ=255)

                                           aCCGTAACGCAGACTCATCCTCCTATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCggg                                                        >  1:1182251/1‑99 (MQ=255)

                                           aCCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCggg                                                        >  2:813455/1‑99 (MQ=255)

                                              gTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCAt                                                     <  2:3855162/99‑1 (MQ=255)

                                              gTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCAt                                                     <  1:1015677/99‑1 (MQ=255)

                                              gTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCAt                                                     >  1:2164627/1‑99 (MQ=255)

                                               tAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATc                                                    >  1:1686178/1‑99 (MQ=255)

                                               tAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATc                                                    >  2:2377914/1‑99 (MQ=255)

                                                 aCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                  >  1:4295502/1‑99 (MQ=255)

                                                 aCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                  >  2:1886694/1‑99 (MQ=255)

                                                   gCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaaa                                                <  2:436/99‑1 (MQ=255)

                                                   gCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaaa                                                <  2:500372/99‑1 (MQ=255)

                                                   gCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaaa                                                <  2:1059917/99‑1 (MQ=255)

                                                   gCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaaa                                                <  2:3232687/99‑1 (MQ=255)

                                                    cAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAat                                               <  2:4184380/99‑1 (MQ=255)

                                                     aGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAata                                              >  1:2680765/1‑99 (MQ=255)

                                                     aGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAata                                              <  1:892009/99‑1 (MQ=255)

                                                     aGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAata                                              >  2:3486347/1‑99 (MQ=255)

                                                      gACTTATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAatat                                             <  2:1431960/99‑1 (MQ=255)

                                                      gACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAatat                                             >  1:3696856/1‑99 (MQ=255)

                                                       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                            <  1:708918/99‑1 (MQ=255)

                                                       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                            <  1:946039/99‑1 (MQ=255)

                                                       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                            >  1:4126236/1‑99 (MQ=255)

                                                        cTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTc                                           >  1:3398501/1‑99 (MQ=255)

                                                         tCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCg                                          <  1:3707443/99‑1 (MQ=255)

                                                           aTCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGcc                                        >  2:3417110/1‑99 (MQ=255)

                                                           aTCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGcc                                        >  1:3446142/1‑99 (MQ=255)

                                                            tCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCa                                       >  1:1518879/1‑99 (MQ=255)

                                                            tCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCa                                       >  2:3143702/1‑99 (MQ=255)

                                                            tCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCa                                       >  1:2370235/1‑99 (MQ=255)

                                                            tCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCa                                       <  2:2464975/99‑1 (MQ=255)

                                                               tCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATAc                                    <  1:1258978/99‑1 (MQ=255)

                                                                 atatGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACtg                                  <  1:749155/99‑1 (MQ=255)

                                                                 atatGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACtg                                  <  2:671211/99‑1 (MQ=255)

                                                                  tatGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACtgt                                 >  1:1908471/1‑99 (MQ=255)

                                                                  tatGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACtgt                                 >  1:542109/1‑99 (MQ=255)

                                                                  tatGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACtgt                                 >  1:2958381/1‑99 (MQ=255)

                                                                    tGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGc                               >  1:3093196/1‑99 (MQ=255)

                                                                            cATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATAcc                       >  1:248185/1‑99 (MQ=255)

                                                                             aTGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCg                      >  1:3855450/1‑99 (MQ=255)

                                                                               gCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTACCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGcc                    <  1:4195328/99‑1 (MQ=255)

                                                                                cTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCa                   >  1:1633326/1‑99 (MQ=255)

                                                                                cTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCa                   >  2:3323473/1‑99 (MQ=255)

                                                                                 ttGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAg                  <  2:2993744/99‑1 (MQ=255)

                                                                                    cGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTaa               <  1:2756350/99‑1 (MQ=255)

                                                                                     ggTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAt              >  1:838386/1‑99 (MQ=255)

                                                                                      gTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtg             >  2:4182394/1‑99 (MQ=255)

                                                                                      gTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtg             >  2:2272793/1‑99 (MQ=255)

                                                                                        tATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgct           <  2:476530/99‑1 (MQ=255)

                                                                                          tGAATTGGTGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgctgc         >  2:1323359/1‑99 (MQ=255)

                                                                                          tGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgctgc         >  2:1104554/1‑99 (MQ=255)

                                                                                          tGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgctgc         >  1:2400364/1‑99 (MQ=255)

                                                                                           gAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATgctgct        <  2:3972751/99‑1 (MQ=255)

                                                                                            aaTTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTc       >  1:1576018/1‑99 (MQ=255)

                                                                                            aaTTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTc       <  2:3081885/99‑1 (MQ=255)

                                                                                              ttGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGt     >  1:1329883/1‑99 (MQ=255)

                                                                                              ttGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGt     >  2:832981/1‑99 (MQ=255)

                                                                                              ttGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGt     >  1:2184181/1‑99 (MQ=255)

                                                                                                ggCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGttt   <  1:749564/99‑1 (MQ=255)

                                                                                                ggCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGttt   <  1:1149887/99‑1 (MQ=255)

GCAGGCCTATCGCATTATTCGCGATGCCGGTATTGAAGCTGGAACCGTAACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTT  >  NC_000913/4542378‑4542573

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