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	99.6%	902.6 / ‑5.0	265	W427R (TGG→CGG)	b3325	YheF
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (119/145); total (119/146)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

Only 100 of 280 total aligned reads displayed.

AAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGC  >  NC_000913/3455203‑3455380

aaaaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCggg                                                                                       <  1:4981068/93‑1 (MQ=255)

 aaaGAAGTGTGGGATCACCCGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCgggg                                                                                      <  2:2838230/93‑1 (MQ=255)

  aaGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCggggg                                                                                     >  1:5318728/1‑92 (MQ=255)

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

   aGAATGGTGGGATCACCAGCGCGTATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAg                                                                                    <  2:4738941/93‑1 (MQ=255)

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

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

       tggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACt                                                                                >  1:1354915/1‑93 (MQ=255)

       tggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACt                                                                                >  2:7915867/1‑93 (MQ=255)

       tggtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACt                                                                                >  2:3395041/1‑93 (MQ=255)

         gtggGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGc                                                                              <  1:3663452/93‑1 (MQ=255)

              gtCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAcc                                                                         <  1:2970619/92‑1 (MQ=255)

              aTCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTAcc                                                                         >  1:6644915/1‑93 (MQ=255)

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

                cACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgc                                                                       >  1:7118177/1‑93 (MQ=255)

                cACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgc                                                                       >  2:7524852/1‑93 (MQ=255)

                 aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgcg                                                                      >  2:8375701/1‑93 (MQ=255)

                 aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgcg                                                                      >  1:4623745/1‑93 (MQ=255)

                 aCCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACcgcc                                                                      >  1:574623/1‑92 (MQ=255)

                  ccAGCGCGGCTCCTTCCTGGGATATGTTTCGCGCCTGCAGTGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCgcgc                                                                     <  1:6701708/93‑1 (MQ=255)

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

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

                    cgcgcGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTg                                                                   <  1:2673659/92‑1 (MQ=255)

                    aGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTg                                                                   >  1:8462534/1‑93 (MQ=255)

                    aGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTg                                                                   >  2:5741149/1‑93 (MQ=255)

                     gcgcGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTgg                                                                  <  2:3084481/93‑1 (MQ=255)

                      cgcgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGc                                                                 <  1:6572122/93‑1 (MQ=255)

                        cgAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCa                                                               >  2:4451362/1‑93 (MQ=255)

                          aaTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGcccga                                                             >  1:2051397/1‑90 (MQ=255)

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

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

                           aTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTa                                                            <  2:6951088/93‑1 (MQ=255)

                           aTCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTa                                                            <  1:3871428/93‑1 (MQ=255)

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

                             cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGCCCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                          <  1:8522873/93‑1 (MQ=255)

                             cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                          >  2:6540868/1‑93 (MQ=255)

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

                             cctgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAAc                                                          >  1:7258128/1‑93 (MQ=255)

                               tgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACaa                                                        <  1:6733869/93‑1 (MQ=255)

                               tgcctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACaa                                                        <  2:5031079/93‑1 (MQ=255)

                                 cctgGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATa                                                      <  2:6232970/93‑1 (MQ=255)

                                    gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                   <  2:1195321/93‑1 (MQ=255)

                                    gggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaa                                                   <  2:5851173/93‑1 (MQ=255)

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

                                     ggATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATaaaaa                                                  <  1:8681588/93‑1 (MQ=255)

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

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

                                       atatGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATg                                                <  2:2857127/93‑1 (MQ=255)

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

                                          tGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAAtacca                                             >  1:3479776/1‑89 (MQ=255)

                                          tGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACa                                             >  1:5455660/1‑93 (MQ=255)

                                          tGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACa                                             >  1:2787986/1‑93 (MQ=255)

                                           gTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACAt                                            >  2:7795968/1‑93 (MQ=255)

                                            tttAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATc                                           >  1:5093669/1‑93 (MQ=255)

                                            tttAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATc                                           <  2:329704/93‑1 (MQ=255)

                                             ttAGCGTCTACAACGGCATGGCCGCAGGCTTCTTCAATGGCGCCCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATcc                                          <  1:6975145/93‑1 (MQ=255)

                                             ttAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATcc                                          <  2:8104217/93‑1 (MQ=255)

                                              tAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCt                                         >  2:7487386/1‑93 (MQ=255)

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

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

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

                                                  gcCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGcc                                     >  1:6302546/1‑93 (MQ=255)

                                                   ccTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGcga                                    >  2:531850/1‑91 (MQ=255)

                                                   ccTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCa                                    <  2:2036914/93‑1 (MQ=255)

                                                    cTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAc                                   <  2:333166/93‑1 (MQ=255)

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

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

                                                       cAATGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCAcccc                                <  2:3970974/93‑1 (MQ=255)

                                                        aaTGGCATGGCCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCa                               >  1:4474633/1‑93 (MQ=255)

                                                         atggcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCaa                              >  2:7211782/1‑93 (MQ=255)

                                                         atggcatggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCaa                              >  1:1864855/1‑93 (MQ=255)

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

                                                         atggcatggcCGCAGGCTTCTTCAATGGAGACCGGGGAGTACTGCTTACCGAGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCaa                              >  2:5714216/1‑93 (MQ=255)

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

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

                                                             catggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCAt                          >  2:1627187/1‑93 (MQ=255)

                                                              atggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATc                         <  1:7508063/93‑1 (MQ=255)

                                                              atggcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATc                         <  2:8308516/93‑1 (MQ=255)

                                                                 gcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCCAGCATCGTa                      >  2:3203088/1‑93 (MQ=255)

                                                                 gcCGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTa                      <  1:7218956/93‑1 (MQ=255)

                                                                  ccGCAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTaa                     <  1:8516321/93‑1 (MQ=255)

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

                                                                   cGCAGGCTTCTTCAATGGCGACCGGGGAGATCTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAAc                    >  1:7687903/1‑93 (MQ=255)

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

                                                                     cAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                  <  2:3726840/93‑1 (MQ=255)

                                                                     cAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                  <  2:621375/93‑1 (MQ=255)

                                                                     cAGGCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGc                  >  1:1155924/1‑93 (MQ=255)

                                                                      aGGCTTCTTCAATGGCGACCGGGGAGAACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTCACGCt                 >  1:3371612/1‑93 (MQ=255)

                                                                       ggCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTg                >  2:7548563/1‑93 (MQ=255)

                                                                       ggCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTg                >  1:7871380/1‑93 (MQ=255)

                                                                        gCTTCTTCAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTgg               >  1:8246022/1‑93 (MQ=255)

                                                                          ttcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGat             <  1:3118483/93‑1 (MQ=255)

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

                                                                           tcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata            >  1:5827627/1‑93 (MQ=255)

                                                                           tcttcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGata            >  1:6505619/1‑93 (MQ=255)

                                                                              tcAATGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGataata         >  1:4444465/1‑93 (MQ=255)

                                                                                aaTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATaaa       >  1:5220760/1‑93 (MQ=255)

                                                                                 aTGGCGACCGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAAc      <  1:1580556/93‑1 (MQ=255)

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

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

AAAAGAATGGTGGGATCACCAGCGCGAATCCTGCCTGGGATATGTTTAGCGCCTACAATGGCATGGCCGCAGGCTTCTTCAATGGCGACTGGGGAGTACTGCTTACCGCGCTGGCCAGTAACAATAAAAATGACATCCTCGCCACCCCAAGCATCGTAACGCTGGATAATAAACTCGC  >  NC_000913/3455203‑3455380

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