BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	54,157	C→T	P182P (CCG→CCA)	surA ←	survival protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	54,157	0	C	T	97.5%	548.7 / ‑0.1	160	P182P (CCG→CCA)	surA	survival protein
Reads supporting (aligned to +/- strand): ref base C (1/3); new base T (90/66); total (91/69)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 167 total aligned reads displayed.

GAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAG  >  NC_000913/54067‑54246

gAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTgg                                                                                         <  2:1759996/93‑1 (MQ=255)

  aaTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGa                                                                                       <  2:2560220/93‑1 (MQ=255)

   aTCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGAt                                                                                      <  1:4583430/93‑1 (MQ=255)

    tCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATc                                                                                     <  1:2427377/93‑1 (MQ=255)

    tCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATc                                                                                     <  2:4544967/93‑1 (MQ=255)

    tCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATc                                                                                     <  2:2875988/93‑1 (MQ=255)

      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGGTCAg                                                                                   >  1:1801236/1‑93 (MQ=255)

      aGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAg                                                                                   >  1:2791898/1‑93 (MQ=255)

         gcCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAt                                                                                >  2:4235250/1‑93 (MQ=255)

           cggtACGCGCCTGATCGACAATGGCGCTCTCCTGTCTTTCCGCTTCTTTCTCCTTATCAGAGGCCGGTTGTTCTTGCAGTGGGAGCAGGAtgg                                                                              >  2:4357189/4‑92 (MQ=255)

           cGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgt                                                                              >  1:1239647/1‑93 (MQ=255)

           cGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgt                                                                              >  1:940643/1‑93 (MQ=255)

           cGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGAtgt                                                                              <  1:2773383/93‑1 (MQ=255)

            gTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgtg                                                                             >  1:1744817/1‑93 (MQ=255)

            gTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgtg                                                                             <  2:1477942/93‑1 (MQ=255)

            gTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgtg                                                                             >  1:528850/1‑93 (MQ=255)

            gTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGAtgtg                                                                             <  1:109513/93‑1 (MQ=255)

                cgcgCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTc                                                                         >  1:3216033/1‑93 (MQ=255)

                cgcgCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTc                                                                         >  2:3540372/1‑93 (MQ=255)

                 gcgcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCa                                                                        >  2:276742/1‑93 (MQ=255)

                 gcgcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCa                                                                        <  1:845723/93‑1 (MQ=255)

                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAgg                                                                      >  1:4221559/1‑93 (MQ=255)

                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAgg                                                                      <  1:2293040/93‑1 (MQ=255)

                   gcCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGGTTTCCGGCAGTGGGATCAGGATGTGGCTCAgg                                                                      >  1:3690582/1‑93 (MQ=255)

                    ccTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGt                                                                     >  2:2966168/1‑93 (MQ=255)

                     cTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGtt                                                                    >  1:1378390/1‑93 (MQ=255)

                      tGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTc                                                                   >  1:4606679/1‑93 (MQ=255)

                      tGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTc                                                                   >  1:2699192/1‑93 (MQ=255)

                         tCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGc                                                                <  1:1963022/93‑1 (MQ=255)

                         tCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGc                                                                <  2:1343344/93‑1 (MQ=255)

                           gACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                              >  1:111990/1‑93 (MQ=255)

                           gACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTc                                                              >  1:2907263/1‑93 (MQ=255)

                            aCAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCa                                                             <  1:147473/93‑1 (MQ=255)

                             cAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAg                                                            <  2:1876727/93‑1 (MQ=255)

                              aaTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGt                                                           >  1:3748961/1‑93 (MQ=255)

                              aaTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGt                                                           >  1:2851288/1‑93 (MQ=255)

                              aaTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGt                                                           >  1:2810270/1‑93 (MQ=255)

                               aTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTTGGATCAGGATGTGGCTCAGGTTCAGCTCAGTg                                                          >  2:4325414/1‑93 (MQ=255)

                               aTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTg                                                          >  1:1784779/1‑93 (MQ=255)

                               aTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTg                                                          >  1:3836407/1‑93 (MQ=255)

                               aTGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTg                                                          >  2:1958228/1‑93 (MQ=255)

                                tGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         <  2:1130548/93‑1 (MQ=255)

                                aggCGCGAGCCTGGGTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGc                                                         <  2:4033595/92‑1 (MQ=255)

                                  gcgcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTg                                                       <  1:1075148/93‑1 (MQ=255)

                                   cgcccgccTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTgg                                                      <  2:939162/93‑1 (MQ=255)

                                    gcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGc                                                     <  1:2104392/93‑1 (MQ=255)

                                    gcgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGc                                                     >  2:2890726/1‑93 (MQ=255)

                                     cgcgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcg                                                    >  2:3157917/1‑93 (MQ=255)

                                       cgcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtc                                                  >  2:3417932/1‑93 (MQ=255)

                                        gcCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGcgtcg                                                 <  2:2010455/93‑1 (MQ=255)

                                          cTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGTATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               <  2:411433/93‑1 (MQ=255)

                                          cTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               >  1:81367/1‑93 (MQ=255)

                                          cTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               >  2:2103609/1‑93 (MQ=255)

                                          cTGGCTTTCCGCTTCGTGCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtt                                               <  1:302850/93‑1 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              <  1:120581/93‑1 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  1:3123440/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  2:1896795/1‑93 (MQ=255)

                                           tGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGttt                                              >  1:4094877/1‑93 (MQ=255)

                                            ggCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGtttt                                             >  2:3749573/1‑93 (MQ=255)

                                               tttCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGt                                          >  1:4513747/1‑93 (MQ=255)

                                               tttCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGt                                          >  2:2377222/1‑93 (MQ=255)

                                                ttCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGtt                                         <  2:1509245/93‑1 (MQ=255)

                                                  ccGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAc                                       <  2:152698/93‑1 (MQ=255)

                                                    gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAccc                                     >  2:637071/1‑93 (MQ=255)

                                                    gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAccc                                     >  1:3458973/1‑93 (MQ=255)

                                                    gCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTAccc                                     >  1:2221010/1‑93 (MQ=255)

                                                      ttCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACTCAc                                   >  2:3639699/1‑93 (MQ=255)

                                                      ttCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc                                   >  2:2970399/1‑93 (MQ=255)

                                                      ttCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCAc                                   >  1:584494/1‑93 (MQ=255)

                                                        cGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACct                                 <  1:941713/93‑1 (MQ=255)

                                                         gTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctg                                >  2:3447266/1‑93 (MQ=255)

                                                          ttCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgc                               <  2:3915733/93‑1 (MQ=255)

                                                          ttCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgc                               <  2:776029/93‑1 (MQ=255)

                                                          ttCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgc                               >  2:2903077/1‑93 (MQ=255)

                                                           tCACCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATAAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACctgct                              >  2:1913827/1‑93 (MQ=255)

                                                             aCCTGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCtgctgc                            <  1:781589/93‑1 (MQ=255)

                                                               cTGATCAGAGGGCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTgcgc                          <  2:2463637/93‑1 (MQ=255)

                                                                tGATCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGTTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGcc                         <  1:2639975/93‑1 (MQ=255)

                                                                  aTCAGAGGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAg                       <  1:3039324/93‑1 (MQ=255)

                                                                     agagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGATCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAgggc                    >  2:4630205/1‑92 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   <  2:1242142/93‑1 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  2:3902773/1‑93 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  1:1298084/1‑93 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  1:3515466/1‑93 (MQ=255)

                                                                      gagGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAt                   >  2:4130889/1‑93 (MQ=255)

                                                                       agGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  >  1:4462182/1‑93 (MQ=255)

                                                                       agGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  >  2:3471169/1‑93 (MQ=255)

                                                                       agGTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGAtt                  >  2:4179944/1‑93 (MQ=255)

                                                                        ggTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTc                 >  1:3301159/1‑93 (MQ=255)

                                                                        ggTCGGGTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGGTTTGGTTACCGCCCTGCGGCGCCAGGGAtgc                 >  1:982501/1‑91 (MQ=255)

                                                                            gggTTTTCCGGGAGTGGGATCAGGATGTGGCTCAGATTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACt             >  2:1585078/1‑93 (MQ=255)

                                                                            gggTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTAGACt             >  1:3719718/1‑93 (MQ=255)

                                                                              gTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTc           >  2:1985698/1‑93 (MQ=255)

                                                                              gTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTc           >  1:1868553/1‑93 (MQ=255)

                                                                              gTTTTCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTAGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTc           <  2:1684670/93‑1 (MQ=255)

                                                                               ttttCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTcc          <  1:4084011/93‑1 (MQ=255)

                                                                               ttttCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTcc          >  2:637232/1‑93 (MQ=255)

                                                                                  tCCGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGc       >  2:4709852/1‑93 (MQ=255)

                                                                                    cGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCgg     <  1:1355768/93‑1 (MQ=255)

                                                                                    cGGCAGTGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCgg     <  2:2067774/93‑1 (MQ=255)

GAAATCAGCGCCGTTACGCGCCTGATCGACAATGGCGCGCGCCTGGCTTTCCGCTTCGTTCACCTGATCAGAGGTCGGGTTTTCCGGCAGCGGGATCAGGATGTGGCTCAGGTTCAGCTCAGTGCTGGCGTCGTTTTGGTTACCCACCTGCTGCGCCAGGGATTCGACTTCCTGCGGCAG  >  NC_000913/54067‑54246

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