BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,479,655	T→C	C241C (TGT→TGC)	ynbD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,479,655	0	T	C	99.2%	390.5 / ‑4.7	122	C241C (TGT→TGC)	ynbD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (53/68); total (54/68)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

Only 100 of 124 total aligned reads displayed.

AGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTACC  >  NC_000913/1479565‑1479744

aGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTg                                                                                         >  1:3708395/1‑93 (MQ=255)

aGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTg                                                                                         <  2:4733379/93‑1 (MQ=255)

  cACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGc                                                                                       <  1:433061/93‑1 (MQ=255)

    cTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCa                                                                                     >  1:1948564/1‑93 (MQ=255)

    cTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCa                                                                                     >  1:2685459/1‑93 (MQ=255)

       cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGta                                                                                  <  1:544326/93‑1 (MQ=255)

       cccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGta                                                                                  <  2:3267912/93‑1 (MQ=255)

        ccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGtat                                                                                 >  1:3661785/1‑93 (MQ=255)

        ccTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGtat                                                                                 >  1:2847646/1‑93 (MQ=255)

          tATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGtatta                                                                               >  1:2846054/1‑93 (MQ=255)

           aTGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGtattat                                                                              >  2:328101/1‑93 (MQ=255)

              tCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGc                                                                           <  1:2956393/93‑1 (MQ=255)

                  aGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACt                                                                       <  1:1031691/93‑1 (MQ=255)

                   ggCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTc                                                                      <  1:1746119/93‑1 (MQ=255)

                    gcgcGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCa                                                                     <  2:266213/93‑1 (MQ=255)

                    gcgcGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCa                                                                     >  1:3996070/1‑93 (MQ=255)

                    gcgcGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCa                                                                     <  1:1573266/93‑1 (MQ=255)

                       cgGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATaa                                                                  <  2:197428/93‑1 (MQ=255)

                       cgGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATaa                                                                  >  1:2956625/1‑93 (MQ=255)

                         ggCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAAtt                                                                >  1:4304586/1‑93 (MQ=255)

                         ggCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAAtt                                                                <  2:4603636/93‑1 (MQ=255)

                         ggCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAAtt                                                                >  2:769697/1‑93 (MQ=255)

                            tCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGc                                                             >  2:2019599/1‑93 (MQ=255)

                              gTGCATTGTATTGATTGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCg                                                           <  1:965538/93‑1 (MQ=255)

                              gTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCg                                                           >  1:3876763/1‑93 (MQ=255)

                              gTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCg                                                           >  1:1031797/1‑93 (MQ=255)

                               tGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGt                                                          >  1:3038187/1‑93 (MQ=255)

                                 cATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTgg                                                        >  1:2668088/1‑93 (MQ=255)

                                 cATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTgg                                                        <  2:1176651/93‑1 (MQ=255)

                                  aTTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTggg                                                       >  2:1859831/1‑93 (MQ=255)

                                  aTTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTggg                                                       >  2:1859985/1‑93 (MQ=255)

                                   ttGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGt                                                      >  1:1681171/1‑93 (MQ=255)

                                    tGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTa                                                     <  2:3831622/93‑1 (MQ=255)

                                     gTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTAc                                                    <  1:2639614/93‑1 (MQ=255)

                                      tATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGTCCGTGGGTACg                                                   >  1:1182126/1‑93 (MQ=255)

                                      tATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACg                                                   >  1:4550271/1‑93 (MQ=255)

                                      tATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACg                                                   >  1:1450785/1‑93 (MQ=255)

                                       aTTGAGGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACgg                                                  <  2:805854/93‑1 (MQ=255)

                                          gATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTgg                                               >  1:2704175/1‑93 (MQ=255)

                                             ggAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCt                                            <  2:3971073/93‑1 (MQ=255)

                                             ggAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCt                                            >  2:122128/1‑93 (MQ=255)

                                             ggAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCt                                            <  1:1461375/93‑1 (MQ=255)

                                              gAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCtt                                           <  2:3355635/93‑1 (MQ=255)

                                              gAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCtt                                           >  1:4588798/1‑93 (MQ=255)

                                               aGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTg                                          <  2:16099/93‑1 (MQ=255)

                                               aGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTg                                          <  1:4301895/93‑1 (MQ=255)

                                                 cTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGgc                                        <  1:910053/93‑1 (MQ=255)

                                                   aaTGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGgcgc                                      >  2:1540223/1‑93 (MQ=255)

                                                     tgatgatgaTTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGa                                    <  1:431641/93‑1 (MQ=255)

                                                      gatgatgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGAt                                   <  2:4235111/93‑1 (MQ=255)

                                                      gatgatgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGAt                                   <  1:522227/93‑1 (MQ=255)

                                                        tgatgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaa                                 >  1:1730634/1‑93 (MQ=255)

                                                         gatgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaac                                <  2:2641766/93‑1 (MQ=255)

                                                         gatgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaac                                <  2:2229070/93‑1 (MQ=255)

                                                          atgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaaca                               <  1:3236648/93‑1 (MQ=255)

                                                           tgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaa                              <  1:2176989/93‑1 (MQ=255)

                                                           tgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaa                              >  1:4415893/1‑93 (MQ=255)

                                                           tgatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaa                              >  1:1783515/1‑93 (MQ=255)

                                                            gatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaac                             <  2:4658329/93‑1 (MQ=255)

                                                            gatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaac                             >  1:1101307/1‑93 (MQ=255)

                                                            gatTCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATaacaac                             >  1:3907764/1‑93 (MQ=255)

                                                              ttCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAg                           >  1:607081/1‑93 (MQ=255)

                                                              ttCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAg                           >  2:2972721/1‑93 (MQ=255)

                                                              ttCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAg                           <  1:3683324/93‑1 (MQ=255)

                                                              ttCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAg                           >  2:1447062/1‑93 (MQ=255)

                                                               tCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAgg                          <  1:3130864/93‑1 (MQ=255)

                                                               tCAGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAgg                          >  1:4741484/1‑93 (MQ=255)

                                                                 aGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGa                        <  2:3774659/93‑1 (MQ=255)

                                                                 aGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGa                        <  2:3766609/93‑1 (MQ=255)

                                                                 aGTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGa                        >  1:1236031/1‑93 (MQ=255)

                                                                  gTTATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTTGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaa                       <  1:2349735/93‑1 (MQ=255)

                                                                   ttATGGTGGTCAGTCTGGTTATGCTTGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGAGCGATAACAACAGGGaaa                      <  2:3169837/93‑1 (MQ=255)

                                                                   ttATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaaa                      <  2:692695/93‑1 (MQ=255)

                                                                   ttATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaaa                      <  1:3706499/93‑1 (MQ=255)

                                                                   ttATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGaaa                      <  1:2801684/93‑1 (MQ=255)

                                                                   ttATGGTGGTCAGTCTGGTTATGCTGGCCAGGATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCGTGGCGCGATAACAACAGGGaaa                      >  1:2010491/1‑93 (MQ=255)

                                                                    tATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAg                     <  2:1266700/93‑1 (MQ=255)

                                                                    tATGGTGGTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAg                     <  2:2809106/93‑1 (MQ=255)

                                                                     aTGGTGGTCAGTATGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGa                    <  1:1787405/93‑1 (MQ=255)

                                                                      tggtggTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGAt                   <  2:578349/93‑1 (MQ=255)

                                                                      tggtggTCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGAt                   <  1:2798083/93‑1 (MQ=255)

                                                                       ggtggtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATa                  <  1:4431430/93‑1 (MQ=255)

                                                                       ggtggtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATa                  <  2:4187604/93‑1 (MQ=255)

                                                                        gtggtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAg                 <  1:2118380/93‑1 (MQ=255)

                                                                        gtggtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAg                 <  2:4611605/93‑1 (MQ=255)

                                                                          ggtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTc               <  2:1338447/93‑1 (MQ=255)

                                                                           gtCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCa              >  2:1821831/1‑93 (MQ=255)

                                                                            tCAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAg             <  1:2123493/93‑1 (MQ=255)

                                                                             cAGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAgg            <  2:4445380/93‑1 (MQ=255)

                                                                              aGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAggg           >  1:934046/1‑93 (MQ=255)

                                                                              aGTCTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGGACGGTGGGCTTGGCGCGATAACGACAGGGAAAGATAGGCAgag           >  1:934047/1‑91 (MQ=255)

                                                                              aGTCTGGTTATGCTGGCAAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAggg           <  2:4171923/93‑1 (MQ=255)

                                                                               gTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAgggg          >  1:3863025/1‑93 (MQ=255)

                                                                                 cTGGTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCCTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGaa        >  1:73250/1‑93 (MQ=255)

                                                                                    gTTATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACt     >  1:1966596/1‑93 (MQ=255)

                                                                                     ttATGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTa    <  2:2631131/93‑1 (MQ=255)

                                                                                       aTGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTAcc  >  2:3393181/1‑93 (MQ=255)

                                                                                       aTGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTAcc  <  1:2074793/93‑1 (MQ=255)

                                                                                       aTGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTAcc  >  2:3168281/1‑93 (MQ=255)

                                                                                       aTGCTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTAcc  <  2:2354202/93‑1 (MQ=255)

AGCACTGCCCTATGTCGTAGGCGCGGGCTCGTGCATTGTATTGATGGAGCTAATGATGATGATTCAGTTATGGTGGTCAGTCTGGTTATGTTGGCCAGTATTATCGCTACTCATAATTGGCCGTGGGTACGGTGGGCTTGGCGCGATAACAACAGGGAAAGATAGTCAGGGGAAACTACC  >  NC_000913/1479565‑1479744

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