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	99.6%	821.1 / ‑3.9	227	coding (236/282 nt)	b4568	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (126/100); total (126/101)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 244 total aligned reads displayed.

GTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005509‑4005687

gTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑c                                                                                        <  2:2742083/93‑1 (MQ=255)

   aTGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑cctg                                                                                     >  2:3739907/1‑90 (MQ=255)

      gTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGCGGCCTGTAt                                                                                   <  2:5230122/93‑1 (MQ=255)

       tttGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGCGCGGTGGCCTGTAtc                                                                                  <  1:2997060/93‑1 (MQ=255)

          gggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGGGGTGGCCTGTATctct                                                                               >  1:4392296/1‑93 (MQ=255)

           ggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                                              <  2:794724/93‑1 (MQ=255)

           ggCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctg                                                                              >  1:5839262/1‑93 (MQ=255)

            gCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgc                                                                             <  1:5104097/93‑1 (MQ=255)

             cTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgct                                                                            >  2:4879192/1‑93 (MQ=255)

              tGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTctgctg                                                                           <  2:4329743/93‑1 (MQ=255)

               gggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTgg                                                                          >  1:4248257/1‑93 (MQ=255)

                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                         <  1:6075749/93‑1 (MQ=255)

                ggATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGa                                                                         >  2:5126022/1‑93 (MQ=255)

                 gATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGAt                                                                        >  1:1758267/1‑93 (MQ=255)

                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                       <  1:376044/93‑1 (MQ=255)

                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                       <  1:3079817/93‑1 (MQ=255)

                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                       <  2:2159850/93‑1 (MQ=255)

                  aTTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATg                                                                       >  1:829387/1‑93 (MQ=255)

                   ttGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATgg                                                                      >  2:438275/1‑93 (MQ=255)

                     gcgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGt                                                                    <  2:1442505/93‑1 (MQ=255)

                      cgcTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGtt                                                                   >  2:2862628/1‑93 (MQ=255)

                       gctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGGGGTGGCCTGTATCTCTGCTGGATGGGTTa                                                                  >  2:693515/1‑93 (MQ=255)

                        ctgctTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTAc                                                                 >  1:5322887/1‑93 (MQ=255)

                            ttGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGa                                                             <  2:894798/93‑1 (MQ=255)

                             tGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGAt                                                            >  1:3637726/1‑93 (MQ=255)

                              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGCTGGGTTACAAGATg                                                           >  1:2165659/1‑93 (MQ=255)

                              ggCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATg                                                           >  1:2601003/1‑93 (MQ=255)

                               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                                          <  1:375124/93‑1 (MQ=255)

                               gCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGc                                                          >  2:2218192/1‑93 (MQ=255)

                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGGGGTGGCCTGTATATCTGCTGGATGGGTTACCAGATGCt                                                         >  2:1136406/1‑93 (MQ=255)

                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         >  1:5941374/1‑93 (MQ=255)

                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         >  2:1194259/1‑93 (MQ=255)

                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCt                                                         >  1:3691002/1‑93 (MQ=255)

                                ccTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGGTACCAGATTCt                                                         >  1:2086386/1‑93 (MQ=255)

                                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        >  2:3878374/1‑93 (MQ=255)

                                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTa                                                        >  2:5106187/1‑93 (MQ=255)

                                 cTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCAGCTGGATGGGTTACCAGATGCTa                                                        >  1:3500820/1‑93 (MQ=255)

                                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  2:5527550/1‑93 (MQ=255)

                                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  2:3449117/1‑93 (MQ=255)

                                  tGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTAc                                                       >  2:4990037/1‑93 (MQ=255)

                                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                                      >  2:1571778/1‑93 (MQ=255)

                                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                                      <  1:2817053/93‑1 (MQ=255)

                                   gCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACg                                                      <  2:3007826/93‑1 (MQ=255)

                                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                     <  1:379429/93‑1 (MQ=255)

                                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                     <  2:4288825/93‑1 (MQ=255)

                                    cATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgt                                                     >  1:1345075/1‑93 (MQ=255)

                                     aTTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtg                                                    >  1:931982/1‑93 (MQ=255)

                                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                   >  1:991840/1‑93 (MQ=255)

                                      tttGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtgg                                                   >  2:1289899/1‑93 (MQ=255)

                                         gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACgtggtga                                                >  1:6127543/1‑92 (MQ=255)

                                         gATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGc                                                >  2:3965784/1‑93 (MQ=255)

                                           ttATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCAc                                              >  2:4618061/1‑93 (MQ=255)

                                            tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             <  2:4253889/93‑1 (MQ=255)

                                            tATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             <  2:2662522/93‑1 (MQ=255)

                                            tATCGAAAAAATGGCCTGGCTACATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTCCTGGATGGGTTACCAGATGCTACGTGGTGCACt                                             <  1:5167213/93‑1 (MQ=255)

                                             aTCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTg                                            <  1:4169752/93‑1 (MQ=255)

                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                           >  1:2690293/1‑93 (MQ=255)

                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGa                                           <  1:1001754/93‑1 (MQ=255)

                                              tCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGCGGCCTGTATCTCTGCTGGATGGGACACCAGATGCTACGTGGCGCACTGa                                           >  1:634097/1‑93 (MQ=255)

                                                gAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaa                                         >  1:1475015/1‑93 (MQ=255)

                                                  aaaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaa                                       <  1:856294/93‑1 (MQ=255)

                                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                      <  2:590716/93‑1 (MQ=255)

                                                   aaaaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGaaaaaa                                      <  1:4019550/93‑1 (MQ=255)

                                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                    <  2:2034476/93‑1 (MQ=255)

                                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                    <  1:3182281/93‑1 (MQ=255)

                                                     aaTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAaga                                    <  2:3833571/93‑1 (MQ=255)

                                                      aTGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAagag                                   <  1:5604024/93‑1 (MQ=255)

                                                         gCCTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCg                                <  1:3397112/93‑1 (MQ=255)

                                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                              <  2:1634052/93‑1 (MQ=255)

                                                           cTGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGt                              >  2:4562110/1‑93 (MQ=255)

                                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCgggt                             >  1:5415244/1‑91 (MQ=255)

                                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             <  2:3249152/93‑1 (MQ=255)

                                                            tGGCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGtt                             >  2:5369969/1‑93 (MQ=255)

                                                             ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt                            <  1:1196861/93‑1 (MQ=255)

                                                             ggCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGttt                            >  1:5065532/1‑93 (MQ=255)

                                                              gCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           >  1:1794470/1‑93 (MQ=255)

                                                              gCTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTc                           <  2:2501090/93‑1 (MQ=255)

                                                               cTGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCt                          <  1:3538590/93‑1 (MQ=255)

                                                                tGCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTg                         >  1:106291/1‑93 (MQ=255)

                                                                 gCATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGc                        >  2:416641/1‑93 (MQ=255)

                                                                  cATACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCa                       >  2:3454761/1‑93 (MQ=255)

                                                                   aTACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAc                      >  2:3118509/1‑93 (MQ=255)

                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                     <  2:2727876/93‑1 (MQ=255)

                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                     >  1:4491564/1‑93 (MQ=255)

                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCAcc                     >  1:4776762/1‑93 (MQ=255)

                                                                    tACGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGGTTCTGCAcc                     >  1:193729/1‑93 (MQ=255)

                                                                      cGCTGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTg                   >  1:1777963/1‑93 (MQ=255)

                                                                         tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc                >  2:1718479/1‑93 (MQ=255)

                                                                         tGATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTgcgc                <  1:1232432/93‑1 (MQ=255)

                                                                          gATTATGGTGGGCGGTGGCCTGTATCTCTGCTGGGTGGGTTACCAGATGGTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTCCGcc               >  2:4895506/1‑93 (MQ=255)

                                                                            ttATGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCcac             <  2:5047659/93‑1 (MQ=255)

                                                                              aTGGTGGGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAg           <  2:443573/93‑1 (MQ=255)

                                                                               tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg          >  2:1956131/1‑93 (MQ=255)

                                                                               tggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAgg          >  1:5022680/1‑93 (MQ=255)

                                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt         <  2:1819882/93‑1 (MQ=255)

                                                                                ggtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGCGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGt         >  2:2689064/1‑93 (MQ=255)

                                                                                 gtggGCGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTc        <  1:5246909/93‑1 (MQ=255)

                                                                                      cGGTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCt   <  1:1974706/93‑1 (MQ=255)

                                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  <  1:4519921/93‑1 (MQ=255)

                                                                                       ggTGGCCTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGCGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTg  >  1:1511635/1‑93 (MQ=255)

GTAATGGTTTGGGCTGGGATTGCGCTGCTTGGCCTGCATTTGATTATCGAAAAAATGGCCTGGCTGCATACGCTGATTATGGTGGGCGGTGG‑CTGTATCTCTGCTGGATGGGTTACCAGATGCTACGTGGTGCACTGAAAAAAGAGGCGGTTTCTGCACCTGCGCCACAGGTCGAGCTG  >  NC_000913/4005509‑4005687

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