BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,463,677	C→T	G166G (GGC→GGT)	hofD →	leader peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,463,677	0	C	T	99.1%	1142.5 / NA	325	G166G (GGC→GGT)	hofD	leader peptidase
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (145/176); minor base G (3/0); total (149/176)
Fisher's exact test for biased strand distribution p-value = 9.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

GGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCA  >  NC_000913/3463587‑3463766

ggTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGgtgc                                                                                         <  1:1523679/93‑1 (MQ=255)

 gTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCa                                                                                        <  1:9645411/93‑1 (MQ=255)

 gTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCa                                                                                        <  1:9039219/93‑1 (MQ=255)

  tGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCAt                                                                                       >  2:6505001/1‑93 (MQ=255)

    cTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATgg                                                                                     >  2:6833584/1‑93 (MQ=255)

     tgtgGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGgt                                                                                    >  2:3926449/1‑93 (MQ=255)

     tgtgGTGTGTTTACTGGGGCGTCTGGGTAGATTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCAgggg                                                                                    >  2:8946088/1‑92 (MQ=255)

      gtggtgTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGgtg                                                                                   >  1:9926289/1‑93 (MQ=255)

      gtggtgTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGgtg                                                                                   <  1:264495/93‑1 (MQ=255)

       tggtgTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCCTGgtgc                                                                                  >  1:6672657/1‑92 (MQ=255)

        ggtgTGTTTGCTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGtgtg                                                                                 <  1:3097924/93‑1 (MQ=255)

         gtgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtgg                                                                                <  2:6177170/93‑1 (MQ=255)

         gtgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtgg                                                                                >  2:3635431/1‑93 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               <  1:2316826/93‑1 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               >  1:2920510/1‑93 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               <  1:6257527/93‑1 (MQ=255)

          tgtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               <  2:1010264/93‑1 (MQ=255)

          tgtgtTTACTGGGGCCGCTGGTTAGTTTGTCACAAACATGTATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggc                                                                               <  1:4860564/93‑1 (MQ=255)

           gtgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggct                                                                              >  1:5013345/1‑93 (MQ=255)

            tgtTTACTGTGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctg                                                                             <  2:4966910/93‑1 (MQ=255)

            tgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctg                                                                             >  2:1889345/1‑93 (MQ=255)

            tgtTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctg                                                                             >  1:8515800/1‑93 (MQ=255)

              tttACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGCGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGtggctggc                                                                           >  2:3853597/1‑93 (MQ=255)

                 tcTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCaaa                                                                        <  1:8372403/92‑1 (MQ=255)

                   tGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCaaaag                                                                      >  2:2008710/1‑91 (MQ=255)

                   tGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACg                                                                      >  1:4628739/1‑93 (MQ=255)

                    ggggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGt                                                                     <  1:8753701/93‑1 (MQ=255)

                    ggggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGt                                                                     >  1:8025364/1‑93 (MQ=255)

                    ggggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGt                                                                     >  1:9673983/1‑93 (MQ=255)

                     gggCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGtt                                                                    <  1:6246514/93‑1 (MQ=255)

                         gTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGGTGGCa                                                                >  1:7369055/1‑93 (MQ=255)

                          tCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGcaca                                                               >  1:5221275/1‑90 (MQ=255)

                           cTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAAt                                                              <  2:7772018/93‑1 (MQ=255)

                            tGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATg                                                             >  2:7242127/1‑93 (MQ=255)

                            tGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATg                                                             >  2:3274276/1‑93 (MQ=255)

                             ggTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGa                                                            <  1:7733597/93‑1 (MQ=255)

                             ggTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGa                                                            <  1:1369227/93‑1 (MQ=255)

                             ggTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGa                                                            <  1:4348456/93‑1 (MQ=255)

                              gTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGAt                                                           <  1:6130552/93‑1 (MQ=255)

                               ttAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATa                                                          >  1:10008695/1‑93 (MQ=255)

                               ttAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATa                                                          <  2:8189841/93‑1 (MQ=255)

                               ttAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATa                                                          <  2:6099238/93‑1 (MQ=255)

                                tAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGGTGCCAATGATAc                                                         >  1:1205356/1‑93 (MQ=255)

                                 tgTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAct                                                        <  2:8833583/92‑1 (MQ=255)

                                  gTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATActg                                                       >  1:2199025/1‑93 (MQ=255)

                                      gTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGa                                                   <  1:586264/93‑1 (MQ=255)

                                         acaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTg                                                >  1:6419118/1‑93 (MQ=255)

                                         acaAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTg                                                >  2:4129247/1‑93 (MQ=255)

                                          caAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGc                                               >  1:5941389/1‑93 (MQ=255)

                                           aaaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGcc                                              >  2:4637289/1‑93 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCCAACGTTGCCCATGATACTGCTGATTTCCt                                             >  2:9785292/1‑93 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             <  1:512274/93‑1 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             >  2:3041603/1‑93 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             <  2:1800275/93‑1 (MQ=255)

                                            aaGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCt                                             <  1:792551/93‑1 (MQ=255)

                                             agagGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTc                                            <  1:1550452/93‑1 (MQ=255)

                                             agagGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTc                                            <  2:4580815/93‑1 (MQ=255)

                                               agGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGc                                          >  1:914907/1‑93 (MQ=255)

                                               agGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGc                                          <  1:300540/93‑1 (MQ=255)

                                                gggATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCt                                         >  1:8910288/1‑93 (MQ=255)

                                                gggATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCt                                         <  1:4487062/93‑1 (MQ=255)

                                                 ggATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTg                                        <  1:9440720/93‑1 (MQ=255)

                                                  gATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTgg                                       <  1:421573/93‑1 (MQ=255)

                                                  gATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTgg                                       >  1:1663506/1‑93 (MQ=255)

                                                    ttGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGGTGCCAATGATACTGCTGATTGCCTCGCTgggg                                     >  2:9025168/1‑92 (MQ=255)

                                                     tGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTg                                    <  1:2510310/93‑1 (MQ=255)

                                                      gggCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGGCCATGATACTGCTGATTGCCTCGCTGGGTgg                                   >  1:822264/1‑93 (MQ=255)

                                                      gggCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTgg                                   >  2:2441020/1‑93 (MQ=255)

                                                        gCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCa                                 >  2:7479503/1‑93 (MQ=255)

                                                        gCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCa                                 >  2:9635683/1‑93 (MQ=255)

                                                        gCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCa                                 <  1:6891871/93‑1 (MQ=255)

                                                        gCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCa                                 <  2:2564460/93‑1 (MQ=255)

                                                          tccGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCAtt                               <  1:2228855/91‑1 (MQ=255)

                                                          tACGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCAtt                               >  2:1211320/1‑93 (MQ=255)

                                                           aCGGTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTg                              >  2:2886375/1‑93 (MQ=255)

                                                             ggTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGt                            <  2:6308993/93‑1 (MQ=255)

                                                              gTGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGtt                           >  1:5640728/1‑93 (MQ=255)

                                                               tGATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTa                          <  1:9335219/93‑1 (MQ=255)

                                                                gATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTAc                         >  2:5235715/1‑93 (MQ=255)

                                                                gATTTCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATAATGCTGATTGCCTCGCTGGGTAGCATTGGTTAc                         >  1:2009128/1‑93 (MQ=255)

                                                                   ttCAAGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGcc                      <  2:6623921/93‑1 (MQ=255)

                                                                    tCATGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCa                     <  2:4561926/93‑1 (MQ=255)

                                                                      aaGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATc                   >  2:1278100/1‑93 (MQ=255)

                                                                       aGCTACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCg                  >  2:840885/1‑93 (MQ=255)

                                                                          tACTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGttt               <  1:3787659/93‑1 (MQ=255)

                                                                           aCTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTc              <  2:5833028/93‑1 (MQ=255)

                                                                            cTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTca             >  1:3562411/1‑93 (MQ=255)

                                                                            cTGGCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTca             <  1:7213297/93‑1 (MQ=255)

                                                                              ggCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTcaca           <  1:45141/93‑1 (MQ=255)

                                                                              ggCCGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTcaca           <  2:4052133/93‑1 (MQ=255)

                                                                                ccGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAAc         >  2:2318338/1‑93 (MQ=255)

                                                                                 cGCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACt        <  1:5485183/93‑1 (MQ=255)

                                                                                  gCTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACtt       <  2:111951/93‑1 (MQ=255)

                                                                                   cTGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTc      <  2:4396226/93‑1 (MQ=255)

                                                                                    tGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCt     >  1:3445797/1‑93 (MQ=255)

                                                                                    tGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCt     <  1:4082505/93‑1 (MQ=255)

                                                                                    tGCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCt     <  2:5484403/93‑1 (MQ=255)

                                                                                     gCAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTg    >  2:7433536/1‑93 (MQ=255)

                                                                                      cAGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGc   >  2:487592/1‑93 (MQ=255)

                                                                                       aGGTGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCa  >  1:8792845/1‑93 (MQ=255)

GGTGCTGTGGTGTGTTTACTGGGGCGTCTGGTTAGTTTGTCACAAAGAGGGATTGGGCTACGGTGATTTCAAGCTACTGGCCGCTGCAGGCGCATGGTGTGGCTGGCAAACGTTGCCAATGATACTGCTGATTGCCTCGCTGGGTGGCATTGGTTACGCCATCGTTTCACAACTTCTGCA  >  NC_000913/3463587‑3463766

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