BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	52,554	A→G	L288P (CTG→CCG)	pdxA ←	PdxA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	52,554	0	A	G	99.5%	670.0 / ‑3.9	213	coding (863/990 nt)	pdxA	PdxA
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (107/105); minor base . (1/0); total (108/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 223 total aligned reads displayed.

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

aGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCggg                                                                                         >  1:729324/1‑93 (MQ=255)

  cgcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGcc                                                                                       <  2:2904691/93‑1 (MQ=255)

   gcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGccc                                                                                      >  1:4720110/1‑93 (MQ=255)

   gcCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGccc                                                                                      >  2:5103107/1‑93 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  1:1288363/1‑93 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   <  2:3410382/93‑1 (MQ=255)

      gTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGc                                                                                   >  1:27433/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  1:4419155/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  2:2476407/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  2:3218056/1‑93 (MQ=255)

       taataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCg                                                                                  >  2:2678749/1‑93 (MQ=255)

        aataaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGt                                                                                 >  1:469077/1‑93 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               >  1:2692505/1‑93 (MQ=255)

          taaAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTaa                                                                               >  1:3937841/1‑93 (MQ=255)

           aaaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAAt                                                                              <  1:4458746/93‑1 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             <  1:3631556/93‑1 (MQ=255)

            aaaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATg                                                                             >  2:570671/1‑93 (MQ=255)

             aaCTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGt                                                                            >  1:6250875/1‑93 (MQ=255)

                tGCCGCCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTca                                                                         <  1:784325/93‑1 (MQ=255)

                tGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTca                                                                         <  1:2302024/93‑1 (MQ=255)

                tGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGCGCCCAGCGTAATGTTca                                                                         <  1:5353355/93‑1 (MQ=255)

                 gCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcac                                                                        <  2:2989751/93‑1 (MQ=255)

                 gCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcac                                                                        >  2:6039621/1‑93 (MQ=255)

                  ccGACTTCCGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcaca                                                                       <  1:980049/93‑1 (MQ=255)

                  ccGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTcaca                                                                       <  1:1610993/93‑1 (MQ=255)

                   cGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGACTAAAGGGCGGGCCCAGCGCAATGTTCacac                                                                      >  1:3337518/1‑93 (MQ=255)

                   cGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCacac                                                                      <  2:3109642/93‑1 (MQ=255)

                    gACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACAcc                                                                     <  1:2444454/93‑1 (MQ=255)

                    gACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACAcc                                                                     >  2:595107/1‑93 (MQ=255)

                     aCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcg                                                                    >  1:81159/1‑93 (MQ=255)

                     aCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcg                                                                    >  2:2926588/1‑93 (MQ=255)

                     aCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcg                                                                    >  1:1135978/1‑93 (MQ=255)

                     aCATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcg                                                                    >  1:3117722/1‑93 (MQ=255)

                       aTCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACcgcg                                                                  <  1:3299772/93‑1 (MQ=255)

                         cGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGcc                                                                <  1:2720496/93‑1 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCg                                                              <  2:4205769/93‑1 (MQ=255)

                           gCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGC‑GGGCCAGCCGAAAGTTCAAACCGGGcccca                                                             >  2:1841580/1‑91 (MQ=255)

                             tttGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGaa                                                            <  1:5525259/93‑1 (MQ=255)

                                 cctcGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  1:5712249/90‑1 (MQ=255)

                                 ccACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGccc                                                        <  2:3455109/93‑1 (MQ=255)

                                   aCGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTg                                                      >  2:5121026/1‑93 (MQ=255)

                                    cGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTgg                                                     <  2:5186251/93‑1 (MQ=255)

                                    cGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTgg                                                     <  1:3737050/93‑1 (MQ=255)

                                    cGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTgg                                                     >  2:4414070/1‑93 (MQ=255)

                                      gcccGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTa                                                   <  1:2558191/92‑1 (MQ=255)

                                          gCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtttt                                               >  2:1297665/1‑93 (MQ=255)

                                          gCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTAtttt                                               >  1:5205988/1‑93 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCCCACCGGGCCCGAAGCCCTGGTATTTTa                                              <  1:4387243/93‑1 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              >  1:1513608/1‑93 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              >  2:86415/1‑93 (MQ=255)

                                           ccAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTa                                              >  1:4950550/1‑93 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            >  2:3863771/1‑93 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            >  1:5454273/1‑93 (MQ=255)

                                             aGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGc                                            >  2:912177/1‑93 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           <  2:1680400/93‑1 (MQ=255)

                                              gTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCa                                           <  2:4823482/93‑1 (MQ=255)

                                               ttCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCAc                                          >  2:5539519/1‑93 (MQ=255)

                                               ttCAAGCGCGGTGCCGTGGTCCACAGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCAc                                          <  1:6300205/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAGAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:1046928/93‑1 (MQ=255)

                                                    gcgcGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGaa                                     <  1:3733350/93‑1 (MQ=255)

                                                       cgGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAc                                  >  2:3535382/1‑93 (MQ=255)

                                                         gTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAccc                                <  2:4187528/93‑1 (MQ=255)

                                                         gTGCCGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGAccc                                <  2:3317621/93‑1 (MQ=255)

                                                             cGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGAt                            <  2:451106/93‑1 (MQ=255)

                                                             cGTGGTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGAt                            <  1:6071453/93‑1 (MQ=255)

                                                                ggTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGt                         <  2:1875138/93‑1 (MQ=255)

                                                                ggTCCACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGt                         <  1:2419597/93‑1 (MQ=255)

                                                                    cACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACTCCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTa                     >  1:3384372/1‑93 (MQ=255)

                                                                    cACTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTa                     >  1:4856519/1‑93 (MQ=255)

                                                                     aCTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTAc                    <  2:3967590/93‑1 (MQ=255)

                                                                      cTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACa                   >  2:3785576/1‑93 (MQ=255)

                                                                      cTGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACa                   <  2:560473/93‑1 (MQ=255)

                                                                       tGATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  >  1:297372/1‑93 (MQ=255)

                                                                       cgATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACAt                  <  1:1709146/92‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 <  1:5170904/93‑1 (MQ=255)

                                                                        gATGTGCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATc                 >  2:1921744/1‑93 (MQ=255)

                                                                          tgtgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGc               >  2:3044533/1‑93 (MQ=255)

                                                                           gtgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGcc              <  2:755606/93‑1 (MQ=255)

                                                                            tgCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCa             <  1:1388560/93‑1 (MQ=255)

                                                                             gCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAg            <  2:5387399/93‑1 (MQ=255)

                                                                             gCGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAg            <  2:5795229/93‑1 (MQ=255)

                                                                              cGAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGc           >  2:5267109/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  1:1300538/1‑93 (MQ=255)

                                                                               gAATAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCa          >  1:5585177/1‑93 (MQ=255)

                                                                                aaTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAc         <  1:5132758/93‑1 (MQ=255)

                                                                                aaTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAc         >  2:4610917/1‑93 (MQ=255)

                                                                                 aTAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcg        >  1:1796878/1‑93 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       >  2:209166/1‑93 (MQ=255)

                                                                                  tAAAGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcgg       <  2:5509032/93‑1 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      <  2:2245289/93‑1 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      >  2:634186/1‑93 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      <  1:1981080/93‑1 (MQ=255)

                                                                                   aaaGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggc      >  1:4040056/1‑93 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    <  2:825716/93‑1 (MQ=255)

                                                                                     aGGGCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgt    >  2:4728522/1‑93 (MQ=255)

                                                                                      gggcgggcCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtc   <  1:2104173/93‑1 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  <  1:5081750/93‑1 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  >  1:6225821/1‑93 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  <  2:4238839/93‑1 (MQ=255)

                                                                                       ggCGGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCAcggcgtcg  <  2:3742645/93‑1 (MQ=255)

AGCGCCGTAATAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGGTCCACTGATGTGCGAATAAAGGGCAGGCCCAGCGTAATGTTCACACCGCGCCCGAAGCCCTGGTATTTTAGCACGGGAAGACCCTGATCGTGGTACATCGCCAGCACGGCGTCG  >  NC_000913/52464‑52643

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