BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,274	T→C	L537P (CTC→CCC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,274	0	T	C	99.6%	914.3 / NA	279	L537P (CTC→CCC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (140/138); minor base A (1/0); total (141/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACACCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCg                                                                                         <  1:6034334/93‑1 (MQ=255)

  gTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCggg                                                                                       <  1:1074544/93‑1 (MQ=255)

  gTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCggg                                                                                       >  1:2703652/1‑93 (MQ=255)

   ttATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgc                                                                                      >  2:2798870/1‑93 (MQ=255)

   ttATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgc                                                                                      <  1:4298677/93‑1 (MQ=255)

    tATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgcg                                                                                     >  2:5004617/1‑93 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   >  1:1439149/1‑93 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   >  2:4626092/1‑93 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   >  2:4259725/1‑93 (MQ=255)

         aGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAt                                                                                <  1:1844945/93‑1 (MQ=255)

         aGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAt                                                                                >  1:1520004/1‑93 (MQ=255)

              aaGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACCTCAaca                                                                           >  1:3957728/1‑93 (MQ=255)

              aaGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAaca                                                                           <  2:1063054/93‑1 (MQ=255)

                gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACt                                                                         <  2:4441241/93‑1 (MQ=255)

                   gAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTa                                                                      <  2:2178845/93‑1 (MQ=255)

                   gAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTa                                                                      >  2:5422901/1‑93 (MQ=255)

                      aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   >  2:3031572/1‑93 (MQ=255)

                      aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   <  2:4228666/93‑1 (MQ=255)

                         cACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAg                                                                >  2:6369816/1‑93 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               <  2:2131644/93‑1 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               >  1:5212036/1‑93 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            >  2:6230184/1‑93 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            >  2:6057848/1‑93 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            >  1:1816853/1‑93 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            >  1:6333744/1‑93 (MQ=255)

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           >  2:5560153/1‑93 (MQ=255)

                                tATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGc                                                         >  2:5111331/1‑93 (MQ=255)

                                   tACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCAt                                                      <  1:2159285/93‑1 (MQ=255)

                                    aCAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                     <  2:4795670/93‑1 (MQ=255)

                                     cAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCa                                                    >  2:5046340/1‑93 (MQ=255)

                                     cAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCa                                                    <  1:660356/93‑1 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  2:3804095/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  1:2146250/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  1:4466025/1‑93 (MQ=255)

                                        ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                                 <  1:5346266/93‑1 (MQ=255)

                                        ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                                 <  1:1596365/93‑1 (MQ=255)

                                        ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                                 >  2:5211938/1‑93 (MQ=255)

                                        ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                                 >  1:4576312/1‑93 (MQ=255)

                                          tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACtt                                               <  1:5875605/93‑1 (MQ=255)

                                           cGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTa                                              <  1:2740169/93‑1 (MQ=255)

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            <  2:35548/93‑1 (MQ=255)

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            <  1:1987904/93‑1 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           >  2:6001351/1‑93 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           >  1:5752210/1‑93 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           >  1:5901206/1‑93 (MQ=255)

                                               tataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                          <  1:6173253/93‑1 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         <  2:6676519/93‑1 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         <  1:3607555/93‑1 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        <  2:3116748/93‑1 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        <  1:4250474/93‑1 (MQ=255)

                                                  aaCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa                                       <  1:2290578/93‑1 (MQ=255)

                                                  aaCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa                                       <  1:5105124/93‑1 (MQ=255)

                                                  aaCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa                                       >  1:5052116/1‑93 (MQ=255)

                                                   aCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGaa                                      <  2:395123/93‑1 (MQ=255)

                                                     aaaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACg                                    >  1:3871553/1‑93 (MQ=255)

                                                     aaaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACg                                    >  2:2374845/1‑93 (MQ=255)

                                                      aaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGa                                   <  1:3365328/93‑1 (MQ=255)

                                                      aaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGa                                   <  2:1131170/93‑1 (MQ=255)

                                                      aaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGa                                   <  2:2021438/93‑1 (MQ=255)

                                                        cgcgGCAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACTCTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGt                                 <  2:4098590/93‑1 (MQ=255)

                                                         gcgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGAGGTAGCCATCAAACTTATTGGGGAACGAGTa                                >  1:4845869/1‑93 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  2:5640851/93‑1 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  1:2414061/93‑1 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              <  1:6303297/93‑1 (MQ=255)

                                                            ggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATg                             <  1:2800904/93‑1 (MQ=255)

                                                            ggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATg                             >  1:5003946/1‑93 (MQ=255)

                                                            ggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATg                             >  1:1355225/1‑93 (MQ=255)

                                                              aaaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           <  2:5169767/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  2:224814/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         >  2:5875893/1‑93 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  1:1924791/93‑1 (MQ=255)

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        >  2:803313/1‑93 (MQ=255)

                                                                  tACATCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                       <  2:6464958/93‑1 (MQ=255)

                                                                  tACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                       >  2:2809960/1‑93 (MQ=255)

                                                                   aCAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGa                      >  2:2877962/1‑93 (MQ=255)

                                                                     aaCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                    >  2:3407257/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGCACGAGTAATGTTGATGAGCa                   >  2:3612910/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                   >  1:856214/1‑93 (MQ=255)

                                                                        tCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAt                 >  1:3058507/1‑93 (MQ=255)

                                                                        tCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAt                 <  1:5563845/93‑1 (MQ=255)

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               >  1:791097/1‑93 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              <  1:1103952/93‑1 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              <  1:487572/93‑1 (MQ=255)

                                                                            cGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCa             <  2:349685/93‑1 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            <  2:905710/93‑1 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            >  2:6233881/1‑93 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          <  2:5271951/93‑1 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          >  1:3151860/1‑93 (MQ=255)

                                                                                aCTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCt         >  1:367545/1‑93 (MQ=255)

                                                                                aCTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCt         <  2:3296559/93‑1 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  2:1774586/1‑93 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  1:170270/1‑93 (MQ=255)

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       >  1:5391323/1‑93 (MQ=255)

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       >  1:1213872/1‑93 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      >  2:1496399/1‑93 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      <  1:4764894/93‑1 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      <  1:2864108/93‑1 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      <  2:4334537/93‑1 (MQ=255)

                                                                                     gCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAtt    >  1:4686244/1‑93 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  >  2:2547952/1‑93 (MQ=255)

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

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