BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,005,601	+C	coding (237/282 nt)	b4568 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,005,600	1	.	C	98.9%	866.0 / 0.6	267	coding (236/282 nt)	b4568	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (159/105); total (160/107)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCT  >  NC_000913/4005508‑4005686

 gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑c                                                                                       >  1:6657839/1‑93 (MQ=255)

 gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑c                                                                                       >  2:6051523/1‑93 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑cc                                                                                      <  2:6291484/93‑2 (MQ=255)

  tAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑cc                                                                                      >  1:5561854/1‑92 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑cct                                                                                     >  2:1562753/1‑91 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑cct                                                                                     >  1:5670371/1‑91 (MQ=255)

   aaTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑cct                                                                                     >  2:5086724/1‑91 (MQ=255)

       gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTAt                                                                                  >  1:4254551/1‑93 (MQ=255)

       gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTAt                                                                                  >  2:4123908/1‑93 (MQ=255)

         ttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTAtct                                                                                <  1:5210057/93‑1 (MQ=255)

          tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTAtctc                                                                               <  2:3536324/93‑1 (MQ=255)

          tGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTAtctc                                                                               <  1:1881745/93‑1 (MQ=255)

              cTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTctgct                                                                           >  2:3767505/1‑93 (MQ=255)

               tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTctgctg                                                                          >  2:2359520/1‑93 (MQ=255)

                gggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTgg                                                                         >  2:669709/1‑93 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTggg                                                                        >  2:1766700/1‑92 (MQ=255)

                 ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGa                                                                        >  1:3814320/1‑93 (MQ=255)

                  gATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGAt                                                                       >  1:5829507/1‑93 (MQ=255)

                   aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGAGGT‑G‑GCCTGTATCTCTGCTGGATg                                                                      >  1:1023851/1‑93 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATgg                                                                     <  1:2223460/93‑1 (MQ=255)

                    ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATgg                                                                     <  2:5984791/93‑1 (MQ=255)

                      gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGt                                                                   >  1:3737431/1‑93 (MQ=255)

                       cgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGtt                                                                  <  2:6026540/93‑1 (MQ=255)

                       cgcTGCTTGGCATGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGtt                                                                  <  1:541797/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTAcc                                                               <  1:2147601/93‑1 (MQ=255)

                          tgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTAcc                                                               <  2:3394531/93‑1 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              >  2:5232792/1‑93 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              >  2:2463021/1‑93 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              <  1:2583431/93‑1 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              >  1:6373030/1‑93 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              >  2:3151255/1‑93 (MQ=255)

                           gctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCa                                                              >  1:2978079/1‑93 (MQ=255)

                             ttGTCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGa                                                            <  2:6186258/93‑1 (MQ=255)

                             ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGa                                                            <  1:2821856/93‑1 (MQ=255)

                              tGGCCTGCATTTGATTATCGAAAGAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                           <  1:6725286/93‑1 (MQ=255)

                              tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                           >  1:3501523/1‑93 (MQ=255)

                                gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                                         >  2:882833/1‑93 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                        >  1:305138/1‑93 (MQ=255)

                                 ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                        <  2:5664256/93‑1 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                       >  1:2489868/1‑93 (MQ=255)

                                  cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                       >  2:39964/1‑93 (MQ=255)

                                   tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                      >  1:5901929/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTCG‑C‑CTGTATCTCTGCTGGATGGGATACCAGATGCTACgt                                                    >  1:6450779/1‑93 (MQ=255)

                                     cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                    >  1:3604643/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  2:2788470/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   <  2:3278758/93‑1 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  2:4447864/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  1:3290363/1‑93 (MQ=255)

                                      aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGACTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                   >  2:4831347/1‑93 (MQ=255)

                                         tGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtg                                                <  1:1478756/93‑1 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                               >  1:6230823/1‑93 (MQ=255)

                                          gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                               >  1:1774094/1‑93 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                             <  1:6503639/93‑1 (MQ=255)

                                            ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                             >  1:4004795/1‑93 (MQ=255)

                                              aTCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                                           <  1:709612/93‑1 (MQ=255)

                                                cGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaa                                         >  1:4573524/1‑93 (MQ=255)

                                                  aaaaaaTGTCCTGGCTGCATACGCTGATTATGGTGGGCGCT‑G‑GCCTGTATATCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaa                                       <  2:3024024/93‑1 (MQ=255)

                                                   aaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaa                                      >  2:2867737/1‑93 (MQ=255)

                                                    aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                     >  1:611449/1‑93 (MQ=255)

                                                    aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                     >  2:2206936/1‑93 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                   >  2:1965478/1‑93 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                   >  2:1457715/1‑93 (MQ=255)

                                                      aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                   >  2:5869492/1‑93 (MQ=255)

                                                       aTGGCGTGGCTGTATACGCTCATAATGCTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTCACCAGATGCTACGTGGTGCACTGAAAAAagag                                  <  1:4664235/93‑1 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                  >  2:4182910/1‑93 (MQ=255)

                                                       aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                  >  1:3365214/1‑93 (MQ=255)

                                                        tGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAgg                                 <  1:936202/93‑1 (MQ=255)

                                                         ggCGTGGCTGCATACGCTGATTATGGTGGGCGGT‑GTC‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGc                                <  1:4705942/93‑1 (MQ=255)

                                                          gCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTCCCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg                               <  1:4329285/93‑1 (MQ=255)

                                                            cTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                             >  2:128740/1‑93 (MQ=255)

                                                            cTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                             <  2:78359/93‑1 (MQ=255)

                                                             tGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                            >  1:2742253/1‑93 (MQ=255)

                                                               gCTGCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                          >  2:4669476/1‑93 (MQ=255)

                                                                  gCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                       >  2:1631224/1‑93 (MQ=255)

                                                                  gCATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                       >  2:4498996/1‑93 (MQ=255)

                                                                   cATACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCa                      <  2:4050077/93‑1 (MQ=255)

                                                                   cATACGCTGATTAGGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCa                      <  2:6463662/93‑1 (MQ=255)

                                                                    aTACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                     >  2:6502511/1‑93 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGTGG‑G‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                    <  1:1969936/93‑1 (MQ=255)

                                                                     tACGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                    <  1:434723/93‑1 (MQ=255)

                                                                       cGCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTg                  >  2:388781/1‑93 (MQ=255)

                                                                        gCTGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgc                 >  2:6239769/1‑93 (MQ=255)

                                                                          tGATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc               >  1:1175524/1‑93 (MQ=255)

                                                                           gATTATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGcc              >  1:1300112/1‑93 (MQ=255)

                                                                             ttATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            <  2:6563641/93‑1 (MQ=255)

                                                                             ttATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac            <  1:5318336/93‑1 (MQ=255)

                                                                              tATGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcaca           <  2:146929/93‑1 (MQ=255)

                                                                               aTGGTGGGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAg          >  1:1732979/1‑93 (MQ=255)

                                                                                tggtggGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  1:4186107/1‑93 (MQ=255)

                                                                                tggtggGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  2:4289290/1‑93 (MQ=255)

                                                                                tggtggGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg         >  1:5024717/1‑93 (MQ=255)

                                                                                 ggtggGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt        <  2:5498175/93‑1 (MQ=255)

                                                                                   tggGCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCg      >  2:3517145/1‑93 (MQ=255)

                                                                                    gggCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa     >  1:5979294/1‑93 (MQ=255)

                                                                                    gggCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGa     <  1:2316249/93‑1 (MQ=255)

                                                                                     ggCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAg    >  2:6445047/1‑93 (MQ=255)

                                                                                     ggCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAg    <  2:4181191/93‑1 (MQ=255)

                                                                                      gCGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGc   <  1:2098401/93‑1 (MQ=255)

                                                                                       cGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCt  <  1:6285347/93‑1 (MQ=255)

                                                                                       cGGT‑G‑GCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCt  <  1:3167576/93‑1 (MQ=255)

CGTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGT‑G‑G‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCT  >  NC_000913/4005508‑4005686

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