BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,801,157	T→C	S14P (TCT→CCT)	b1721m →	ankyrin repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,157	0	T	C	88.4%	493.2 / 44.1	172	S14P (TCT→CCT)	b1721m	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base T (10/10); new base C (77/75); total (87/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 181 total aligned reads displayed.

AATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACT  >  NC_000913/1801068‑1801245

aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGcctc                                                                                       >  1:2958288/1‑93 (MQ=255)

  ttaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGcctcct                                                                                     <  2:5766144/93‑1 (MQ=255)

    aatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCa                                                                                   <  1:6170036/93‑1 (MQ=255)

      tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                                                 <  1:313984/93‑1 (MQ=255)

      tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATa                                                                                 <  1:6120696/93‑1 (MQ=255)

        atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGa                                                                               >  1:5432409/1‑93 (MQ=255)

         tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                                              >  1:3103324/1‑93 (MQ=255)

         tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGAt                                                                              <  1:471717/93‑1 (MQ=255)

            tcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGc                                                                           >  1:3797938/1‑93 (MQ=255)

              cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAATCCTCCTCATAGATTGcac                                                                         >  2:4532282/1‑93 (MQ=255)

              cTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGcac                                                                         <  1:5830524/93‑1 (MQ=255)

               tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGcaca                                                                        >  1:4050156/1‑93 (MQ=255)

                cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACAt                                                                       <  2:1907694/93‑1 (MQ=255)

                  aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCg                                                                     >  1:2044940/1‑93 (MQ=255)

                   aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGa                                                                    <  2:4380777/93‑1 (MQ=255)

                   aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGa                                                                    >  1:2618981/1‑93 (MQ=255)

                   aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGa                                                                    <  2:5776121/93‑1 (MQ=255)

                     ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATa                                                                  <  1:1461691/93‑1 (MQ=255)

                     ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATa                                                                  <  1:4177063/93‑1 (MQ=255)

                      cTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAg                                                                 <  1:2675914/93‑1 (MQ=255)

                       tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGc                                                                <  2:4866696/93‑1 (MQ=255)

                       tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGc                                                                >  1:1002665/1‑93 (MQ=255)

                       tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGc                                                                <  1:6326482/93‑1 (MQ=255)

                         cAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGa                                                              <  2:1101749/93‑1 (MQ=255)

                          aGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                                             <  2:5104039/93‑1 (MQ=255)

                          aGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGAc                                                             >  2:5046181/1‑93 (MQ=255)

                           gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACa                                                            <  2:5697324/93‑1 (MQ=255)

                             tatatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATa                                                          <  1:2295901/93‑1 (MQ=255)

                               tatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATAcc                                                        <  1:4975526/93‑1 (MQ=255)

                                atataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCa                                                       <  1:884313/93‑1 (MQ=255)

                                 tataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCaa                                                      >  1:584652/1‑93 (MQ=255)

                                  ataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                                     >  1:3127587/1‑93 (MQ=255)

                                  ataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                                     <  2:6185575/93‑1 (MQ=255)

                                    aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAacac                                                   >  1:460360/1‑93 (MQ=255)

                                     aGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                                                  >  2:2068527/1‑93 (MQ=255)

                                     aGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                                                  >  1:6427752/1‑93 (MQ=255)

                                      ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                                                 <  1:4830168/93‑1 (MQ=255)

                                      ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                                                 <  2:1809563/93‑1 (MQ=255)

                                      ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                                                 >  2:1597923/1‑93 (MQ=255)

                                       gCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                                                <  2:2021635/93‑1 (MQ=255)

                                        cGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                                               <  1:563180/93‑1 (MQ=255)

                                        cGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                                               >  2:2547311/1‑93 (MQ=255)

                                         gAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                                              <  2:4054638/93‑1 (MQ=255)

                                         gAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                                              <  1:255772/93‑1 (MQ=255)

                                         gAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                                              <  2:3102671/93‑1 (MQ=255)

                                          aaCAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                                             <  1:2093513/93‑1 (MQ=255)

                                          aaCAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                                             <  2:3173727/93‑1 (MQ=255)

                                           aCAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                                            >  2:1141091/1‑93 (MQ=255)

                                           aCAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                                            <  2:4202046/93‑1 (MQ=255)

                                              gATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                                         >  2:5194987/1‑93 (MQ=255)

                                              gATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                                         >  2:3337026/1‑93 (MQ=255)

                                              gATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                                         >  1:718077/1‑93 (MQ=255)

                                              gATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                                         >  1:1023688/1‑93 (MQ=255)

                                              gATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                                         <  2:6532579/93‑1 (MQ=255)

                                               aTTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                                        >  2:2236954/1‑93 (MQ=255)

                                                ttATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                                       >  2:641096/1‑93 (MQ=255)

                                                 tATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                                      <  1:6260069/93‑1 (MQ=255)

                                                   tGAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTa                                    >  2:4604467/1‑93 (MQ=255)

                                                    gAGTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt                                   <  2:6624075/93‑1 (MQ=255)

                                                      gTCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTa                                 >  1:4732360/1‑93 (MQ=255)

                                                       tCAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTaa                                >  2:1520420/1‑93 (MQ=255)

                                                        cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt                               >  2:3693640/1‑93 (MQ=255)

                                                        cAAAACGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt                               <  1:3289820/93‑1 (MQ=255)

                                                         aaaaCGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATc                              <  1:5886753/93‑1 (MQ=255)

                                                          aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa                             >  2:3678371/1‑93 (MQ=255)

                                                          aaaCGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa                             <  2:5607859/93‑1 (MQ=255)

                                                          aaaCGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa                             >  2:1344820/1‑93 (MQ=255)

                                                             cGATATCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt                          <  2:5438673/93‑1 (MQ=255)

                                                               atatCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg                        <  1:5256677/93‑1 (MQ=255)

                                                               atatCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg                        >  2:3445678/1‑93 (MQ=255)

                                                                tatCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc                       >  2:973353/1‑93 (MQ=255)

                                                                tatCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc                       >  1:1269402/1‑93 (MQ=255)

                                                                 atCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg                      >  1:5727666/1‑93 (MQ=255)

                                                                 atCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg                      >  2:1722233/1‑93 (MQ=255)

                                                                 atCATTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg                      <  2:1019055/93‑1 (MQ=255)

                                                                    aTTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcgc                   >  1:2244796/1‑93 (MQ=255)

                                                                    aTTATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcgc                   >  1:529262/1‑93 (MQ=255)

                                                                     ttATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGcc                  >  2:6330271/1‑93 (MQ=255)

                                                                     ttATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGcc                  <  1:5817180/93‑1 (MQ=255)

                                                                     ttATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGcc                  <  2:5001844/93‑1 (MQ=255)

                                                                      tATCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCa                 <  1:1144499/93‑1 (MQ=255)

                                                                       aTCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAg                <  2:2652997/93‑1 (MQ=255)

                                                                        tCAGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGc               <  2:4844909/93‑1 (MQ=255)

                                                                          aGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCtc             >  1:393866/1‑93 (MQ=255)

                                                                          aGAACTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCtc             >  1:4381205/1‑93 (MQ=255)

                                                                            aaCTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCtcat           >  1:2791432/1‑93 (MQ=255)

                                                                             aCTCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCtcatc          >  2:2937730/1‑93 (MQ=255)

                                                                               tCATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAc        <  1:796174/93‑1 (MQ=255)

                                                                                cATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAcc       >  2:5830162/1‑93 (MQ=255)

                                                                                cATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAcc       <  2:2575742/93‑1 (MQ=255)

                                                                                cATTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCATATTAATCAATTTGCGCGCCAGCTCATCAcc       >  2:3049762/1‑93 (MQ=255)

                                                                                 aTTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAccc      <  2:5942375/93‑1 (MQ=255)

                                                                                 aTTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAccc      >  1:4649447/1‑93 (MQ=255)

                                                                                 aTTATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCAccc      <  1:5014207/93‑1 (MQ=255)

                                                                                  ttATAAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCt     >  1:4185166/1‑93 (MQ=255)

                                                                                   tataAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTa    <  2:1072662/93‑1 (MQ=255)

                                                                                   tataAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTa    >  1:6144393/1‑93 (MQ=255)

                                                                                     taAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACt  <  2:6038125/93‑1 (MQ=255)

                                                                                     taAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACt  <  2:2704888/93‑1 (MQ=255)

                                                                                     taAGCCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACt  <  2:289994/93‑1 (MQ=255)

AATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACT  >  NC_000913/1801068‑1801245

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