BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,563,585	(C)_5→6	coding (1320/1404 nt)	b2451 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,563,580	1	.	C	98.4%	386.8 / ‑1.0	128	coding (1325/1404 nt)	b2451	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (59/66); total (61/67)
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.58e-01

Only 100 of 147 total aligned reads displayed.

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

gCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑t                                                                                           >  2:3341484/1‑93 (MQ=255)

gCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑t                                                                                           <  2:1476658/93‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑c                                                                                       <  1:2694948/93‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑c                                                                                       <  1:2961940/93‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑c                                                                                       <  1:3850909/93‑1 (MQ=255)

  gAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑c                                                                                       <  1:461758/93‑1 (MQ=255)

   aaTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑cc                                                                                      >  1:2752318/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑ccc                                                                                     >  2:4069247/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑ccc                                                                                     >  2:409971/1‑93 (MQ=255)

    aTACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑ccc                                                                                     >  1:450756/1‑93 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑cccc                                                                                    >  2:4262589/1‑93 (MQ=255)

     tACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑cccc                                                                                    >  2:2220400/1‑93 (MQ=255)

       ccGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑c‑cccccg                                                                                  <  1:1894969/93‑3 (MQ=255)

           aTTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCcgcgcg                                                                               >  2:221505/1‑93 (MQ=255)

           aTTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCcgcgcg                                                                               >  2:3449407/1‑93 (MQ=255)

            ttCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCgg                                                                              >  1:2379768/1‑93 (MQ=255)

            ttCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCgg                                                                              >  2:4024216/1‑93 (MQ=255)

             tCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGa                                                                             >  1:3023516/1‑93 (MQ=255)

               aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACa                                                                           >  2:1042694/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACaa                                                                          >  2:1423733/1‑93 (MQ=255)

                ggaaggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACaa                                                                          >  1:1398619/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAACAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCa                                                                       >  1:1123510/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCa                                                                       >  1:3401166/1‑93 (MQ=255)

                   aggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑CTCCCCGCGCGGACAATCa                                                                       >  2:540115/1‑93 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAACAGGGCGTACCAATGCCGATATA‑G‑T‑CACACCGCGCGGACAATCAc                                                                      >  1:2046515/1‑93 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCAc                                                                      <  2:1663164/93‑1 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCAc                                                                      <  1:2214150/93‑1 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCAc                                                                      <  2:826316/93‑1 (MQ=255)

                    ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCAc                                                                      <  2:524563/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACt                                                                     <  2:2692829/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACt                                                                     <  2:2192249/93‑1 (MQ=255)

                     gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACt                                                                     >  1:2347934/1‑93 (MQ=255)

                      aaaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACtt                                                                    <  2:1292087/93‑1 (MQ=255)

                       aaTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTc                                                                   <  1:2661017/93‑1 (MQ=255)

                        aTGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCg                                                                  >  2:2501643/1‑93 (MQ=255)

                         tGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGt                                                                 <  1:4021629/93‑1 (MQ=255)

                          gCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCCCGGACAATCACTTCGTc                                                                >  1:2025508/1‑93 (MQ=255)

                            gantgantTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCaa                                                              <  1:1459683/93‑1 (MQ=255)

                             aGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAAt                                                             <  2:2533582/93‑1 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGAc                                                          <  1:3541283/93‑1 (MQ=255)

                                gATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGAc                                                          <  1:2655079/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         <  1:2540568/93‑1 (MQ=255)

                                 aTTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACt                                                         <  2:1804937/93‑1 (MQ=255)

                                      cccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCaa                                                    <  1:2544472/92‑1 (MQ=255)

                                      accgtcaccgACACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCaa                                                    <  2:2986661/93‑1 (MQ=255)

                                       ccgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAAc                                                   >  1:354144/1‑93 (MQ=255)

                                        cgtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACg                                                  <  1:893348/93‑1 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 >  2:3311129/1‑93 (MQ=255)

                                         gtcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACgg                                                 <  2:2511881/93‑1 (MQ=255)

                                          tcaccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACggg                                                >  2:4330510/1‑93 (MQ=255)

                                           caccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGa                                               <  1:3891217/93‑1 (MQ=255)

                                           caccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TAC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGa                                               >  1:4252023/1‑93 (MQ=255)

                                            cccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAg                                              <  1:1406623/92‑1 (MQ=255)

                                             ccgGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgt                                             <  1:4459647/93‑1 (MQ=255)

                                                gCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                          <  2:806518/93‑1 (MQ=255)

                                                 cacaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgtt                                         <  1:2122889/93‑1 (MQ=255)

                                                 cacaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgtt                                         <  2:3335583/93‑1 (MQ=255)

                                                  acaACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGttgttg                                        >  2:1656091/1‑93 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑GCC‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  1:2878176/93‑1 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  2:310372/93‑1 (MQ=255)

                                                   caACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTtgttgt                                       <  1:1496996/93‑1 (MQ=255)

                                                    aaCCGCGCCCCCAAAAAGAGGGCGACGAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCTCTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      <  2:3790402/93‑1 (MQ=255)

                                                    aaCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGG‑C‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTa                                      <  2:4675942/93‑1 (MQ=255)

                                                     aCCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAg                                     <  2:4320735/93‑1 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  <  1:3706586/93‑1 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  >  1:636398/1‑93 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  >  1:1128528/1‑93 (MQ=255)

                                                        gATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtg                                  <  2:1678367/93‑1 (MQ=255)

                                                         gtCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgt                                 <  2:4116716/92‑1 (MQ=255)

                                                          tCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  2:4074537/93‑1 (MQ=255)

                                                          accGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtg                                <  1:2357519/92‑1 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                               >  2:3790561/1‑93 (MQ=255)

                                                           ccGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                               <  2:1516282/93‑1 (MQ=255)

                                                             gCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgc                             >  1:2225589/1‑93 (MQ=255)

                                                             gCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgc                             <  1:2844112/93‑1 (MQ=255)

                                                              ccAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcg                            >  1:4251238/1‑93 (MQ=255)

                                                               cAAAAAGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGgcgc                           >  1:3952984/1‑93 (MQ=255)

                                                                aaaaaGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGca                          >  2:2167409/1‑93 (MQ=255)

                                                                 aaaaGAGGCGTACCAATGTCAATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGcaa                         >  1:2283779/1‑93 (MQ=255)

                                                                   aaGAGGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGAGcaaca                       >  2:2009764/1‑93 (MQ=255)

                                                                     gagGCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCaacaac                     >  2:897342/1‑93 (MQ=255)

                                                                        gcatACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  1:1040325/90‑1 (MQ=255)

                                                                        gCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  2:701904/93‑1 (MQ=255)

                                                                        gCGTACCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATg                  <  1:940886/93‑1 (MQ=255)

                                                                            aCCAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCa              <  1:3692496/93‑1 (MQ=255)

                                                                             ccAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCGATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             >  1:2505490/1‑93 (MQ=255)

                                                                             ccAAGGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCag             <  1:2505771/93‑1 (MQ=255)

                                                                              cAATGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGGTGTAGCTGTGGGCGCAACAACATGCCCaga            >  1:504856/1‑93 (MQ=255)

                                                                               aaTGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCagag           >  2:3483216/1‑93 (MQ=255)

                                                                               aaTGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCagag           >  1:3080601/1‑93 (MQ=255)

                                                                                aTGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGc          >  2:296371/1‑93 (MQ=255)

                                                                                 tGTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCt         >  1:2841968/1‑93 (MQ=255)

                                                                                  gTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCtt        <  2:2803755/93‑1 (MQ=255)

                                                                                  gTCGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCtt        >  1:597452/1‑93 (MQ=255)

                                                                                    cGATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTa      <  2:14242/93‑1 (MQ=255)

                                                                                     gATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAc     >  2:1995397/1‑93 (MQ=255)

                                                                                     gATATA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAc     >  1:1886784/1‑93 (MQ=255)

                                                                                      atatA‑G‑TCC‑CCCCGCGCGGACCATCACTTAGTCAATCACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCCACAACATGACCAGATCTTTAcc    >  2:917797/1‑93 (MQ=255)

                                                                                      atatA‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTAcc    >  2:1195033/1‑93 (MQ=255)

                                                                                       tata‑G‑TCC‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCa   >  2:4517181/1‑93 (MQ=255)

GCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATA‑G‑T‑C‑CCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAA  >  NC_000913/2563488‑2563667

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