BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,455,292	T→C	W427R (TGG→CGG)	b3325 →	YheF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,455,292	0	T	C	85.1%	620.1 / 92.3	268	W427R (TGG→CGG)	b3325	YheF
Reads supporting (aligned to +/- strand): ref base T (21/19); new base C (100/128); total (121/147)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 270 total aligned reads displayed.

AAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCC  >  NC_000913/3455202‑3455384

aaaaaGAATGGTGGGATCACCAGCGCGGATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgg                                                                                            >  2:6422025/1‑93 (MQ=255)

aaaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTgg                                                                                            <  2:2254181/93‑1 (MQ=255)

aaaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgg                                                                                            <  1:4888306/93‑1 (MQ=255)

aaaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgg                                                                                            <  2:1259933/93‑1 (MQ=255)

aaaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgg                                                                                            >  2:528247/1‑93 (MQ=255)

 aaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCggg                                                                                           >  2:3779638/1‑93 (MQ=255)

   aaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGa                                                                                         <  1:1458081/93‑1 (MQ=255)

   aaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGa                                                                                         <  2:6002723/93‑1 (MQ=255)

   aaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGa                                                                                         <  1:7376622/93‑1 (MQ=255)

    aGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAg                                                                                        >  1:2216653/1‑93 (MQ=255)

    aGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAg                                                                                        <  2:3075175/93‑1 (MQ=255)

     gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGt                                                                                       >  2:5725982/1‑93 (MQ=255)

     gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGt                                                                                       <  1:2158892/93‑1 (MQ=255)

     gAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGt                                                                                       >  2:3653089/1‑93 (MQ=255)

      aaTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTa                                                                                      <  2:4988070/93‑1 (MQ=255)

       aTGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTAc                                                                                     >  2:3384798/1‑93 (MQ=255)

        tggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACt                                                                                    <  1:3067244/93‑1 (MQ=255)

         ggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTg                                                                                   >  2:1030180/1‑93 (MQ=255)

          gtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGc                                                                                  <  2:2961612/93‑1 (MQ=255)

            gggATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCtt                                                                                <  2:924767/93‑1 (MQ=255)

              gATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAc                                                                              <  2:6170875/93‑1 (MQ=255)

                tCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcg                                                                            >  2:5839207/1‑93 (MQ=255)

                 cACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgc                                                                           <  2:2085397/93‑1 (MQ=255)

                    cAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCt                                                                        <  1:473269/93‑1 (MQ=255)

                     aGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTg                                                                       <  1:6090721/93‑1 (MQ=255)

                      gcgcGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTCCCGCGCTgg                                                                      <  1:2372512/93‑1 (MQ=255)

                           aaTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                                 >  1:1554283/1‑93 (MQ=255)

                           aaTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGt                                                                 <  1:3560078/93‑1 (MQ=255)

                            aTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTa                                                                >  1:1056268/1‑93 (MQ=255)

                             tCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTaa                                                               <  1:3286537/93‑1 (MQ=255)

                             tCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTaa                                                               <  2:4641781/93‑1 (MQ=255)

                             tCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTaa                                                               <  2:5446896/93‑1 (MQ=255)

                             tCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTaa                                                               <  1:6346015/93‑1 (MQ=255)

                              cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                              <  2:2639639/93‑1 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                           <  1:533002/93‑1 (MQ=255)

                                 gcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAAt                                                           <  1:871571/93‑1 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                         <  1:5888758/93‑1 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                         <  2:5668613/93‑1 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                         >  2:2500263/1‑93 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                         >  1:7417471/1‑93 (MQ=255)

                                   ctgGGATATGTTTAGCGCCTACAATAGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaa                                                         >  1:2890205/1‑93 (MQ=255)

                                    tgGGATATGTTTAGCGCCTACAATGGCATGTCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaa                                                        <  2:3240545/93‑1 (MQ=255)

                                    tgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaa                                                        >  1:731267/1‑93 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                       >  1:290644/1‑93 (MQ=255)

                                     gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                       >  1:5715885/1‑93 (MQ=255)

                                      ggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaaa                                                      >  2:1177101/1‑93 (MQ=255)

                                       gATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAAt                                                     >  1:3929180/1‑93 (MQ=255)

                                        atatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                    <  1:1370883/93‑1 (MQ=255)

                                        atatGTTTAGCCCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGGCCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                    <  1:3717005/93‑1 (MQ=255)

                                         tatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGa                                                   <  2:4815516/93‑1 (MQ=255)

                                         tatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGa                                                   <  1:3230718/93‑1 (MQ=255)

                                         tatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGa                                                   >  2:1841356/1‑93 (MQ=255)

                                          atGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGAc                                                  <  1:496269/93‑1 (MQ=255)

                                          atGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGAc                                                  <  1:7000764/93‑1 (MQ=255)

                                          atGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGAc                                                  <  1:5207570/93‑1 (MQ=255)

                                            gTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACAt                                                <  1:4264345/93‑1 (MQ=255)

                                            gTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACAt                                                <  2:3272884/93‑1 (MQ=255)

                                               tAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCt                                             <  2:3465366/93‑1 (MQ=255)

                                                aGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTc                                            <  2:6902688/93‑1 (MQ=255)

                                                 gcgcCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCg                                           >  2:1909803/1‑93 (MQ=255)

                                                 gcgcCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCg                                           >  2:7117342/1‑93 (MQ=255)

                                                  cgcCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGc                                          <  2:5230404/93‑1 (MQ=255)

                                                     cttcAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGTTGGCCAGTAACAATAAAAATGACATCCTCGCCAc                                       <  1:4322683/90‑1 (MQ=255)

                                                      tACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAcc                                      >  2:6309529/1‑93 (MQ=255)

                                                       aCAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAccc                                     >  2:2420716/1‑93 (MQ=255)

                                                          atggcatggcCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCaa                                  >  2:6777971/1‑93 (MQ=255)

                                                          atggcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCaa                                  <  1:1333205/93‑1 (MQ=255)

                                                             gcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCa                               <  1:1762726/93‑1 (MQ=255)

                                                             gcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCa                               >  2:1731057/1‑93 (MQ=255)

                                                              catggcCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                              <  2:800652/93‑1 (MQ=255)

                                                              catggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                              <  2:3894060/93‑1 (MQ=255)

                                                              catggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                              <  2:4477635/93‑1 (MQ=255)

                                                                tggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCg                            >  1:6032096/1‑93 (MQ=255)

                                                                  gcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTa                          >  2:4164869/1‑93 (MQ=255)

                                                                    cGCAGGCTTCTTCAATGGGGACCGGGGAGTACTGCTTACCGAGCTGGCCAGTAACAATAAAAATCACATCCTCGCCGCCCCCAGCATCGTAAc                        >  1:6174796/1‑93 (MQ=255)

                                                                    cGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAAc                        <  1:4551488/93‑1 (MQ=255)

                                                                     gCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                       >  1:3997256/1‑93 (MQ=255)

                                                                     gCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                       >  1:1843618/1‑93 (MQ=255)

                                                                     gCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACg                       >  1:574390/1‑93 (MQ=255)

                                                                       aGGCTTCTACAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCt                     <  2:5628351/93‑1 (MQ=255)

                                                                        ggCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTg                    <  2:856624/93‑1 (MQ=255)

                                                                          cttcttCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGa                  <  2:6612303/93‑1 (MQ=255)

                                                                            tcttcAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata                >  2:2108913/1‑93 (MQ=255)

                                                                            tcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata                <  2:437915/93‑1 (MQ=255)

                                                                            tcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata                <  1:2663058/93‑1 (MQ=255)

                                                                             cttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataa               <  2:132503/93‑1 (MQ=255)

                                                                              ttcAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataat              >  1:6495452/1‑93 (MQ=255)

                                                                              ttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataat              <  2:2772101/93‑1 (MQ=255)

                                                                              ttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataat              >  1:2435207/1‑93 (MQ=255)

                                                                                 aaTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATaaa           <  1:1180997/93‑1 (MQ=255)

                                                                                 aaTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATaaa           <  2:829330/93‑1 (MQ=255)

                                                                                  aTGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc          >  2:2557065/1‑93 (MQ=255)

                                                                                  aTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc          >  2:7295077/1‑93 (MQ=255)

                                                                                  aTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc          >  2:453242/1‑93 (MQ=255)

                                                                                   tGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACt         <  1:5554533/93‑1 (MQ=255)

                                                                                   tGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACt         <  1:988728/93‑1 (MQ=255)

                                                                                     gCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcg       <  1:4146592/93‑1 (MQ=255)

                                                                                      cGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgc      <  1:6935697/93‑1 (MQ=255)

                                                                                      cGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTcgc      <  2:7159234/93‑1 (MQ=255)

                                                                                          tGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTcc  >  1:5801598/1‑93 (MQ=255)

AAAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGCGTCC  >  NC_000913/3455202‑3455384

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