BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,621,256	G→A	E65K (GAA→AAA)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,621,256	0	G	A	87.2%	738.0 / 64.9	251	E65K (GAA→AAA)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand): ref base G (18/14); new base A (108/111); total (126/125)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 269 total aligned reads displayed.

AAAGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTAC  >  NC_000913/2621164‑2621348

  aGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                            <  2:4375698/93‑1 (MQ=255)

  aGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                            <  1:3244827/93‑1 (MQ=255)

  aGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaa                                                                                            <  1:6792081/93‑1 (MQ=255)

   ggATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaag                                                                                           >  1:3646107/1‑93 (MQ=255)

    gATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaaga                                                                                          >  1:82986/1‑93 (MQ=255)

     aTGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCaaagaa                                                                                         <  2:384856/93‑1 (MQ=255)

      tGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAAc                                                                                        >  2:5882945/1‑93 (MQ=255)

      tGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAAc                                                                                        <  2:5121857/93‑1 (MQ=255)

       gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACa                                                                                       <  2:6121019/93‑1 (MQ=255)

       gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                       <  2:1754084/93‑1 (MQ=255)

       gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                       >  1:2851590/1‑93 (MQ=255)

       gCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACa                                                                                       <  2:2946993/93‑1 (MQ=255)

        cGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACaa                                                                                      <  2:1051357/93‑1 (MQ=255)

        cGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACaa                                                                                      <  2:5352790/93‑1 (MQ=255)

         gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAATAACAAg                                                                                     <  2:4724300/93‑1 (MQ=255)

         gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                     <  2:3713997/93‑1 (MQ=255)

         gTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAg                                                                                     <  2:2037283/93‑1 (MQ=255)

          tttCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAgg                                                                                    <  2:115553/93‑1 (MQ=255)

           ttCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGc                                                                                   <  2:6559981/93‑1 (MQ=255)

           ttCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGc                                                                                   <  1:1538866/93‑1 (MQ=255)

             ccTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTc                                                                                 <  2:3627914/93‑1 (MQ=255)

                ggggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGGCGCATCATTATTAGCAAAGAACAAGGCTCCaa                                                                              <  1:7084989/93‑1 (MQ=255)

                 gggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAAc                                                                             >  2:2812903/1‑93 (MQ=255)

                 gggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAAc                                                                             <  1:1097587/93‑1 (MQ=255)

                  ggATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAAct                                                                            >  2:453669/1‑93 (MQ=255)

                    aTCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAActct                                                                          >  2:3356083/1‑93 (MQ=255)

                    aTCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAActct                                                                          >  1:3274888/1‑93 (MQ=255)

                     tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACtctc                                                                         >  1:7445944/1‑93 (MQ=255)

                     tCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACtctc                                                                         >  1:4747300/1‑93 (MQ=255)

                       tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCAt                                                                       >  2:1644755/1‑93 (MQ=255)

                       tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCAt                                                                       <  1:2422604/93‑1 (MQ=255)

                       tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                       >  2:1470559/1‑93 (MQ=255)

                       tATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCAt                                                                       <  1:6407776/93‑1 (MQ=255)

                          tCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTcc                                                                    >  2:4133114/1‑93 (MQ=255)

                          tCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTcc                                                                    <  2:1860536/93‑1 (MQ=255)

                           ccAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTccc                                                                   <  2:425025/93‑1 (MQ=255)

                           ccAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTccc                                                                   <  2:7090316/93‑1 (MQ=255)

                              aTAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGcc                                                                >  1:1230799/1‑93 (MQ=255)

                               tAACCTTGATGATTTCTATAAAGCCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGAAAAGAACAAGGCCCCAACTCTCATTCCCGCCa                                                               <  1:7328748/93‑1 (MQ=255)

                                 aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAAGTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                             <  1:6658599/93‑1 (MQ=255)

                                 aCCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCAtt                                                             <  1:5758987/93‑1 (MQ=255)

                                   cTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCGACTCTCATTCCCGCCATTTa                                                           >  2:7426045/1‑93 (MQ=255)

                                   cTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTa                                                           >  1:90957/1‑93 (MQ=255)

                                    ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                          <  1:5963937/93‑1 (MQ=255)

                                    ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                          >  1:5005965/1‑93 (MQ=255)

                                    ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                          <  2:3881642/93‑1 (MQ=255)

                                    ttGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTAc                                                          >  1:2980459/1‑93 (MQ=255)

                                     tgatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACt                                                         <  1:1932186/93‑1 (MQ=255)

                                      gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                        >  1:7160811/1‑93 (MQ=255)

                                      gatgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTg                                                        <  1:20706/93‑1 (MQ=255)

                                       atgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTgg                                                       <  1:995290/93‑1 (MQ=255)

                                       atgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTgg                                                       <  1:1568480/93‑1 (MQ=255)

                                        tgaGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTggg                                                      >  2:4647246/1‑93 (MQ=255)

                                        tgaGTTCTATAAAGACCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTggg                                                      >  2:484299/1‑93 (MQ=255)

                                            ttCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaaa                                                  >  2:2449554/1‑93 (MQ=255)

                                            ttCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCaaa                                                  >  2:799163/1‑93 (MQ=255)

                                              cTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAt                                                >  1:1353073/1‑93 (MQ=255)

                                              cTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAAt                                                >  2:6281309/1‑93 (MQ=255)

                                                ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                              >  2:7244668/1‑93 (MQ=255)

                                                ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                              >  1:2712333/1‑93 (MQ=255)

                                                ataAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTc                                              >  2:4056669/1‑93 (MQ=255)

                                                 taAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCa                                             <  2:5506957/93‑1 (MQ=255)

                                                 taAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCa                                             <  1:4391342/93‑1 (MQ=255)

                                                  aaaGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAg                                            >  2:7414163/1‑93 (MQ=255)

                                                  aaaGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAg                                            <  1:6552701/93‑1 (MQ=255)

                                                   aaGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGt                                           <  2:4079701/93‑1 (MQ=255)

                                                    aGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTc                                          <  1:18542/93‑1 (MQ=255)

                                                     gTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTcc                                         <  2:1263865/93‑1 (MQ=255)

                                                      tCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTccc                                        <  1:4731762/93‑1 (MQ=255)

                                                      tCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTccc                                        >  1:247845/1‑93 (MQ=255)

                                                      tCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTccc                                        <  1:2685778/93‑1 (MQ=255)

                                                         gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCggg                                     <  1:572292/93‑1 (MQ=255)

                                                         gCTTCGCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCggg                                     >  2:3177821/1‑93 (MQ=255)

                                                          cTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGt                                    >  1:4137278/1‑93 (MQ=255)

                                                              gCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTg                                >  1:7251982/1‑93 (MQ=255)

                                                              gCTGAACTGAAGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTg                                <  1:3800622/93‑1 (MQ=255)

                                                               ctgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGa                               <  1:5055581/93‑1 (MQ=255)

                                                                tgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaa                              >  2:7107928/1‑93 (MQ=255)

                                                                tgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaa                              >  2:6484282/1‑93 (MQ=255)

                                                                tgaactgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGaa                              >  2:2252785/1‑93 (MQ=255)

                                                                   actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                           >  1:4007881/1‑93 (MQ=255)

                                                                   actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                           <  1:5434211/93‑1 (MQ=255)

                                                                   actgaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGc                           >  1:756363/1‑93 (MQ=255)

                                                                     tgaaGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCg                         <  2:7473255/93‑1 (MQ=255)

                                                                      gaaGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGa                        <  1:4566390/93‑1 (MQ=255)

                                                                        aGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATc                      >  2:1215676/1‑93 (MQ=255)

                                                                        aGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATc                      >  2:6863474/1‑93 (MQ=255)

                                                                        aGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATc                      <  2:217382/93‑1 (MQ=255)

                                                                        aGCGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATc                      >  1:3551748/1‑93 (MQ=255)

                                                                          cGACGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAg                    <  2:2339152/93‑1 (MQ=255)

                                                                             cGCATCATTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGaa                 >  1:5509093/1‑93 (MQ=255)

                                                                                 tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg             >  2:3210046/1‑93 (MQ=255)

                                                                                 tcatTATTAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCg             >  1:417511/1‑93 (MQ=255)

                                                                                   attattAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGAc           <  2:5258853/93‑1 (MQ=255)

                                                                                   attattAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGAc           <  1:4646909/93‑1 (MQ=255)

                                                                                   attattAGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGAc           >  1:6337661/1‑93 (MQ=255)

                                                                                     tattaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACgg         >  1:7050143/1‑93 (MQ=255)

                                                                                      attaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc        >  1:215347/1‑93 (MQ=255)

                                                                                       ttaGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGcc       <  1:3910484/93‑1 (MQ=255)

                                                                                        taGCAAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCct      <  2:2573579/93‑1 (MQ=255)

AAAGGATGCGTTTCCTGGGGATCTATTCCAATAACCTTGATGAGTTCTATAAAGTCCGCTTCGCTGAACTGAAGCGACGCATCATTATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTAC  >  NC_000913/2621164‑2621348

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