BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	568,690	T→A	15.4%	intergenic (+43/‑212)	emrE → / → ybcK	DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,690	0	T	A	15.4%	432.8 / 64.0	162	intergenic (+43/‑212)	emrE/ybcK	DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand): ref base T (69/68); new base A (15/10); total (84/78)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 187 total aligned reads displayed.

TTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATT  >  NC_000913/568591‑568789

ttGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatat                                                                                                     >  2:2572364/1‑100 (MQ=255)

ttGATTTGTGCCGGTGTGTTGATTATTAATTTAGTGTCACGAAGCACACCACATTAAAATAATTTGTTTCCAAACGACTAAAATATGGAGGCTCttatat                                                                                                     >  1:2568350/1‑100 (MQ=255)

 tGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatt                                                                                                    <  2:2700074/100‑1 (MQ=255)

  gATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattt                                                                                                   >  2:1588447/1‑100 (MQ=255)

  gATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattt                                                                                                   >  2:624174/1‑100 (MQ=255)

   aTTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGCAGGCTCttatattta                                                                                                  >  2:4205735/1‑100 (MQ=255)

      tgtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                               >  2:3932903/1‑100 (MQ=255)

       gtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                              >  2:369159/1‑100 (MQ=255)

       gtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                              >  2:3908485/1‑100 (MQ=255)

         gCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAg                                                                                            <  2:1741296/100‑1 (MQ=255)

          ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGc                                                                                           >  1:744667/1‑100 (MQ=255)

          ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGc                                                                                           >  1:4105988/1‑100 (MQ=255)

           cGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGcc                                                                                          <  1:2525637/100‑1 (MQ=255)

            ggTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCt                                                                                         >  1:2553394/1‑100 (MQ=255)

             gtgtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTc                                                                                        <  2:1921699/100‑1 (MQ=38)

              tgtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTcc                                                                                       <  1:3884728/100‑2 (MQ=38)

                tgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtt                                                                                     >  1:4206728/1‑100 (MQ=255)

                 gtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGttt                                                                                    <  2:1169923/100‑1 (MQ=255)

                  ttGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtttt                                                                                   <  1:704877/100‑1 (MQ=255)

                  ttGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGttt                                                                                    <  1:4023486/99‑1 (MQ=255)

                   tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                  <  2:1036163/100‑1 (MQ=255)

                   tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                  <  2:1126058/100‑1 (MQ=255)

                   tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                  >  2:788850/1‑100 (MQ=255)

                   tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                  >  1:233681/1‑100 (MQ=255)

                    gATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTAt                                                                                 <  2:1857260/100‑1 (MQ=255)

                      ttattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGc                                                                               >  1:1760802/1‑100 (MQ=255)

                       tattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCt                                                                              >  1:3401260/1‑100 (MQ=255)

                        attaATTTATTGTCACGAAGCAAACCACATTAAAATAATTTGTTTCTAAACGACAAAAATATGAAGGCTCTTATATTTATATGAACCACGTTTTATGc                                                                               >  1:2665862/1‑98 (MQ=255)

                         ttaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCttt                                                                            >  1:1966351/1‑100 (MQ=255)

                         ttaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCttt                                                                            >  2:2155789/1‑100 (MQ=255)

                          taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCtttt                                                                           <  1:116791/100‑1 (MQ=255)

                           aaTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCttttt                                                                          <  1:1217662/100‑1 (MQ=255)

                              ttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                       <  1:2746453/100‑1 (MQ=255)

                               tATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTa                                                                      >  1:1935508/1‑100 (MQ=255)

                               tATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTa                                                                      >  1:876027/1‑100 (MQ=255)

                                 ttGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAAt                                                                    <  2:2519109/100‑1 (MQ=255)

                                  tGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATg                                                                   <  1:3111587/100‑1 (MQ=255)

                                   gTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGt                                                                  >  2:3137844/1‑100 (MQ=255)

                                    tCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTc                                                                 <  2:4052774/100‑1 (MQ=255)

                                    tCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTCTTTGTTAATGTc                                                                 <  2:2801897/100‑1 (MQ=255)

                                       cGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCttt                                                              >  2:1086739/1‑100 (MQ=255)

                                        gAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTa                                                             <  2:2966435/100‑1 (MQ=255)

                                            cacaCCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttt                                                         <  2:2833715/100‑1 (MQ=255)

                                             acacCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAttttt                                                        >  2:206888/1‑100 (MQ=255)

                                              caccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttttt                                                       <  2:4140279/100‑1 (MQ=255)

                                              caccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttttt                                                       >  1:2716618/1‑100 (MQ=255)

                                              caccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttttt                                                       <  1:4233060/100‑1 (MQ=255)

                                                ccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                     <  1:1604143/100‑1 (MQ=255)

                                                   caTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTa                                                  >  1:3531296/1‑100 (MQ=255)

                                                    aTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTAt                                                 >  2:2138908/1‑100 (MQ=255)

                                                      tAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTc                                               >  1:1036510/1‑100 (MQ=255)

                                                       aaaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCt                                              >  1:2304126/1‑100 (MQ=255)

                                                        aaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                             >  2:332597/1‑100 (MQ=255)

                                                        aaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                             >  1:1639106/1‑100 (MQ=255)

                                                           taatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                          <  1:3915159/100‑1 (MQ=255)

                                                           taatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtg                                          <  2:3185439/100‑1 (MQ=255)

                                                             atTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgtg                                        <  2:3451378/100‑1 (MQ=255)

                                                             atTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgtg                                        >  2:444640/1‑100 (MQ=255)

                                                                 gTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTTTGCttt                                    >  1:3128297/1‑100 (MQ=255)

                                                                 gTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                    >  2:81391/1‑100 (MQ=255)

                                                                 gTTTCTAAACGACGAAAATATGGAGGCTCTTCTATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                    >  2:2874497/1‑100 (MQ=255)

                                                                  tttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                   >  1:2816307/1‑100 (MQ=255)

                                                                  tttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                   <  1:1438605/100‑1 (MQ=255)

                                                                   ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                   >  2:388217/1‑99 (MQ=255)

                                                                   ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCa                                  <  1:640858/100‑1 (MQ=255)

                                                                     cTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                <  2:1293274/100‑1 (MQ=255)

                                                                       aaaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAt                              >  2:875153/1‑100 (MQ=255)

                                                                       aaaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGAt                              >  1:4163749/1‑100 (MQ=255)

                                                                         aCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTa                            <  2:2747176/100‑1 (MQ=255)

                                                                         aCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTa                            <  2:2459746/100‑1 (MQ=255)

                                                                          cGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTAt                           >  1:470891/1‑100 (MQ=255)

                                                                            aCTAAAAGATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                         <  2:362485/100‑1 (MQ=255)

                                                                             cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                        >  2:1984654/1‑100 (MQ=255)

                                                                             cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                        >  2:1003279/1‑100 (MQ=255)

                                                                             cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                        >  2:3581148/1‑100 (MQ=255)

                                                                              gcataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAATATTATGTCg                       >  1:2051994/10‑100 (MQ=21)

                                                                               aaaaTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                      <  2:470893/100‑1 (MQ=255)

                                                                               aaaaTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTATTTTGTGCTTTCAAGATTATGGCGt                      >  2:535006/1‑100 (MQ=255)

                                                                                 taaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTaa                    >  2:4013516/7‑100 (MQ=37)

                                                                                  aaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTTAAGATTATGGCGTAAg                   >  1:31769/6‑100 (MQ=37)

                                                                                   tatGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGa                  >  1:2981964/1‑100 (MQ=255)

                                                                                   aaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGa                  >  2:3273829/5‑100 (MQ=37)

                                                                                    atGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaa                 >  1:2933029/1‑100 (MQ=255)

                                                                                    aacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATTTCGTAAGa                  >  1:2157193/4‑99 (MQ=37)

                                                                                     tGGAGGATCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATATATTATATTGTGCTTTCAAGATTATGGCGTAAGa                  <  1:1717175/98‑1 (MQ=255)

                                                                                     acgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTTCAGATTATGGCGTTAGa                  >  1:1011939/3‑98 (MQ=37)

                                                                                      ggAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaaa               >  1:3028615/1‑100 (MQ=255)

                                                                                       gAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAAt              >  1:279751/1‑100 (MQ=255)

                                                                                         ggCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTg            >  2:3802936/1‑100 (MQ=255)

                                                                                          gCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc           <  1:1542986/100‑1 (MQ=255)

                                                                                           ctctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCa          >  1:3591786/1‑100 (MQ=255)

                                                                                            tctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc           <  1:4139459/98‑1 (MQ=255)

                                                                                             ctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAAt        >  2:761190/1‑100 (MQ=255)

                                                                                                 tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa    >  1:2773254/1‑100 (MQ=255)

                                                                                                 tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa    <  2:4163771/100‑1 (MQ=255)

                                                                                                 tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa    >  2:3760917/1‑100 (MQ=255)

                                                                                                 tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa    >  2:1561809/1‑100 (MQ=255)

                                                                                                  atttatatGAGCCTCGTTTTATGCTTTTTTTTAATGTATTTATTTTTTATGTATTCTTTTGTGCTTTCAATATTATGGCGTAAGAAAATTGCAATACGat   >  1:1420529/1‑100 (MQ=255)

                                                                                                  atttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGat   <  1:2840168/100‑1 (MQ=255)

                                                                                                   tttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGatt  >  1:3260866/1‑100 (MQ=255)

TTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATT  >  NC_000913/568591‑568789

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