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.7%	1162.0 / NA	366	L537P (CTC→CCC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (172/193); minor base A (1/0); total (173/193)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCg                                                                                         <  2:88679/93‑1 (MQ=255)

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

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

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

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

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

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

         aGGTTAAGCCAAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAt                                                                                <  2:3403549/93‑1 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               >  2:5318968/1‑93 (MQ=255)

            ttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCaa                                                                             >  2:1146713/1‑93 (MQ=255)

             tAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAac                                                                            <  1:7732459/93‑1 (MQ=255)

             tAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAac                                                                            >  2:3531264/1‑93 (MQ=255)

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

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

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

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAacac                                                                          >  1:6660633/1‑93 (MQ=255)

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

                 ccGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTg                                                                        >  1:1767511/1‑93 (MQ=255)

                  cGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGt                                                                       <  2:1561851/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    <  2:4419206/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    >  1:1912868/1‑93 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    <  2:6218533/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    <  1:4972351/93‑1 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  <  2:5999211/93‑1 (MQ=255)

                        tCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGa                                                                 <  1:5952757/93‑1 (MQ=255)

                         cACCGACTATTACAACCTCGCTTATAACAATCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTgcg                                                                <  2:634754/93‑3 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             <  2:161399/93‑1 (MQ=255)

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

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           <  1:1006773/93‑1 (MQ=255)

                               cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAg                                                          <  2:4928433/93‑1 (MQ=255)

                               cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAg                                                          >  1:582428/1‑93 (MQ=255)

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

                                tATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGc                                                         <  1:2584987/93‑1 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        <  1:7887462/93‑1 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  1:3402529/1‑93 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  2:2006746/1‑93 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  2:3032620/1‑93 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        <  2:4627452/93‑1 (MQ=255)

                                  ttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCa                                                       <  2:1352982/93‑1 (MQ=255)

                                  ttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCa                                                       <  2:5765396/93‑1 (MQ=255)

                                  ttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCa                                                       >  2:472845/1‑93 (MQ=255)

                                  ttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCa                                                       <  1:1047337/93‑1 (MQ=255)

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

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

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

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

                                        ccTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCACACATCAACACTGTATTTGAGTGGTAGCCATCAAAc                                                 <  1:6618397/93‑1 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                >  2:2446717/1‑93 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                <  1:2130407/93‑1 (MQ=255)

                                          tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACtt                                               >  1:5592828/1‑93 (MQ=255)

                                          tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACtt                                               >  1:3739711/1‑93 (MQ=255)

                                            gCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAt                                             <  2:638505/93‑1 (MQ=255)

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

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            >  2:4556758/1‑93 (MQ=255)

                                             cTTATAACAAACACGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGATTGGTAGTCATCAAACTTAtt                                            >  2:7636156/1‑93 (MQ=255)

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

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:5294325/93‑1 (MQ=255)

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

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  1:3546189/93‑1 (MQ=255)

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

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:3313027/93‑1 (MQ=255)

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

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  2:1674210/93‑1 (MQ=255)

                                              ttATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTg                                           <  1:6487771/93‑1 (MQ=255)

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

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

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

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

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGAGGTAGCCATCAAACTTATTgggg                                        >  2:5413413/1‑93 (MQ=255)

                                                    cAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAAc                                     <  1:861494/93‑1 (MQ=255)

                                                    cAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAAc                                     >  1:523357/1‑93 (MQ=255)

                                                     aaaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACg                                    <  2:180149/93‑1 (MQ=255)

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

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

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

                                                             gAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGt                            <  1:2785490/93‑1 (MQ=255)

                                                              aaaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           >  2:7864003/1‑93 (MQ=255)

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

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

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

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        <  1:3747359/93‑1 (MQ=255)

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

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        <  1:2322376/93‑1 (MQ=255)

                                                                  tACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                       <  1:476669/93‑1 (MQ=255)

                                                                   aCAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGa                      <  2:2640464/93‑1 (MQ=255)

                                                                    cAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAg                     <  1:6890047/93‑1 (MQ=255)

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

                                                                       cTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                  <  2:7122679/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                              ttACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAgg           >  2:1496180/1‑93 (MQ=255)

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

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        <  2:2291254/93‑1 (MQ=255)

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       <  1:1447037/93‑1 (MQ=255)

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

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

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