BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,896,196	T→C	C182C (TGT→TGC)	yieJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,896,196	0	T	C	86.7%	687.7 / 89.7	256	C182C (TGT→TGC)	yieJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (18/16); new base C (114/107); total (133/123)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 262 total aligned reads displayed.

GAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCT  >  NC_000913/3896104‑3896287

  agaagaTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGcct                                                                                           >  2:1915983/1‑93 (MQ=255)

     agaTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCa                                                                                        >  1:6520418/1‑93 (MQ=255)

      gaTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCAc                                                                                       >  1:2264210/1‑93 (MQ=255)

      gaTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCAc                                                                                       >  1:754783/1‑93 (MQ=255)

       aTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACt                                                                                      <  1:5006494/93‑1 (MQ=255)

        ttGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTg                                                                                     <  1:5070427/93‑1 (MQ=255)

           tgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCgg                                                                                  <  2:2604107/93‑1 (MQ=255)

            gagaATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGc                                                                                 >  2:2638297/1‑93 (MQ=255)

             agaATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCa                                                                                >  1:4519865/1‑93 (MQ=255)

               aaTTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAg                                                                              <  2:1415890/93‑1 (MQ=255)

                aTTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGc                                                                             >  2:3856609/1‑93 (MQ=255)

                  ttCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTg                                                                           <  2:5099101/93‑1 (MQ=255)

                   tCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTga                                                                          <  1:42855/93‑1 (MQ=255)

                   tCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTga                                                                          >  2:748661/1‑93 (MQ=255)

                   tCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTga                                                                          <  2:6348893/93‑1 (MQ=255)

                     ggTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTgaga                                                                        >  1:3090228/1‑93 (MQ=255)

                      gTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGAc                                                                       >  2:3624135/1‑93 (MQ=255)

                      gTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGAc                                                                       <  2:3656618/93‑1 (MQ=255)

                      gTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGAc                                                                       <  2:5231445/93‑1 (MQ=255)

                       tATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACt                                                                      >  2:5461489/1‑93 (MQ=255)

                       tATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACt                                                                      >  2:541893/1‑93 (MQ=255)

                        aTTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACtg                                                                     <  1:214583/93‑1 (MQ=255)

                        aTTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACtg                                                                     <  1:5289025/93‑1 (MQ=255)

                        aTTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACtg                                                                     <  1:4742746/93‑1 (MQ=255)

                          tAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACtgtg                                                                   <  1:4428889/93‑1 (MQ=255)

                          tAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACtgtg                                                                   <  2:4719705/93‑1 (MQ=255)

                           aGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTgg                                                                  <  1:4551704/93‑1 (MQ=255)

                            gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTggg                                                                 <  1:5389544/93‑1 (MQ=255)

                              aaTTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGt                                                               <  2:575208/93‑1 (MQ=255)

                              aaTTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGt                                                               <  2:2331026/93‑1 (MQ=255)

                                ttCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGa                                                             >  1:4011801/1‑93 (MQ=255)

                                ttCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGa                                                             >  2:3302755/1‑93 (MQ=255)

                                ttCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGa                                                             >  1:382654/1‑93 (MQ=255)

                                    gATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGAtttt                                                         >  1:1132617/1‑93 (MQ=255)

                                    gATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGAtttt                                                         >  2:2816853/1‑93 (MQ=255)

                                      tCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTc                                                       >  2:20711/1‑93 (MQ=255)

                                       cTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCg                                                      <  1:2546086/93‑1 (MQ=255)

                                       cTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCg                                                      >  2:4481234/1‑93 (MQ=255)

                                        tAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGt                                                     >  1:1561898/1‑93 (MQ=255)

                                          gCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAg                                                   >  2:2483784/1‑93 (MQ=255)

                                           cTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGt                                                  <  2:2926132/93‑1 (MQ=255)

                                            tAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGtt                                                 >  2:2538037/1‑93 (MQ=255)

                                            tAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGtt                                                 >  1:874214/1‑93 (MQ=255)

                                             aaaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTa                                                >  2:325545/1‑93 (MQ=255)

                                              aaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTAt                                               >  2:1832070/1‑93 (MQ=255)

                                              aaTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTAt                                               >  2:2615513/1‑93 (MQ=255)

                                                tGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTAttt                                             <  2:6545835/93‑1 (MQ=255)

                                                tGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTAttt                                             <  2:6052976/93‑1 (MQ=255)

                                                tGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTAttt                                             <  2:6292280/93‑1 (MQ=255)

                                                  ccTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTaa                                           <  2:6315641/93‑1 (MQ=255)

                                                  ccTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTaa                                           >  2:3063271/1‑93 (MQ=255)

                                                  ccTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTaa                                           >  2:2844606/1‑93 (MQ=255)

                                                    tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaat                                         <  2:6625968/93‑1 (MQ=255)

                                                     gCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaata                                        >  1:3474176/1‑93 (MQ=255)

                                                      cAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaataa                                       <  2:5047951/93‑1 (MQ=255)

                                                      cAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaataa                                       <  1:4823182/93‑1 (MQ=255)

                                                       aaaaGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaataat                                      >  1:3740053/1‑93 (MQ=255)

                                                       aaaaGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAaataat                                      >  2:5287919/1‑93 (MQ=255)

                                                        aaaGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATg                                     >  1:5623997/1‑93 (MQ=255)

                                                        aaaGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATg                                     >  1:5316295/1‑93 (MQ=255)

                                                         aaGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATga                                    <  2:3217641/93‑1 (MQ=255)

                                                          aGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgag                                   <  2:3805596/93‑1 (MQ=255)

                                                          aGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgag                                   >  1:1076239/1‑93 (MQ=255)

                                                          aGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgag                                   >  2:2285017/1‑93 (MQ=255)

                                                           gggTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgaga                                  >  1:2994071/1‑93 (MQ=255)

                                                           gggTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgaga                                  <  1:1120008/93‑1 (MQ=255)

                                                           gggTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgaga                                  >  2:4861353/1‑93 (MQ=255)

                                                           gggTGGTCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATgaga                                  <  2:3462996/93‑1 (MQ=255)

                                                            ggtggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGaa                                 <  2:2963141/93‑1 (MQ=255)

                                                             gtggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAAc                                <  2:318478/93‑1 (MQ=255)

                                                             gtggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAAc                                <  1:4872233/93‑1 (MQ=255)

                                                              tggtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACa                               <  1:4235979/93‑1 (MQ=255)

                                                              tggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACa                               >  1:466583/1‑93 (MQ=255)

                                                               ggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAg                              >  2:1246597/1‑93 (MQ=255)

                                                               ggtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAg                              <  1:312629/93‑1 (MQ=255)

                                                                gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAgg                             >  1:1002259/1‑93 (MQ=255)

                                                                gtCATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAgg                             >  1:1410306/1‑93 (MQ=255)

                                                                  cATTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAggac                           >  1:6152192/1‑91 (MQ=255)

                                                                   aTTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCg                          <  1:3996372/93‑1 (MQ=255)

                                                                   aTTGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGGGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCg                          >  1:1679150/1‑93 (MQ=255)

                                                                    ttGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCgg                         >  1:5126569/1‑93 (MQ=255)

                                                                     tGTCAGGGTTATCTCTTCCGCTGCCTCCGCTGCGGCAAGCTTAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGa                        <  2:1592246/93‑1 (MQ=255)

                                                                     tGTCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGa                        <  1:6385664/93‑1 (MQ=255)

                                                                      gTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAg                       >  1:3220873/1‑93 (MQ=255)

                                                                       tCAGGGTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGc                      >  1:283093/1‑93 (MQ=255)

                                                                        cAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCg                     >  1:355475/1‑93 (MQ=255)

                                                                         aGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGt                    >  1:4190421/1‑93 (MQ=255)

                                                                            gTTATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAAt                 >  1:5551712/1‑93 (MQ=255)

                                                                              tATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTc               >  1:6086687/1‑93 (MQ=255)

                                                                              tATCTCTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTc               >  2:1010526/1‑93 (MQ=255)

                                                                                tctcTTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTcac             >  2:4339058/1‑93 (MQ=255)

                                                                                  tctTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTcacac           <  1:2901291/93‑1 (MQ=255)

                                                                                   ctTCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTcacaca          <  2:4110672/93‑1 (MQ=255)

                                                                                     tCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc        <  2:5601094/93‑1 (MQ=255)

                                                                                     tCCGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGCACGGAGCGTAATTcacacacc        >  2:4619600/1‑91 (MQ=255)

                                                                                      ccGCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCTGAGCGTAATTCACACATcc       >  2:6113060/1‑93 (MQ=255)

                                                                                        gCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCgg     >  2:2394947/1‑93 (MQ=255)

                                                                                        gCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCgg     >  2:400848/1‑93 (MQ=255)

                                                                                        gCTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCgg     <  1:3350063/93‑1 (MQ=255)

                                                                                         cTGCCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGc    >  1:596357/1‑93 (MQ=255)

GAAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCT  >  NC_000913/3896104‑3896287

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