BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	139,413	T→A	I605F (ATT→TTT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,413	0	T	A	99.6%	995.0 / NA	281	I605F (ATT→TTT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (138/142); minor base C (1/0); total (139/142)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATG  >  NC_000913/139322‑139502

 tGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACg                                                                                         >  2:1764718/1‑93 (MQ=255)

  gCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGc                                                                                        >  2:4481271/1‑93 (MQ=255)

    cTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAGCGGCAGTGCCATTGGGTTGGCAAACGCCa                                                                                      <  1:5292846/93‑1 (MQ=255)

    cTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCa                                                                                      <  2:576521/93‑1 (MQ=255)

     tttGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCAc                                                                                     <  2:5000283/93‑1 (MQ=255)

     tttGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCAc                                                                                     <  1:343254/93‑1 (MQ=255)

        ggCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCCTTGGGTTGGCAACCCCCACtgc                                                                                  >  1:2096363/1‑91 (MQ=255)

           aTCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACg                                                                               >  2:2079185/1‑93 (MQ=255)

            tCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGa                                                                              >  1:6641482/1‑93 (MQ=255)

            tCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGa                                                                              <  1:7291945/93‑1 (MQ=255)

            tCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGa                                                                              >  2:2766240/1‑93 (MQ=255)

            tCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGAAAACGCCACTTCACgt                                                                              >  1:2312277/1‑92 (MQ=255)

             cTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGAt                                                                             <  1:2485150/93‑1 (MQ=255)

             cTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGAt                                                                             <  1:4368577/93‑1 (MQ=255)

              tttCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGAtt                                                                            >  2:2213333/1‑93 (MQ=255)

               ttCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGAttt                                                                           <  1:3247381/93‑1 (MQ=255)

                 cGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTgg                                                                         <  1:2222678/93‑1 (MQ=255)

                   gctgctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGAt                                                                       <  2:7201293/93‑1 (MQ=255)

                   gctgctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGAt                                                                       <  1:3982893/93‑1 (MQ=255)

                    ctgctCCAGCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATc                                                                      <  2:3982585/93‑1 (MQ=255)

                     tgctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCa                                                                     <  2:6535588/93‑1 (MQ=255)

                     tgctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCa                                                                     <  2:118157/93‑1 (MQ=255)

                     tgctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCa                                                                     >  2:3479401/1‑93 (MQ=255)

                     tgcgCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCa                                                                     <  1:6253349/93‑1 (MQ=255)

                      gctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAACCGCCACTTCACGATTTGGATCaa                                                                    >  1:3904428/1‑93 (MQ=255)

                      gctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCaa                                                                    >  1:5898302/1‑93 (MQ=255)

                      gctCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCaa                                                                    >  1:1290839/1‑93 (MQ=255)

                        tCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACg                                                                  >  2:762706/1‑93 (MQ=255)

                         ccATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACgg                                                                 >  1:6858412/1‑93 (MQ=255)

                         ccATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACgg                                                                 >  2:4585242/1‑93 (MQ=255)

                          cATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGa                                                                <  1:2586639/93‑1 (MQ=255)

                          cATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGa                                                                <  2:6421569/93‑1 (MQ=255)

                           aTCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGaa                                                               <  1:3788840/93‑1 (MQ=255)

                           aTCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGaa                                                               >  2:1888444/1‑93 (MQ=255)

                            tCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGaaa                                                              >  1:2238750/1‑93 (MQ=255)

                            tCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGaaa                                                              <  1:5540722/93‑1 (MQ=255)

                              gggTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATc                                                            <  2:3545155/93‑1 (MQ=255)

                               ggTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATcc                                                           <  2:6743119/93‑1 (MQ=255)

                               ggTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATcc                                                           <  1:7797342/93‑1 (MQ=255)

                                gTTGCCAGGACCACGCGGGATCTGTTTCGAAACAAACGGCAGGGCCATTGGGTTGGCAAACGCCACTTTACGAGTTGGATCAACGGAAATcct                                                          >  2:6997494/1‑92 (MQ=255)

                                gTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATccc                                                          <  1:5802153/93‑1 (MQ=255)

                                gTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATccc                                                          >  1:3713556/1‑93 (MQ=255)

                                gTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATccc                                                          <  1:1787900/93‑1 (MQ=255)

                                 ttGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCg                                                         <  1:5250799/93‑1 (MQ=255)

                                  tGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGc                                                        <  1:1671103/93‑1 (MQ=255)

                                  tGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGc                                                        <  2:3106061/93‑1 (MQ=255)

                                  tGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGc                                                        <  2:6156194/93‑1 (MQ=255)

                                  tGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGc                                                        >  1:4712161/1‑93 (MQ=255)

                                    caaagACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGccc                                                      <  2:1475246/89‑1 (MQ=255)

                                       ggACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCAt                                                   <  2:3649429/93‑1 (MQ=255)

                                       ggACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCAt                                                   >  1:2818460/1‑93 (MQ=255)

                                        gACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCAtt                                                  <  1:5221149/93‑1 (MQ=255)

                                         gcaaCGAGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTc                                                 <  1:7453378/90‑1 (MQ=255)

                                         aCCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTc                                                 <  2:7345845/93‑1 (MQ=255)

                                          ccACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCAATGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATTCCGCCCCATTCg                                                >  2:7047539/1‑93 (MQ=255)

                                           cACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                               <  2:1266618/93‑1 (MQ=255)

                                           cACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTAACGATTTGGATCAACGGAAATCCCGCCCCATTCGa                                               <  2:3294791/93‑1 (MQ=255)

                                             cgcgGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAAc                                             >  1:5653045/1‑93 (MQ=255)

                                               cgGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACAt                                           >  2:1049480/1‑93 (MQ=255)

                                                gggATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATc                                          <  2:1719418/93‑1 (MQ=255)

                                                 ggATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATcc                                         >  2:4307604/1‑93 (MQ=255)

                                                  gATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATccc                                        <  1:3069727/93‑1 (MQ=255)

                                                   aTCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATcccc                                       >  1:1152165/1‑93 (MQ=255)

                                                   aTCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATcccc                                       >  1:3628139/1‑93 (MQ=255)

                                                   aTCAGTTTCAAAAAAAACGGCAGTGCCATTGGGTTGGCAAGCGCCACTCCCGGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATcccc                                       <  2:5708544/93‑1 (MQ=255)

                                                    tCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGACCATCCCCa                                      >  1:1049951/1‑93 (MQ=255)

                                                    tCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCa                                      <  2:4848854/93‑1 (MQ=255)

                                                    tCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCa                                      >  2:4819837/1‑93 (MQ=255)

                                                    tCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCa                                      >  2:3818918/1‑93 (MQ=255)

                                                    tCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCa                                      <  2:6256655/93‑1 (MQ=255)

                                                     cAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAg                                     <  2:5937196/93‑1 (MQ=255)

                                                        tttCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGtt                                  >  1:6044483/1‑93 (MQ=255)

                                                        tttCGAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGtt                                  >  2:5359676/1‑93 (MQ=255)

                                                            gAAACAAACGGCAGTGCCATTGGTTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTAccc                              <  2:4113428/93‑1 (MQ=255)

                                                            gAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTAccc                              <  1:837324/93‑1 (MQ=255)

                                                            gAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTAccc                              >  1:4662143/1‑93 (MQ=255)

                                                            gAAACAAACGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTAccc                              <  2:3175605/93‑1 (MQ=255)

                                                             aaacaaacGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCg                             >  1:1504158/1‑93 (MQ=255)

                                                                caaacGGCAGTGCCATTGGGTTGGCTAACGCCACTTCACGTTTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGa                          <  1:4179086/93‑1 (MQ=255)

                                                                caaacGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGa                          >  2:2370995/1‑93 (MQ=255)

                                                                 aaacGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGaa                         <  1:3374124/93‑1 (MQ=255)

                                                                  aacGGCAGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAg                        >  2:6920982/1‑93 (MQ=255)

                                                                      gCACTGCCATTGGGTTGGCAAACGCCACTCAACGATTTGGATAAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTTCCCGGGAAGACCa                    <  2:7138777/93‑1 (MQ=255)

                                                                        aGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGc                  <  2:5188092/93‑1 (MQ=255)

                                                                        aGTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGc                  <  1:1601315/93‑1 (MQ=255)

                                                                         gTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCg                 >  1:987101/1‑93 (MQ=255)

                                                                         gTGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCg                 >  1:6684294/1‑93 (MQ=255)

                                                                          tGCCATTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGt                <  1:7112598/93‑1 (MQ=255)

                                                                              aTTGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTAccc            >  1:2469667/1‑93 (MQ=255)

                                                                               ttGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCt           >  2:243936/1‑93 (MQ=255)

                                                                               ttGGGTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCt           >  1:6486518/1‑93 (MQ=255)

                                                                                 gggTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGt         >  1:5212323/1‑93 (MQ=255)

                                                                                  ggTTGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGtt        <  1:6602423/93‑1 (MQ=255)

                                                                                    ttGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCa      <  2:4241139/93‑1 (MQ=255)

                                                                                    ttGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCa      <  1:5830106/93‑1 (MQ=255)

                                                                                    ttGGCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCa      >  1:5457256/1‑93 (MQ=255)

                                                                                      ttcAAACGCCTCTTCACGATTAGGATCAACGTAAATCCTGCCCAATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGa    <  2:5488289/91‑1 (MQ=255)

                                                                                      ggCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGa    <  2:3370240/93‑1 (MQ=255)

                                                                                      ggCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGa    <  2:6836020/93‑1 (MQ=255)

                                                                                      ggCAAACGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGa    <  1:2738324/93‑1 (MQ=255)

GTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAGGACCACGCGGGATCAGTTTCGAAACAAACGGCAGTGCCATTGGGTTGGCAATCGCCACTTCACGATTTGGATCAACGGAAATCCCGCCCCATTCGAACATCCCCAGGTTACCCGGGAAGACCAGCGTACCCTGTTCAGATG  >  NC_000913/139322‑139502

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