BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,292	T→A	S101R (AGT→AGA)	mutL →	MutL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,292	0	T	A	98.8%	549.5 / ‑1.7	163	S101R (AGT→AGA)	mutL	MutL
Reads supporting (aligned to +/- strand): ref base T (0/2); new base A (85/76); total (85/78)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 173 total aligned reads displayed.

GCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAA  >  NC_000913/4395202‑4395381

gctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAAt                                                                                         <  2:2384030/93‑1 (MQ=255)

 ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATc                                                                                        <  2:1700365/93‑1 (MQ=255)

 ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATc                                                                                        >  2:3458677/1‑93 (MQ=255)

 ctggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATc                                                                                        >  2:4460299/1‑93 (MQ=255)

  tggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCa                                                                                       <  1:1554940/93‑1 (MQ=255)

  tggcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCa                                                                                       >  2:1889404/1‑93 (MQ=255)

   ggcTCGTCATGCCACAAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAg                                                                                      >  2:3333316/1‑93 (MQ=255)

    gcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGt                                                                                     <  2:1286350/93‑1 (MQ=255)

    gcTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGt                                                                                     <  2:2352786/93‑1 (MQ=255)

     cTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGGCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGtt                                                                                    <  2:2216548/93‑1 (MQ=255)

      tCGTCATGCCACCAGTAAAATCGCCTCTCTGGCCGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGCGGCGCTGGCGAGCATCAGTTc                                                                                   >  1:902244/1‑93 (MQ=255)

      tCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTc                                                                                   <  2:3009291/93‑1 (MQ=255)

      tCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTc                                                                                   <  2:1375744/93‑1 (MQ=255)

      tCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTc                                                                                   >  1:4301924/1‑93 (MQ=255)

       cGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCg                                                                                  >  2:4127476/1‑93 (MQ=255)

       cGTCATGCCACCAGTAAAATCGCACCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCg                                                                                  <  1:4632810/93‑1 (MQ=255)

        gTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCgg                                                                                 >  2:2418274/1‑93 (MQ=255)

         tCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGt                                                                                >  1:814180/1‑93 (MQ=255)

         tCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGt                                                                                <  2:3045879/93‑1 (MQ=255)

         tCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGt                                                                                >  2:3765597/1‑93 (MQ=255)

         tCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGt                                                                                >  1:3958201/1‑93 (MQ=255)

          cATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGtt                                                                               >  2:2698641/1‑93 (MQ=255)

          cATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGtt                                                                               >  1:1100993/1‑93 (MQ=255)

             gCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTcc                                                                            <  1:1756912/93‑1 (MQ=255)

               caccaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCg                                                                          >  1:1604070/1‑93 (MQ=255)

               caccaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCg                                                                          >  1:4268731/1‑93 (MQ=255)

                accaGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGc                                                                         <  2:4612503/93‑1 (MQ=255)

                 ccaGTAGAATCGCCTCTCTGGGCGGTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCAGcc                                                                        <  2:2499798/93‑1 (MQ=255)

                 ccaGTAAAAGCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGcc                                                                        <  2:2914657/93‑1 (MQ=255)

                  caGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCt                                                                       <  2:222162/93‑1 (MQ=255)

                  caGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCt                                                                       <  2:222163/93‑1 (MQ=255)

                     tAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGAc                                                                    <  1:2830668/93‑1 (MQ=255)

                     tAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGAc                                                                    >  1:2259262/1‑93 (MQ=255)

                     tAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGAc                                                                    >  1:2259261/1‑93 (MQ=255)

                     tAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGAc                                                                    <  1:1433939/93‑1 (MQ=255)

                     aaaaaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGTTGAGGCTCTGCCGAGAATCAGTTCGGTTTCCCGCCTGAc                                                                    <  2:2478908/92‑1 (MQ=255)

                      aaaaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACg                                                                   <  1:2577282/93‑1 (MQ=255)

                       aaaTCGGCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGc                                                                  >  2:4669124/1‑93 (MQ=255)

                        aaTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCt                                                                 <  1:3695651/93‑1 (MQ=255)

                         aTCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTc                                                                >  2:1162005/1‑93 (MQ=255)

                              ctctctGGAGGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTc                                                           >  1:4674299/1‑93 (MQ=255)

                              ctctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTc                                                           >  2:4577888/1‑93 (MQ=255)

                              ctctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTc                                                           >  2:1863116/1‑93 (MQ=255)

                              ctctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTc                                                           >  2:2958402/1‑93 (MQ=255)

                               tctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCa                                                          <  2:939447/93‑1 (MQ=255)

                               tctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCa                                                          >  1:3118299/1‑93 (MQ=255)

                               cctcgGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCa                                                          <  1:4066770/92‑1 (MQ=255)

                                ctctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCAc                                                         >  2:1922033/1‑93 (MQ=255)

                                 tctGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACg                                                        <  2:4401871/93‑1 (MQ=255)

                                   tGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCa                                                      >  1:2656124/1‑93 (MQ=255)

                                   tGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCa                                                      >  1:3838338/1‑93 (MQ=255)

                                    ggACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCAc                                                     >  2:2275347/1‑93 (MQ=255)

                                     gACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCAcc                                                    <  2:3115769/93‑1 (MQ=255)

                                       cGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGc                                                  >  2:1026620/1‑93 (MQ=255)

                                         aTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAg                                                <  1:67033/93‑1 (MQ=255)

                                         aTCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAg                                                >  1:2145352/1‑93 (MQ=255)

                                          tctcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGa                                               >  2:4285702/1‑93 (MQ=255)

                                           ctcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGaa                                              <  1:2662746/93‑1 (MQ=255)

                                            tcGAAGCCATTATCAGCCTGGGCTTTCGCGTTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAc                                             <  1:3901044/93‑1 (MQ=255)

                                            tcGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAc                                             <  2:2610290/93‑1 (MQ=255)

                                              gAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCACAAcag                                           >  1:2692056/1‑93 (MQ=255)

                                                aGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagca                                         <  1:1824772/93‑1 (MQ=255)

                                                aGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagca                                         <  1:1993982/93‑1 (MQ=255)

                                                aGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagca                                         <  2:2023000/93‑1 (MQ=255)

                                                 gCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagcag                                        <  2:6129/93‑1 (MQ=255)

                                                 gCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagcag                                        >  2:1524911/1‑93 (MQ=255)

                                                 gCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAcagcag                                        >  2:1066816/1‑93 (MQ=255)

                                                   cATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGa                                      >  1:1194184/1‑93 (MQ=255)

                                                    aTTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGaa                                     <  1:4413917/93‑1 (MQ=255)

                                                     ttATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAg                                    >  1:1759074/1‑93 (MQ=255)

                                                     ttATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAg                                    <  1:1392295/93‑1 (MQ=255)

                                                      tATCAGCCTGGGCTTGCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACGTCACGCACCGCAGAACAGCAGGAAGc                                   <  2:3274219/93‑1 (MQ=255)

                                                        tCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAAAAGAaagaagacc                                 >  2:3213414/1‑90 (MQ=255)

                                                         cAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTg                                >  2:737278/1‑93 (MQ=255)

                                                          atccTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTgg                               <  1:1466594/91‑1 (MQ=255)

                                                           gCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGc                              <  1:4554434/93‑1 (MQ=255)

                                                           gCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGc                              <  1:4165994/93‑1 (MQ=255)

                                                           gCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGc                              <  2:3644521/93‑1 (MQ=255)

                                                            ccTGGGCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCGGAAAAGCAGCAAGCCTGGCa                             >  2:1131345/1‑93 (MQ=255)

                                                               gggCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGc                          >  1:3945949/1‑93 (MQ=255)

                                                                 gCTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCt                        >  2:894621/1‑93 (MQ=255)

                                                                  cTTTCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTa                       <  2:2437489/93‑1 (MQ=255)

                                                                   tttCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTAt                      <  1:1729569/93‑1 (MQ=255)

                                                                     tCGCGGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGc                    <  2:4034957/93‑1 (MQ=255)

                                                                      cgcgGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGcc                   <  2:734213/93‑1 (MQ=255)

                                                                      cgcgGTGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGcc                   <  1:3059106/93‑1 (MQ=255)

                                                                           tGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAgg              >  1:547237/1‑93 (MQ=255)

                                                                           tGAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAgg              >  2:502164/1‑93 (MQ=255)

                                                                            gAGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAggg             >  2:39586/1‑93 (MQ=255)

                                                                             aGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgc            <  2:3629928/93‑1 (MQ=255)

                                                                             aGGCGCTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGgc            <  2:2179483/93‑1 (MQ=255)

                                                                                 gcTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGa        <  2:419043/93‑1 (MQ=255)

                                                                                  cTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGat       <  1:1848656/93‑1 (MQ=255)

                                                                                  cTGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGat       <  1:3607351/93‑1 (MQ=255)

                                                                                   tGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGata      <  2:3135080/93‑1 (MQ=255)

                                                                                   tGGCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGata      <  1:71749/93‑1 (MQ=255)

                                                                                    ggCGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCCGGCAGGCCTATGCCGAAGGGCGCGatat     <  1:3200678/93‑1 (MQ=255)

                                                                                      cGAGAATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGa   <  1:2355331/93‑1 (MQ=255)

                                                                                       gagaATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGaa  >  2:4384887/1‑93 (MQ=255)

                                                                                       gAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGaa  <  2:1556626/93‑1 (MQ=255)

GCTGGCTCGTCATGCCACCAGTAAAATCGCCTCTCTGGACGATCTCGAAGCCATTATCAGCCTGGGCTTTCGCGGTGAGGCGCTGGCGAGTATCAGTTCGGTTTCCCGCCTGACGCTCACTTCACGCACCGCAGAACAGCAGGAAGCCTGGCAGGCCTATGCCGAAGGGCGCGATATGAA  >  NC_000913/4395202‑4395381

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