BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,170	T→C	I502I (ATT→ATC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,170	0	T	C	99.5%	714.5 / ‑3.8	217	coding (1506/2637 nt)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (116/100); minor base . (1/0); total (117/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

Only 100 of 224 total aligned reads displayed.

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATATAACATACAGTCGAATGAATGGCTACAACATCGa                                                                                         >  1:4571280/1‑93 (MQ=255)

 aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaa                                                                                        >  2:1873725/1‑93 (MQ=255)

 aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaa                                                                                        >  2:3080669/1‑93 (MQ=255)

  atatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaaa                                                                                       <  2:796775/93‑1 (MQ=255)

   tatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAac                                                                                      >  1:2670738/1‑93 (MQ=255)

    atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                                     >  2:3043594/1‑93 (MQ=255)

     ttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                                    >  1:2774813/1‑93 (MQ=255)

         gTTAGTGGGTTCCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                <  1:3156328/93‑1 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                >  2:4283310/1‑93 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                >  2:1397416/1‑93 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACTCAGGAc                                                                               <  2:4668699/93‑1 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               >  2:58919/1‑93 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               <  2:2425802/93‑1 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              <  1:620138/93‑1 (MQ=255)

            aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACgg                                                                             >  1:898546/1‑93 (MQ=255)

            aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACgg                                                                             >  2:1952030/1‑93 (MQ=255)

             gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGa                                                                            >  2:3281734/1‑93 (MQ=255)

               gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGt                                                                          >  2:4336969/1‑93 (MQ=255)

                ggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                                         >  1:1394785/1‑93 (MQ=255)

                  ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                                       >  2:4312958/1‑93 (MQ=255)

                   tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                                      <  2:1652228/93‑1 (MQ=255)

                    aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                     <  2:2097621/93‑1 (MQ=255)

                    aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                     <  2:2643705/93‑1 (MQ=255)

                     ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCa                                                                    <  1:2185631/93‑1 (MQ=255)

                       gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                                  <  2:1423382/93‑1 (MQ=255)

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                                                 <  1:4369884/93‑1 (MQ=255)

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                                                 >  1:2069544/1‑93 (MQ=255)

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAAGGAGTTATTCAGGt                                                                 >  1:3909978/1‑93 (MQ=255)

                         tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                                                >  2:3704091/1‑93 (MQ=255)

                          aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                                               <  1:1770920/93‑1 (MQ=255)

                            tCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAg                                                             <  1:2621923/93‑1 (MQ=255)

                              gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                                           <  1:1086647/93‑1 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        >  2:724929/1‑93 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        >  1:411431/1‑93 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        >  1:1986002/1‑93 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        <  1:4405108/93‑1 (MQ=255)

                                  aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                                       <  2:368495/93‑1 (MQ=255)

                                    cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                     <  1:4080465/93‑1 (MQ=255)

                                    cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                     <  2:3068753/93‑1 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                    >  1:3406818/1‑93 (MQ=255)

                                      gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                   >  2:4180333/1‑93 (MQ=255)

                                      gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                   >  2:1913384/1‑93 (MQ=255)

                                      gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                   >  2:127673/1‑93 (MQ=255)

                                      gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                   <  1:1241378/93‑1 (MQ=255)

                                        tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                 >  2:3402467/1‑93 (MQ=255)

                                        aatTTTAACTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                 <  1:2249360/92‑1 (MQ=255)

                                         atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                <  2:397611/93‑1 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              >  1:3382472/1‑93 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAACCGGAAATTAAccaat                                              >  1:1263111/1‑90 (MQ=255)

                                             tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGATTTATTCAGGTTAAGCCGAAATTCACCGACTa                                            <  2:4116809/93‑1 (MQ=255)

                                              aaTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                                           <  2:407669/93‑1 (MQ=255)

                                               aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                          >  2:3684066/1‑93 (MQ=255)

                                                  tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACa                                       <  1:1622216/93‑1 (MQ=255)

                                                   cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                                      <  1:4033775/93‑1 (MQ=255)

                                                   cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                                      <  1:2300537/93‑1 (MQ=255)

                                                     cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAAAAcc                                    >  2:1016706/1‑93 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  1:234002/93‑1 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  1:1647586/93‑1 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  2:3384034/93‑1 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                <  1:829138/93‑1 (MQ=255)

                                                          acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCt                               >  2:367824/1‑93 (MQ=255)

                                                          acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCt                               <  1:4461603/93‑1 (MQ=255)

                                                           caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt                              >  2:4723606/1‑93 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat                            <  1:3444939/93‑1 (MQ=255)

                                                              caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           >  1:2952457/1‑93 (MQ=255)

                                                              caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           >  2:1591282/1‑93 (MQ=255)

                                                              caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           >  1:2983269/1‑93 (MQ=255)

                                                              caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           >  1:1549127/1‑93 (MQ=255)

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          <  1:4272533/93‑1 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         >  2:302039/1‑93 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         >  2:4579024/1‑93 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         >  2:4089008/1‑93 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         >  1:3569228/1‑93 (MQ=255)

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      >  2:4363802/1‑93 (MQ=255)

                                                                    gTCGCATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTCCAACCTCGCTTATAACAAAc                     <  1:3603710/93‑1 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     <  2:3514713/93‑1 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     <  1:2486486/93‑1 (MQ=255)

                                                                      cGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgc                   <  1:3309111/93‑1 (MQ=255)

                                                                       gaatgaatGGCTACATC‑‑‑GAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGa               >  2:3105458/1‑93 (MQ=255)

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  >  1:3748396/1‑93 (MQ=255)

                                                                       gaatgaatGGCT‑‑‑ACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGa               >  2:1311812/1‑93 (MQ=255)

                                                                        aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                 >  1:1446024/1‑93 (MQ=255)

                                                                         atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCggg                <  2:3763547/93‑1 (MQ=255)

                                                                           gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaa              <  2:2729368/93‑1 (MQ=255)

                                                                             atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt            <  1:4641839/93‑1 (MQ=255)

                                                                                gCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc         >  2:4542192/1‑93 (MQ=255)

                                                                                gCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc         <  2:688812/93‑1 (MQ=255)

                                                                                gCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc         <  1:3103648/93‑1 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  1:505964/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  2:433610/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       <  2:1060924/93‑1 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  2:2811926/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       <  2:2468402/93‑1 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      >  1:3567747/1‑93 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      <  1:1347808/93‑1 (MQ=255)

                                                                                     aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc    <  2:2330371/93‑1 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  1:2140713/1‑93 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  1:2915670/1‑93 (MQ=255)

                                                                                       caTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc  >  2:4271757/1‑93 (MQ=255)

                                                                                       caTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc  >  2:4453716/1‑93 (MQ=255)

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

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